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The x-ray lines La, and Laz have been examined for ac- 
companying satellites in the atomic number range from 
Ta(73) to U(92). La is found to have two satellites, herein 
called La‘** and La*. The former (longer wave-length) ex- 
tends from Au(79) to U(92); the latter from Os(76) to 
Bi(83). Even over the short range of atomic numbers over 
which they occur, these two satellites change rapidly as to 
relative intensity and appearance. They do not seem to be 
connected with the prominent group of satellites of La 


which occur in the region 36<Z<54. Lf, is found to have 
two satellites, Lf,’ and Lf,"’. LB’’ extends from Ta(73) to 
U(92). LB,’ begins several atomic numbers below Ta(73) 
and extends to U(92). The satellite structure accompanying 
Lf, in this high atomic number range is somewhat similar 
to that found in the range 40<Z<54. But the absence of 
any satellites of Lf, in the range 54<Z<70 suggests that 
these high atomic number satellites may be new lines. 


NE of the most significant characteristics of 
x-ray satellites is the fact their atomic 
number range is not at all coincident with that of 
the parent lines with which the satellites are 
respectively associated. For example, although 
Le is first observed at Va(23), the satellites of La 
are not found until Ni(28)' or Cu(29)? is reached. 
For higher atomic numbers the satellites become 
more prominent, reaching a maximum of in- 
tensity, relative to the parent line for elements in 
the neighborhood of Rh(45). Beyond Rh(45) the 
satellite structure becomes less prominent. At 
Sn(50) only a broad, unresolved band remains.* 
This band is very faint at Te(52). No trace of 
satellite structure accompanying La is observable 
for Ba(56). 


* The senior author wishes to express his best thanks to 
the Heckscher Research Council of Cornell University for 
a grant in aid of this investigation. 

i A. Karlsson, Ark. Math. Astr. o. Fys. (Upsala) 22, 

30. 

* J. H. Van der Tuuk, Dissertation, Groningen, 1928. 

*F. K. Richtmyer, J. Franklin Inst. 208, 325 (1929). 


Similarly it has been shown‘ that the promi- 
nent satellites of LS. are not observed for the 
elements immediately above I(53). 

There have been reported in the literature, 
however, numerous satellites of LZ lines in the 
spectra of the elements of higher atomic number, 
particularly from Ta(73) to U(92). Since there 
has been no systematic and searching study of 
satellites in this atomic number range, the 
present work was undertaken primarily to ascer- 
tain whether these high-atomic-number satellites 
are in fact, as has been assumed,° identical with 
the satellites found associated with elements of 
lower atomic number; or whether they are to be 
regarded as new lines. 

The spectrometer used was of the Siegbahn- 
Thoraeus high-vacuum type,® the distance be- 


*R. D. Richtmyer, Phys. Rev. 38, 1802 (1931). 

5 These data are most readily found in tables in Sieg- 
bahn's Spektroskopie der Réntgenstrahlen, second edition, 
(1931), Siegbahn uses the term “spark lines” instead of 
“satellites.” 

* J. Opt. Soc. Am. and Rev. Sci. Inst. 13, 235 (1926). 
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tween the entrance slit and the photographic 
plate being 54.9 cm for some of the plates; and 
211.4 cm for others. The higher resolution 
attainable at the greater distance, with suitable 
slit widths, is an important factor in studying the 
satellite structure. 


SATELLITES OF La 


In the atomic number range 29<Z<50, La 
has five (in some elements more) satellites, which, 
beginning with In(49), decrease rapidly in in- 
tensity and are not observed above Xe(54). 
Contrary to the observations of Coster? and 
Dauvillier® no trace of satellites accompanying 
La has been found by the senior author for 
Ba(56), although spectrum plates have been 
made and examined with great care. If satellites 
do accompany Ba La, they must be exceedingly 
faint. From La(57) to Yb(70) inclusive, no 
observer has reported satellites of La. Dauvillier® 
reports one satellite (La;, Siegbahn’s notation) 
for Lu(71), although no such line was found by 
Wennerléf" or by Coster.” Two observers, namely 
Wennerléf'® and Auger and Dauvillier™ report a 
satellite (Las, Siegbahn) for Ta(73), in dis- 
agreement with Idei,!* whose work appears to 
have been most carefully done. Of the 14 
observers quoted by Siegbahn as having meas- 
ured La; of W(74) only Rogers® reports a 
satellite for La. The present authors have made 
many plates of Ta(73) and of W(74) under a 
variety of different conditions, and have failed 
to find the slightest trace of a satellite of La for 
either element. One concludes, therefore, either 
that no satellites exist for La between Cs(55) and 
W(74) inclusive, or, if there is a satellite in this 
atomic number range, it is so very faint as to 
have escaped detection by most of the observers. 

For atomic numbers higher than W(74) there 
is ample evidence that La is accompanied by one 
(or more) satellites, except possibly in the case of 
Rh(75), for which Beuthe" alone reports La;. 


7 Coster, Phil. Mag. 43, 1070; 44, 545 (1922). 

_ * Dauvillier, C. R. 174, 1347 (1922). 
® Dauvillier, J. d. Physique et le Radium 3, 230 (1922). 
” Wennerléf, Ark. Math. Astr. 0. Fys. 22, No. 8 (1930). 
" Auger and Dauvillier, C. R. 176, 1297 (1923). 
12 Tdei, Tohoku Imp. Univ., Sci. Rep. 19, 651 (1930). 
8 Rogers, Proc. Camb. Phil. Soc. 21, 430 (1923). 
4 Beuthe, Zeits. f. Physik 46, 873 (1928). 


In the present work, in addition to the 
observations on Ta(73) and W(74) mentioned 
above, we have studied the satellite structure 
accompanying La for the elements Os(76), Ir(77), 
Pt(78), Au(79), T1(81), Pb(82), Bi(83), Th(90) 
and U(92). Two satellites were found which we 
shall call La* and La*.t La* is found in the 
atomic number range Os(76)"4 to Pb(83); La‘ in 
the range of Au(79) to U(82), although we failed 
to find definite evidence of this line for T1(81), 
since on our plates the Lf, line of W(74), 
appearing because of our use of a tungsten 
filament in the x-ray tube, comes at exactly the 
same wave-length as would La* of T1(81). Data 
regarding these two satellites are given in Table I. 


TABLE I. Satellites of La for the range Os(76) to U(93). 


Lo® 
dv dv 

Element|A(X.U.)* (a,;—a'*) R (a,;—a*) R R 
Os 76 1383.07 5.52 2.62 1.62 
Ir 77 1343.17 5.30 2.66 1.63 
Pt 78 1305.13 5.20 2.77 1.66 
Au 79 | 1270.54 3.23 1.82 1.35 | 1268.63 5.14 2.89 1.70 
Tl 81 ** 1200.00 4.93 3.10 1.76 
Pb 82 | 1169.29 3.26 2.16 1.47 | 1167.85 4.73 3.13 1.77 
Bi 83 | 1138.14 3.36 2.35 1.53 | 1136.95 4.55 3.18 1.78 
Th 90} 950.80 3.25 3.26 1.81 
U 92 | 905.35 3.39 3.72 1.93 


*Siegbahn’s values of \ for La, for each element are 
taken as standard. 
** Obscured by Lf, of W(74). 


The values of (Av/R)! are shown graphically in 
Fig. 1 as a function of atomic number for these 
satellites, La** and La*. The circles represent the 
data of the present authors. The crosses are from 
the observations of Idei.'*? For comparison, there 
are plotted on the same graph the (six) satellites 
of La in the atomic number range 36<Z<50, 
from the measurements of Richtmyer and 
Richtmyer.'® 

If one relies solely upon the graphical relations 
shown in Fig. 1, one might conclude that the 


t Idei (reference 12) recognizes the existence of these two 
lines, although he did not find both of them for any one 
element. He designates them La’ and La’’, respectively. 
Since it is obvious from Fig. 1 that these lines are not to be 
identified with the satellites Le’ and La” which occur in the 
atomic number range Rb(37)—Te(52), we have adopted the 
symbols La‘ and Le* to avoid confusion. Siegbahn 
(reference 5) does not distinguish between these two lines 
and calls them Lag. 

18 Richtmyer and Richtmyer, Phys. Rev. 34, 574 (1929). 
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Fic. 1. (Av/R)! as a function of atomic number for satellites of La. (Av =difference 
in frequency between satellite and parent line.) The crosses are from the measure- 
ments of Idei.'* The circles, for elements 76-92, are from data by present authors; for 
elements 37-50, from data by Richtmyer and Richtmyer.* 


satellite designated La’ is to be identified with 
La” and La’”’, perhaps unresolved for these high 
atomic number elements, which latter are the 
most intense satellites of the lower atomic 
number range. 

There are, however, weighty arguments against 


this conclusion. First of all, there is no evidence 
that additional satellites accompany La‘ and 
La?*, similar to La‘* and La’ which are associated 
with La’’ and La’” in the lower range. The 
intensities of La‘® and La’ are such that, if 
present with La’, they should be readily de- 
tected. Second, the combination La* and La? 
presents a structure very different from La’ 
—La’', as may be seen from Fig. 1. Third, and 
most important, La* appears very abruptly at 
Os(76), and can therefore hardly be regarded as 
an extension of La’’+’’’, which, as above noted, 
are not found above I(53). 

Of greater importance than the wave-length 
and atomic number range of La“ and La’, is the 
rapid change in character of these lines as one 
advances up the atomic number scale above 
Os(76). At Os(76) La? is a sharp, fairly prominent 
line, well separated from La;. At Ir(77) the 
satellite is still sharp, but the “‘valley’’ between 
it and the parent line is less pronounced. At 
Pt(78) this ‘“‘valley’’ has become almost as 
intense as the satellite itself, which now appears 
almost like a sharp, prominent edge to a narrow 
band. At Au(79) there is unmistakable evidence 
of the presence of ‘wo lines, the new line, La®, 


appearing in what was formerly the valley 
between La? and La. Very likely the “filling up” 
of this valley in Ir(77) and Pt(78) foreshadows 
the appearance of La‘. For T1(81), Pb(82) and 
Bi(83) there is clear evidence of the presence of 
both La* and La’, although, as above noted, 
La* for T1(81) is probably covered up by Ly of 
W(74). For Th(90) and U(92) our plates are not 
as intense as they should be for best observation, 
but unmistakably there is only one satellite 
present, namely La‘. 

One concludes tentatively, therefore, that La* 
and La? are new satellites, having a different 
origin from the group of satellites La’—La”. 
Furthermore, the rapid change in the character 
of these two satellites with atomic number after 
passing Os(76) suggests that their origin is in 
some way connected with the changes which take 
place in the Oyy and the several P electron shells 
in the region of Os(76) to U(92), there being over 
this atomic number range no change in the L and 
M shells in which Lay, originates. 

It is not easy to reconcile these data with the 
Wentzel-Druyvesteyn theory, according to which 
satellites have their origin in atoms which are 
doubly ionized as to their inner electron shells. 
W(74) differs from Os(76) only as regards the 
number of electrons in the Ory shell: W(74) has no 
electrons in this shell while Os(76) has two. Yet in 
Os(76) there is a prominent satellite while in 
W(74) there is none. 
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Fic, 2. (Av/R)! as a function of atomic number for satellites of Lz. The crosses 
are from the measurements of Idei."* The circles, for elements 73-92, are from data 
by present authors; for elements 40-54, from data by Richtmyer and Richtmyer. 


SATELLITES FOR 


In the atomic number range 40<Z<54, LBy 
has four satellites, well separated from the 
parent line. Two of these, L82’ and LB,” are quite 
sharp and, for satellites, are relatively intense. 
Associated with these two lines is a_ short 
continuous spectrum extending toward shorter 
wave-lengths from well past Midway 
between L8,’ and LB.” is a faint, diffuse satellite!® 
called L62(c); beyond Lf.” is a fourth line, 
similarly faint and diffuse, called'® LB2(d). This 
type of structure persists substantially un- 
changed, except as to intensity relative to parent 
line, over the above atomic number range. 

These satellites disappear, abruptly, at 1(53), 
as has been shown by R. D. Richtmyer,! who 
pointed out that the lines designated by Coster'® 
as satellites, in the spectrum of Ba(56) are, in 
reality, diagram lines (L8; and L8,o respectively). 
There are no other satellites of L82 reported in the 
literature for elements above Ba(56) until Er(68) is 
reached. For this element Coster’? reports a 
single satellite, L8:’. This one line is reported by 
various observers for the elements Yb(70) to 
Ir(77). For the remaining elements two satellites 
LB’ and Lf,” are reported. 

Our own measurements begin at Ta(73) and 
include all elements readily available up to 
U(92). For each element studied, two satellites, 
LB,’ and LB,” have been found,f as is shown in 
Table IT. 
~ 18 Coster, Phil. Mag. 44, 545 (1922). 


17 Coster, Phil. Mag. 43, 1070 (1922). 
t To save confusion in terminology Siegbahn's designa- 


(Av/R)* as a function of Z is shown in Fig. 2, in 
which, for comparison there is also shown the 
four satellites of LB2 in the lower atomic number 
range. The circles are the observations of the 
present authors. The crosses are from the work of 
Idei." 


TABLE II. Satellites of LB, for the range Ta(73) to U(92). 


Lp,’ 


Av ‘Av\4 Av 
Element/A(X.U.)* (82—82’) R R A(X.U.) (B2—-B2"") R R 


3.46 1.86 | 1273.57 8.33 4.62 2.15 
3.43 1.85 | 1233.53 8.50 5.02 2.24 
3.65 1.91 | 1161.13) 7.71 5.15 2.27 
36 «1.95 | 1125.85 7.12 5.07 
3.86 
4.12 
4.17 


Ir(77) | 1127.61 5 2.25 
Pt(78) | 1094.64 5.10 1.96 | 1092.93 681 5.15 2.27 
Au(79) | 1062.86 5.15 2.03 | 1061.28 6.73 5.39 2.32 
TI(81) | 1003.58 464 4. 2.04 | 1002.07 6.15 5.53 2.35 
Pb(82) | 976.09 974.50 

Bi(83) | (masked by 181) 947.62 5.62 5.67 2.38 
Th(90) | 788.47 3.45 5.02 2.24 | 787.36 4.56 6.66 2.58 
U(92) 749.92 3.15 5.07 2.25 | 748.93 4.14 6.66 2.58 


*Siegbahn’s values for LS, for each element taken as 
standard. 


As in the case of the satellites of La;, the wave- 
length data, either in tabular or graphical form, 
do not tell the whole story. There is no doubt 
concerning the absence of both 82’ and 82” over a 
considerable range of atomic numbers in the rare 
earth region. There is agreement in the literature 
that 2’ is the first to reappear. It seems to 
increase in intensity until Os(76) is reached, 


tion of these high atomic number satellites of LA, namely 


LB,’ and Lf,"’, has been retained, although it is by no 
means certain that these lines are to be identified with 
those similarly designated in the lower atomic numbet 
range. 
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beyond which its intensity seems to be constant. 
Our plates show clearly that, although 8,’ is 
unmistakably present at Ta(73) it is much 
weaker than £2’, and it grows in intensity less 
rapidly than does 82’. Not until T1(81) or Pb(82) 
is reached do the two lines have the same 
intensity, which equality is maintained up to 
U(92). 

The behavior of these two satellites in the 
high-atomic number range is thus, as was found 
to be the case with the satellites of La, quite 
different from that observed in the lower range in 
which the two lines maintained the same 
intensity relative to each other throughout. 

A further point of difference is that while LS.’ 
and L62"" are equally sharp for the range 
410<Z<54, LB,” is noticeably “fuzzy” for the 
elements Ta(73) to Pt(78). 

Again, as in the case of La, there is no change 
in the electron configurations in the shells which 
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give rise to LB, over the atomic number range for 
which the satellites are observed. 

The following two conclusions therefore, seem 
to be justified by the data above reported: First, 
the high atomic number satellites of La and LB. 
are not to be identified with the satellites which 
accompany these lines at the lower atomic 
number range. Second, the peculiar behavior of 
these satellites in passing from element to 
element in this high atomic number region, seems 
to indicate a close connection of some kind 
between the satellites and the peripheral electron 
structures of the atom. 

Since it is the purpose of this paper to present 
experimental results, no further comments will be 
made concerning the bearing of these data on 
current theories of the origin of satellites, beyond 
the obvious remark that any acceptable theory 
must be in agreement with the observations 
herein reported! 
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Soft X-Rays from (100) and (111) Faces of Copper Single Crystals 


MorcGan L. Brown University 
(Received July 13, 1933) 


Between 90 and 220 volts 14 critical potentials were 
found for both the (111) and the (100) faces. Each break is 
the same within 1.5 volts for the two targets. Agreement 
with published results for polycrystalline copper is fair, 
although there are fewer discontinuities than for the 
ordinary metal, as in other experiments with single 
crystals. The critical potentials may be grouped to show 
level differences, attributed by Richardson to structure 
electrons in the crystal lattice. Comparison with other 
results for metals of similar structure casts doubt on this 
theory. Kronig’s theoretical interpretation of fine structure 
in x-ray absorption presents a possible explanation of soft 


x-ray critical potentials. Assuming 4 to 7 volts inner 
potential, 10 of the 14 observed critical potentials agree 
with Kronig’s calculated energy discontinuities and only 
one calculated discontinuity was not observed. A strong 
break at 104 volts, not predicted by this theory, is within 2 
volts of Kurth’s value for the N series convergence 
frequency, Similar good agreement exists between Thomas’ 
critical potentials for polycrystalline iron and nickel and 
the energy discontinuities calculated from Kronig’s theory. 
From these results it appears that the theory of forbidden 
energy zones may account for soft x-ray critical potentials. 


INTRODUCTION 


OFT x-ray emission measurements by critical 
potential methods show many discontinuities 

not explained by the simple Bohr theory. 
Andrewes, Davies and Horton' suggested that 
some of these critical potentials might be due to 
the arrangement of atoms at the target surface, 
since the M electrons are not very deeply seated 
within the atoms of lighter elements. This sug- 
gestion was supported by Farnsworth’'s? obser- 
vation that the position and shape of many of the 
discontinuities in secondary electron emission 
curves depend on the surface arrangement of 
atoms in the target or the orientation of single 
crystal targets. Using an imperfect natural 
graphite crystal, Richardson and Andrewes* 
found a decided reduction in the number of 
discontinuities, compared with those previously 
found for carbon targets. Richardson and Rao* 
reported similar results for the (100) face of a 
nickel single crystal. Using the same (100) face 
of nickel Rao® found the changes in slope of 
secondary electron curves to be in good agree- 


1 Andrewes, Davies and Horton, Proc. Roy. Soc. A117, 
649 (1928). 

? Farnsworth, Phys. Rev. 25, 41 (1925); 31, 419 (1928). 

* Richardson and Andrewes, Proc. Roy. Soc. A128, 1 
(1930). 

* Richardson and Rao, Proc. Roy. Soc. A128, 57 (1930). 

5 Rao, Proc. Roy. Soc. Al28, 57 (1930). 


ment with the positions of the soft x-ray dis- 
continuities. With a (111) face of a _ nickel 
crystal, Rupp® found electron diffraction maxima 
and soft x-ray breaks to agree as either the 
primary voltage or the angle of incidence was 
varied. He concluded that some, at least, of the 
soft x-ray critical potentials are connected with 
properties of the crystal grating. 

The object of the present investigation was to 
determine whether there is a real difference in the 
soft x-ray critical potentials for different faces of 
the same metal crystal. 


DESCRIPTION OF APPARATUS 


Fig. 1 shows the electrical connections and a 
cross section of the experimental tube which is of 
Pyrex in the form of a cross consisting of a main 
tube G, 7 cm in diameter, with side arms 
A-—A of 3 cm diameter placed near the top of the 
main tube. Targets 7, and 72, mounted on a pair 
of quartz rods R, may be moved by magnets 
acting on —M. The targets were both cut from 
the same single crystal made in this laboratory.’ 
The faces were prepared as described by Farns- 
worth. Ammonium persulphate solution was 
used for final etching of the (100) target. The 


® Rupp, Naturwiss. 18, 880 (1930). 
7 Farnsworth, Phys. Rev. 40, 690 (1932). 
8 Farnsworth, Phys. Rev. 34, 679 (1929). 
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Amplifier 275 volts volts 
F022 


Fic. 1. Experimental tube and electrical connections. 


(111) target was etched electrolytically in slightly 
acid copper sulphate solution. Low current 
density of about 3 m.a. per cm? produced best 
development of the (111) face. Final dimensions 
of the targets were about 8.07.5 mm face 
and 10.5 mm height. A cylindrical nickel shield C, 
24 mm in diameter, surrounded the targets and 
with collar D supported the rods carrying them. 
An opening //7, 18 X19 mm, permitted exposure of 
only one target at a time to the filament F. Both 
targets were connected to the tungsten seal L 
by a tungsten spring £. A filament B of 8 mil 
tungsten was placed so that the targets could be 
withdrawn to position 7,’ and heated by electron 
bombardment of their molybdenum’ supports 
which fitted into grooves cut in the targets. The 
x-ray filament F of oxide coated nickel, placed 7 
mm below the target face, was usually operated 
at 2.2 amperes, so that the drop across the 
emitting portion was less than 0.5 volt. A large- 
capacity storage battery supplied the heating 
current, which was controlled by a tapped oil 
immersed chromel resistance W and clips sliding 
on nickel wires for fine adjustment. The ion trap 
consisted of two sets of screens and vanes S,V; 
and S:Ve2, of 1/16-inch mesh nickel screening, 
insulated by Pyrex tubes. The potentials, —9.2 
volts for the top screen S,V, and cylinder C, and 
+275 volts for the lower screen S2V2, were 
supplied by dry cells shielded in sheet iron boxes. 
The photoelectric plate P of sheet nickel rested 
on the bottom of the main tube G and was 
insulated from the screens by 3 mm rods Q of 
clear fused quartz. The insulation resistance of P 
was higher than 10'° ohms. The photoelectric 
current was measured by a one-tube balanced 
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amplifier of the type described by Soller,’ 
slightly modified, having a sensitivity of 7.33 
amp./mm with grid resistance 3.74 10° 
ohms (pencil mark on amber). Part of the plate 
resistance R, (Soller’s figure) was made variable, 
so the plate current could be balanced out for 
any grid input in the operating range. The effect 
was the same as in Compton and Thomas’ 
method" of balancing out part of the electrome- 
ter current, allowing use of the full sensitivity to 
measure changes of photoelectric current over a 
small range of target voltages. The variable 
resistance used in this compensating circuit was a 
400 ohm Centralab potentiometer. Large B 
batteries supplied the x-ray target voltage which 
could be adjusted to any value by taps and 
potentiometers. It was measured by a special 
high-resistance Weston voltmeter, which could 
be read accurately to less than 0.1 volt on the 
300 volt range. The thermionic current J;, of the 
order 0.5 m.a., was measured to 1 part in 5000 by 
a null circuit using a galvanometer to indicate 
balance of a standard cell against the potential 
drop in a resistance carrying J,. Balance was 
secured and J; held constant by changing the 
filament current. The experimental tube was 
shielded with Aquadag, and all wires entering 
the room were fully shielded. 

The vacuum system, consisting of two mercury 
vapor pumps in series, backed by a Hyvac pump, 
was connected to the experimental tube through 
two traps immersed in carbon dioxide snow and 
acetone. Liquid air was used until experiments on 
contact potential difference in this laboratory 
had shown that results were the same with either 
carbon dioxide or liquid air. Pressure was 
measured by a General Electric F P62 ionization 
gauge placed near the x-ray chamber. The 
experimental tube and trap nearest to it were 
outgassed by baking for 3 periods totalling 50 
hours at 430 to 480°C. Between bakings the 
filaments were flashed and the targets heated 
several hours at red heat. After this outgassing 
process, pressures while bombarding the targets 
did not rise above 10~*° mm of mercury. Pressures 
while taking readings, 2 to 5 hours after heating 
the targets, were 1 to 210-7? mm. 


® Soller, Rev. Sci. Inst. 3, 416 (1932). 
10 Compton and Thomas, Phys. Rev. 28, 601 (1926). 


|_| 
A 
= 
i an 

e 
f | 
a 
n 
e 
r 
n 


MORGAN L. 


EXPERIMENTAL PROCEDURE 


Before each run, the target was heated 2 to 10 
hours. As outgassing progressed, the target 
temperature was gradually increased but evapo- 
ration was prevented until the last stages. At the 
conclusion of bombardment the target was 
immediately moved to the x-ray position and the 
target voltage adjusted at the highest value to 
be used in the run. As soon as the thermionic 
current J; could be held at the desired value, the 
amplifier was adjusted by means of the variable 
resistance in the plate compensating circuit, so 
that its reading 7,’ would be near the higher end 
of the scale. With the target voltage and /, 
constant, readings of J,’ were taken at intervals 
until it ceased to change. During this time, 2 to 5 
hours, J,’ decreased, rapidly at first, then more 
slowly. This is attributed to changing photo- 
electric characteristics of the nickel plate as it 
was outgassed by action of the soft x-rays. 
Nakaya" observed similar effects on outgassing a 
photoelectric plate by bombardment. 

When /,’ had reached a steady state, readings 
were taken from which the soft x-ray critical 
potential curves were plotted. The amplifier 
sensitivity was adjusted by an Ayrton shunt on 
the galvanometer so that J,’ would cover nearly 
the whole scale for the range of voltages to be 
used in the run. While the target voltage was 
decreased by steps of 1.0 or 2.0 volts and J; was 
kept constant, J,’ was read for each voltage. All 
runs were taken in the order of decreasing 
voltage, since J, reached equilibrium more 
quickly than in the reverse order. Check runs 
taken with increasing target voltage showed 
breaks in the same positions. At the conclusion 
of each run the constancy of the amplifier was 
checked by repeating readings at certain voltages. 
A slight change in target voltage, too small to be 
observed on the voltmeter, could be detected by 
its effect on J,’. Since J; was held constant to 1 
part in 5000, the limiting factor in the sensitivity 
of the apparatus was the accuracy of measuring 
the voltage, less than 0.1 volt. The oxide coated 
filament used at low temperature and the small 
potential drop along the emitting part of the 
filament make the speeds of the electrons more 


" Nakaya, Proc. Roy. Soc. Al24, 616 (1929). 
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uniform than in most soft x-ray experiments to 
date. 


RESULTs | 


Fig. 2 shows a typical soft x-ray critical 
potential curve. Since J; was held constant 


TARGET VOLTAGE-= 
780 200 


22 


Fic. 2. Typical critical potential curve. Run No. 240, 
(100) target. The uncorrected voltage of each break is 
shown above the curve, followed by a superscript denoting 
the intensity or certainty of the break. 


during each run, the amplifier reading J,’ is 
proportional to J,/J; minus a constant and is 
plotted directly against the voltage. The in- 
tensity or certainty of the break is estimated 
from the change in slope of the curves. Intensity 
(1) signifies the least change of slope which can be 
detected or a break dependent on the positions of 
only a few points and (3) indicates the most 
pronounced breaks. 

The uncorrected voltages of all breaks found in 
20 runs with the (111) target, in the range 90-220 
volts, are tabulated in Table I. Table II shows the 
critical potentials found in 18 runs with the (100) 
target. Readings were started 1.5 to 8 hours after 
heating the targets. The pressure while taking 
readings varied from 1.1 to 2.2X10-? mm in 
different runs, although it was nearly constant 
during each run. No significant changes in either 
intensity or position of the breaks can be 
attributed to these variations of pressure or to 
the amount of outgassing or time since outgassing 
the target. Runs taken during the earlier stages 
of outgassing were so irregular, from other causes, 
that no conclusions can be drawn from them. In 
an attempt to determine whether changed surface 
conditions on the target cause changes in the soft 
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SOFT X-RAYS IN COPPER 


TABLE I. Critical potentials for (111) target (uncorrected). 


Groups including nearly all of the breaks are indicated at the heads of the columns. Breaks falling less than two volts 
outside of these groups are included in the column, in iialics. They are included in the averages and the same weight 
is given to each value regardless of intensity. Breaks not obviously belonging to any group are written between the 
columns. A break missing in any run is indicated by an x where it should appear and is considered as zero intensity when 
computing the mean value. Superscripts denote intensity of the breaks. Numbers 243-257 are the final series of readings 
taken under most carefully controlled conditions. The remaining 6 runs which were taken before certain improvements in 
apparatus were made and were noted as exceptionally steady at time of reading, are seen to agree with the final series. 
Target was heated a total of 157.5 hours for run 243 and 243 hours for the last run 257. For the earlier runs given the 
total heating was from 20 to 74 hours. 
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x 
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1378 
116 


116? 137 
116? 136 


108% 
108 
108° 


1077-117? 


{or } 


voltage 93.3 107.8 116.8 


159.7 177.4 193.1 203.7 


Mean 
int. 2.2 : 2.6 1.5 1.1 


2.3 1.1 


x-ray characteristics, before runs 256 and 257 the 
(111) target was heated to a higher temperature 
than before, until a considerable amount of 
copper evaporated. No definite changes were 
observed in these runs. 

The construction of the x-ray tube did not 
permit an experimental determination of the 
contact potential difference between target and 
filament. Farnsworth has found that for a copper 
target and an oxide coated filament the cor- 
rection is negative and of the order 3 volts. 
Adding 1 volt for drop in the null resistances, we 
find 4 volts, to be subtracted from the voltmeter 
readings, as the most probable correction for the 
soft x-ray critical potentials. 


DISCUSSION 


The average critical potentials for the (111) 
and (100) faces of copper single crystals from 
Tables I and II are shown in the first two 
columns of Table III. Here the most probable 
correction of —4.0 volts is applied. It is evident 
that the same critical potentials appear for both 
targets. In the region 90-220 volts, 14 breaks are 
found and, in all cases, the voltage of the break 
agrees within 1.5 volts for the two targets. 

Critical potentials for polycrystalline copper 
targets published by Compton and Thomas'® and 
by.Richardson and Rao" are shown, respectively, 
in columns 3 and 4. In both cases the corrected 


® Richardson and Rao, Proc. Roy. Soc. A128, 16 (1930), 
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Notation is the same as for Table I. Numbers 233-241 are the final series with total heating of target between 220 and 
250 hours, Others were obtained in earlier runs under best conditions after 100 hours of heating. 


TABLE II. Critical potentials for (100) target (uncorrected). 
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Run 92- 98- 


108- 116- 122- 130- 137- 144- 159- 176- 186- 193- 204- 212- 
No. 93 99 109 118 123 131 138 146 160 178 188 196 206 214 


Final series of readings 


234 

235 

236 130 138% 
237 98? 108° 116 122? 130? ©6138 


98? 1098 123? 1318 


92? 122? 


240 130 


241 932 { oo} 109) 15 


1292 
121 138 


Earlier runs 


203 124 130-137 
2 2 

206 92 { 108? 118% 123 «131136 
207 

208 138? 
221 1383 
222 902 108? 116 123% 1308138? 
223 932 99% 108% 116% 123° 130 

227 922 98% 108% 116 122% 132 138% 


228 122? 130° 


160° 176 1872, 1944 214" 

178% 188 x 206 212? 

178° 187? 212? 
1452 159: 176? 
144 «176? 


143 161 1777 187? 1948 206%) 


1445, «160? 187 x 204 212? 


144° 159 178? 186% 193? 204 213 
1445 160? 177% 187% 193% 204% 214 
1462, 160? 1195? 206 — 214" 
1467-159? 


1457159 176 


Mean 


voltage 92.5 98.2 108.2 116.6 122.5 1304 137. 


7 144.6 159.7 176.9 186.9 194.0 205.2) 213.1 


Mean 
int. 


1.8 2.3 2.6 1.4 1.9 1.7 


1.8 


voltages are given. Both groups of experimenters 
found more breaks for polycrystalline copper 
than were found by the author with single 
crystal targets. This is in agreement with the 
results of others’: ‘ using single crystals. Agree- 
ment between the author’s results and published 
results for polycrystalline copper is only fair 
since, on considering the small separation of the 
breaks found for the ordinary metal in this 
voltage range, it is evident that a considerable 
number of points might accidentally agree within 
two volts. However, it may be said that agree- 
ment of the critical potentials found for the 
single crystals with certain of the critical po- 
tentials for the polycrystalline metal is as good as 
is usually found in soft x-ray measurements of 
different experimenters. 

All of the 14 breaks observed in this experi- 
ment may be separated into 3 groups and 
arranged in a table showing horizontal and 
vertical differences of the order of magnitude 


predicted by an approximate calculation based on 
Richardson’s theory of structure electrons.” 
Some doubt is cast upon this apparent corre- 
lation by comparison of results of other observers 
for metals of similar crystal structure. For 
example, nickel has the same crystal structure as 
copper and the atomic numbers and _ lattice 
constants of these metals differ only slightly, 
hence correction terms neglected in the approxi- 
mate calculation should be of the same order of 
magnitude. If we use the true lattice constants in 
the development of Richardson's formula and 
substitute the atomic number for z (neglecting 
screening effects) the calculated level spacings 
for copper and nickel are very nearly the same; 
but Richardson's" tabulation of Rao’s® critical 
potentials of secondary electron emission from 4 
nickel single crystal shows level spacing con- 
siderably larger than that found for the soft x-ray 


8 Richardson, Proc. Roy. Soc, A128, 63 (1930). 
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TABLE III. Comparison of critical potentials for copper. 


1.5 
13. 
1.9 
1.5 


182.9 
190.0 
201.2 
209.1 


2.1 
2.2 
1.7 
1.9 


181.6 
189.1 
199.7 


209.8 


1 2 3 4 5 
(111) Face (100) Face . 
Corrected Corrected P.C.Cu P.C.Cu Absorption 
Author Author Cc. R. & Kronig" 
Volts Int. Volts Int. Volts Int Volts Int Volts r 
89.3 2.2 88.5 1.8 86.7 1.2 93 6 
92.2 1.6 
94.4 1.8 94.2 2.3 94.5 0.9 101 24) 
98.8 1.6 104 15/ 
103.8 2.6 104.2 2.6 104.5 2.0 — -= 
108.2 1.8 107.3 1.0 
112.8 1.5 112.6 1.4 111.0 1.4 111.0 0.7 116 12 
115.2 1.8 
119.0 1.1 118.5 1.9 117.8 1.5 117.8 0.6 {iar 3} 
126.5 1.5 126.4 1.7 127.1 1.4 128.0 1.5 — — 
132.5 1.3 133.7 1.8 132.9 0.9 133.5 1.8 139 4 
140.9 2.1 140.6 2.2 140.6 1.5 139.5 0.8 { 148 a 
148.4 1.4 148.3 0.9 150 12 
155.7 2.3 155.7 1.9 157.1 1.4 157.4 0.9 162 24 
164.6 1.6 
168.3 1.1 166.3 1.1 170* 36 
173.4 1.1 172.9 1.8 


critical potentials of copper in the present 
experiment. 

It thus seems necessary to seek other possible 
explanations of the soft x-ray discontinuities. 
One such possibility, suggested to the author by 
Dr. Farnsworth, is presented by Kronig’s" recent 
interpretation of fine structure of x-ray ab- 
sorption spectra by means of the quantum 
mechanics of electrons in crystal lattices. This 
theory is based on the idea introduced by Bloch 
that influence of the crystal lattice can be 
replaced by a periodic potential field having the 
same period as the lattice. The theory shows 
that not every value of velocity is possible for 
electrons moving through a crystal lattice but 
rather that an energy spectrum exists, consisting 
of alternating zones of permitted and forbidden 
energies. Kronig has calculated, from the wave 
equation, the perturbed energy of an electron 
travelling in a crystal lattice. It is shown that 


“ Kronig, Zeits. f. Physik 75, 191 (1932). 


there are discontinuities in the energy of the 
perturbed electron when its momentum param- 
eters are related in a certain way to parameters 
a, B, y, which may be considered as defining 
planes in the crystal lattice. Although these 
discontinuities or forbidden energy zones depend 
on direction of propagation of the electron, 
Kronig has shown that there are still finite 
variations in the absorption coefficient corre- 
sponding to each discontinuity, after integration 
over all directions of propagation. In terms of the 
parameters, the energies of these discontinuities 
are where Ais 
Planck's constant, m the electron mass, and a the 
distance between potential barriers equal to the 
lattice constant of the absorbing substance. 
Using values of a, 8, y determined by the 
sequence of planes in the metal lattice, Kronig 
calculated the energies, in electron volts, of these 
discontinuities for copper. Results of these 
calculations are shown in the fifth column of 
Table III. Here r represents the abundance of the 


615 
178.6 1.5 178.4 0.7 
188.8 1.7 188.0 0.6 { 194 
197 24 
ve — — 197.1 1.3 208 18 
{ 204.8 1.4 205.3 0.6 
217 28 
216.8 1.4 217.1 0.5 
) TO 
TS 
or 
as 
ce 
y, 
of 
in 
nd 
ng 
gs 
1e; 
-al 
la 
n- 
ay 


616 MORGAN L. 


discontinuity planes considering permutations 
and changes of sign of the parameters. 

In 4 cases, pairs of Kronig’s calculated values 
are separated by 3 volts or less. If we group these 
pairs as indicated by the brackets in Table III, it 
can be seen that 10 of the 14 critical potentials 
found by the author fall 4 to 7 volts below a 
calculated value or group. This uses all but one 
of the calculated values—the one at 170 volts 
marked with a star. This calculated discontinuity 
is due to the (553) and (731) planes of the 
crystal, planes in which atomic density is very 
small. The four observed critical potentials which 
do not correspond to a calculated value are 104, 
126.4, 173.2, and 182.2 volts. All of these were 
found by other observers with polycrystalline 
copper targets, so it is improbable that they 
might be due to impurities. The break at 104 
volts is very strong, in fact, for both targets the 
intensity is the highest of all the breaks observed. 
Kurth,'® gives 116A, or 106 volts, as the con- 
vergence frequency of the N series for copper, 
which checks closely with this observed strong 
break. 

Similar correlation can be shown between 
Kronig’s table of energy values for iron, which 
has different crystal structure, and soft x-ray 
breaks observed by Thomas’ for the polycrystal- 
line metal. Strikingly good agreement is found on 
comparing Thomas’ critical potentials for nickel 
with the energy discontinuities, which may be 
calculated from Kronig’s table for copper by 
correcting for the difference in lattice constants. 
Thomas found one break not predicted by the 
theory and the calculations show two discon- 
tinuities close together which were not reported 
by Thomas. All except 2 of the 12 breaks reported 
by Thomas fall within 4 volts below a calculated 
value. It was noted above that critical potentials 
for copper found in the present experiment also 
fall below the calculated values by 4 to 7 volts. 
This might be explained by an inner potential 
of the metal which causes acceleration of elec- 
trons on entering the target. Farnsworth’ found 
that the experimental voltage for electron diffrac- 
tion beams from a copper crystal falls below the 
theoretical voltage, so that an inner potential @ 
of the order 25 volts is required. Since it appears 


'® Kurth, Phys. Rev, 18, 461 (1921). 
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that the full value of ¢ is not effective unless the 
velocity of incident electrons exceeds a certain 
value, Farnsworth suggests that ¢ depends on 
depth of penetration of the surface ‘ayer. The 
fact that a smaller inner potential is required to 
correlate soft x-ray data is in accord with this 
suggestion, since the high absorption coefficient 
for soft x-rays requires that they be generated 
near the surface of the target. 

The observed agreement between soft x-ray 
critical potentials and Kronig’s calculated ab- 
sorption discontinuities may be explained as 
follows: The absorption discontinuities are calcu- 
lated for forbidden zones of energy of electrons 
moving through the crystal lattice. Electrons 
having energies falling within these forbidden 
zones cannot enter the target so they are 
reflected, either diffusely or as diffraction beams, 
Consequently, electrons having these energies 
are not available for production of x-rays and a 
minimum in the soft x-ray curves results. These 
minima in the continuous radiation, super- 
imposed on the general upward trend as voltage 
is increased, appear as discontinuities in the 
critical potential curves and account for the 
breaks not predicted by the Bohr theory. Thus 
many of the observed soft x-ray critical potentials 
are not due to true characteristic radiation of the 
atoms, but are characteristic rather of the 
lattice structure of the target. 

The observed agreement between the breaks 
found for the two different crystal faces remains 
to be considered. Since Kronig’s values are 
obtained by integrating over all directions of 
motion of electrons in the lattice, it appears at 
first sight that some of the discontinuities might 
not appear when bombarding electrons fall 
within 15 to 20° of normal incidence on a single 
crystal target, as in the present experiment. 
Obviously one would expect different breaks to be 
missing for different orientations of the targets. 
Failure to observe this difference may be ac- 
counted for by random changes of direction of 
the incident electrons on first striking the 
targets. In Kikuchi's'® experiments on diffraction 
of high-speed electrons from mica, certain lines in 
the photographs are explained as due to electrons 
scattered on first entering the crystal and 


“ Kikuchi, Japanese J. Phys. 5, 83 (1928). 
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subsequently reflected from crystal planes at the 
Bragg angle. This requires an amount of scat- 
tering, on first striking the target, sufficient to 
account for the agreement of the soft x-ray 
breaks for the two targets. It is known from 
diffraction experiments with low-speed electrons 
that only a small fraction of incident electrons is 
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returned as diffraction beams and consequently a 
considerable number must also be scattered 
diffusely. 

In conclusion, the author wishes to express his 
sincere thanks to Dr. H. E. Farnsworth, for his 
many helpful suggestions and criticisms during 
the course of this research. 
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Decay of Luminescence and Light Absorption in Phosphorescent Materials 


D. H. KaBAKJIAN, University of Pennsylvania 
(Received May 1, 1933) 


Effect of continuous irradiation with a, 8 and y-rays 
on the light absorption coefficient of certain phosphor- 
escent materials. With continuous irradiation the light 
absorption coefficient of each sample reaches a character- 
istic saturation value which is a function of the intensity 
of the radiation and the temperature at which the irradi- 
ation takes place, but is independent of the variations in 


the light emission of the material and also of the total 
amount of energy absorbed by it. The experimental curves 
of the rise and decay of phosphorescence and of the changes 
in light absorption coefficient of these materials cannot be 
explained by the generally accepted theories of these 
phenomena. 


HEORIES of luminescence produced in cer- 
tain materials by absorption of various 
forms of radiation can be classified in two groups: 
(1) the “excitation” theory which considers the 
emission of light as due to displacement of elec- 
trons and their subsequent return to normal 
energy levels without any permanent change in 
the material; and (2) the “‘destruction”’ theory 
which considers the emission of light as resulting 
from a chemical change and therefore gradual dis- 
appearance of luminescent centers in the material. 
The second theory was adopted to explain the 
decay of luminescence, in some cases, as a result 
of continuous irradiation. In the case of scintilla- 
tions produced by alpha-particles in phosphor- 
escent zinc sulfide Rutherford’s destruction 
theory seems to be generally accepted.? The 
theory postulates that each active center upon 
being hit by an alpha-particle is destroyed and 
emits a flash of light. Hence the brightness of the 
material, when bombarded by a constant source 
of alpha-particles, decays exponentially with 
time. This law of decay, however, is not verified 
experimentally. 

In studying the decay of self-luminous zinc 
sulfide paints, extending over a period of ten 
years, Walsh came to the conclusion that the 
discrepancy between the destruction theory and 
his experimental curves could be explained if the 


1 Bulletin of Research Council, No. 30, pp. 7-27, 1923. 

?E. Rutherford, Proc. Roy. Soc. A85, 561 (1909-10); 
J. Charitan and C. A. Lea, Proc. Roy. Soc. Al22, 335 
(1929). 


change in the light absorption coefficient of the 
material was taken into consideration. By making 
certain assumptions as to the nature of light ab- 
sorption in an irradiated material he succeeded in 
deriving a theoretical formula which agrees with 
his curves within the limits of experimental error, 
thereby apparently confirming the original de- 
struction theory of Rutherford.* 

According to this modified theory the destruc- 
tion of active centers which results in the emission 
of light, also results in increasing the light ab- 
sorption coefficient of the material. This assumes 
a close relationship between the two phenomena, 
since both result from the destruction of the same 
active centers, the contribution of the non-active 
centers to this change being considered as negli- 
gible. Also, in correcting for the absorption of the 
glass containers, Walsh assumes the change in 
absorption coefficient of glass to follow a simple 
exponential curve of the type e~”', where 
is the transmission factor, pa constant depending 
upon the intensity of the radiation, and ¢ the time. 

In a previous paper attention has been drawn 
to the inadequacy of this modified destruction 
theory to explain luminescence phenomena in a 
number of materials, including zinc sulfide, and 
doubt was expressed as to the validity of the as- 
sumptions on which the theory is based. It will 
be shown further that the luminescence of these 
materials is of a complex nature and a simple ex- 


3]. T. Walsh, Proc. Phys. Soc. London 38, 318 (1926-27). 
*D. H. Kabakjian, Phys. Rev. 37, 1120 (1931). 
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DECAY OF LUMINESCENCE 


ponential decay curve cannot be expected except 
as an approximation. 

The coloration of transparent substances re- 
sulting from irradiation has been extensively 
studied but there is no general agreement as to 
the nature of the phenomenon. In some lumin- 
escent substances this coloration is quite intense 
and would certainly affect the brightness of the 
samples. The variation of brightness due to this 
cause alone has never been studied to my knowl- 
edge. Therefore it was considered desirable to 
make an experimental study of this effect and of 
any relationship that might exist between the 
brightness decay and progressive changes in the 
absorption coefficient of a number of substances. 


MATERIALS AND METHODS OF PROCEDURE 


Several types of glass, some of known composi- 
tion, furnished by Corning Glass Company, were 
used for the absorption measurements. Data 
given in this paper refer to an ordinary micro- 
scope slide glass and to one containing 15 percent 
BaO, 25 percent KO, and 60 percent SiO, 
(Corning Glass No. 6). No luminescence meas- 
urements were made on these. Both luminescence 
and absorption measurements were made on 
crystals of fluorite, calcite, and anhydrous cal- 
cium sulfate activated with manganese. The 
latter was prepared by crystallizing calcium sul- 
fate from a solution in concentrated sulfuric acid. 
These were cut into more or less transparent 
plates. Similar measurements were also made on 
thin films of crystalline powders of willemite and 
zinc sulfide. The powder was mixed with a dilute 
solution of Canada balsam and deposited on very 
thin glass or mica. The thickness of the ZnS films 
varied from 0.03 to 0.05 mm. Changes in absorp- 
tion of very thin glass and a small thickness of 
balsam irradiated with a-rays was found to be 
negligible. 

These materials were exposed to radiation of 
known intensity and changes in their brightness 
as well as in their absorption coefficients were 
made at various intervals of time. 

The sources of radiation consisted of films of 
crystalline radium sulfate applied with a varnish 
on gold plates. The varnish was thin enough to 
absorb only a small part of the a-rays. Although 
there was a slight loss of emanation from the sur- 
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face, after equilibrium had been reached with the 
occluded emanation and the short-life products of 
radium these plates could be considered as con- 
stant sources of a, 8 and y-rays. The radium 
contents mentioned in this paper refer to y-ray 
determinations of radium C in the plates. For 
purely a-ray sources, polonium deposits on silver 
plates were used. A correction for the decay of 
polonium was made in the computations, by us- 
ing 136 days as the half-life period of polonium. 

Absorption measurements were made with an 
ordinary Lummer-Brodhun photometer. The 
transmission of each plate before irradiation and 
after measured intervals of irradiation were de- 
termined, and the absorption coefficient computed 
from the known thickness of the plates. In the 
case of powder films, these measurements do not 
give the absorption coefficient of the material of 
which the film is made but of the film itself. This 
was considered sufficient for the present purpose. 

Luminosity measurements were made with a 
special type of photometer somewhat similar to 
the one described by Meacock and Lambert® for 
the measurement of luminous paints, with, 
however, a Lummer-Brodhun cube for compari- 
son of brightness. The photometer was calibrated 
with a Nutting polarization photometer. This 
calibration introduces an additional experimental 
error in the absolute values of the brightness of 
the samples. It does not affect, however, the rela- 
tive values of brightness with which I was con- 
cerned. The experimental error in these measure- 
ments was about 5 percent. 


RESULTS 


Variation of absorption coefficient 


Fig. 1 shows the effect of irradiating with dif- 
ferent amounts of radium, ordinary microscope 
slide glass; Fig. 2, glass containing 15 percent 
BaO, 25 percent K,O, and 60 percent SiO,. It 
may be seen that in both cases the coefficient of 
absorption, in defined by the relation 
I/Ij=e-** (where I/J is the fraction of the light 
transmitted and x the thickness of the medium 
in cm), increases with time and approaches a 
saturation value. It is also quite evident that this 


°H. F. Meacock and G, E, V, Lambert, J. Sci. Inst. 8, 
214 (1913). 
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Fic. 1. Change in the coefficient of absorption for 
samples of ordinary microscope slide glass exposed to 3.3, 
0.52, and 0.15 mg/cm? of radium. 


| | 
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IN DAYS 


Fic. 2. Change in coefficient of absorption for samples 
of glass containing 60 percent SiO», 25 percent K,O, and 
— BaO when exposed to 4.6 and 0.52 mg/cm* of 
radium. 


7 


COEF FICIENT 


MILLIGR AM DAYS / Cu‘ 


Fic. 3. Same samples as in Fig. 3. Coefficient of ab- 
sorption plotted against the energy absorbed, the energy 
being expressed in milligram-days per square centimeter. 


saturation value is a function of the intensity of 
the radiation rather than of the amount of energy 
absorbed by the material. In Fig. 3 the variation 
of uw is plotted against the energy absorbed in 
milligram-days/cm?* for one of these samples. 

Variation of absorption coefficient of strongly 
luminescent materials was quite similar to that 
observed in glass, with different constants char- 
acteristic of each material. There was always an 
initial rapid rise with an approach to a saturation 
value. An absorption curve for zinc sulfide is 
shown in Fig. 7. 


The effect of temperature 


The maximum value of yu in irradiated glass is 
fairly constant at ordinary room temperature, 
If, however, the temperature is raised the value of 
u decreases reaching different constant values for 
different temperatures, until at a critical tem- 
perature the effect of irradiation completely dis- 
appears. Irradiating glass at this temperature 
does not produce any coloration. Fig. 4A shows 


- 


F ABSORPTION 


IN DAYS 


COEFFICIENT 


Fic. 4. Effect of temperature on the coloration for two 
samples of ordinary microscope slide glass. Both were 
irradiated at 20°C. Before irradiation « was 1.5. Tem- 
perature of sample A was raised to 100°C and kept con- 
stant; that of sample B was first 50°C and later 100°C. 


the recovery at 100°C and B the recovery first at 
50°C and later at 100°C of the transmission of 
glass which had been irradiated at 20°C. It is 
seen that both glasses have completely recovered 
their original transmission. No attempt was made 
to determine the critical temperature at which 
complete recovery takes place. Qualitative ob- 
servations indicate, however, that the critical 
temperature for recovery is different for different 
materials and for glasses of different composition. 
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Brightness decay and light absorption 


Simultaneous determinations of change in 
brightness and in absorption coefficient on the 
same samples show that these two phenomena 
are not directly related to each other. In Fig. 5 is 


+ | 
a © 
: Fluorite exposed to Radium 4 
A- Luminescence ‘ti 
> B - Transmission 

* al 


IN DAYS 


Fic. 5. Luminescence and transmission of a sample of flu- 
orite 0.09 cm thick exposed to 11.7 mg/cm? of radium. 


shown the change in luminescence and in trans- 
mission for a clear plate of fluorite exposed to 
radium. There is a very large initial rise in bright- 
ness with no corresponding change in the trans- 
mission. Fig. 6 shows the rise of brightness in 
calcite followed by a gradual decay. There was no 
measurable change in transmission in this case. 


| 

/_| 
| 
2 
Calcite exposed to Radwm 2 


Tiwi IN DAYS 


Fic. 6. Luminescence and transmission for a sample of cal- 
cite exposed to 11.7 mg/cm? of radium. 


Because of the photosensitivity of zinc sulfide, 
brightness and absorption measurements could 
not be made on the same plate, as in the case of 
fluorite and calcite. They were made on two 
similar plates prepared from the same material 
and irradiated with approximately equal quanti- 
ties of polonium. In Fig. 7, curve A represents the 


~ > 
2 3 
Pi z 
g 

< 

2 

IN DAYS 


Fic. 7. Luminescence and absorption for a crystalline 
deposit of ZnS 0.005 cm thick, exposed to polonium. 


rise and decay of brightness in zine sulfide and 
curve B the increase in its absorption. It is seen 
that the absorption has reached a saturation 
value in about 8 days while the decay of lumin- 
escence continues over 100 days. Therefore the 
relationship between the decay of brightness and 
light absorption assumed by Walsh for this ma- 
terial is not borne out by these observations. 

It is possible to check the modified destruction 
theory in another way. From the curve B, Fig. 7, 
and the known thickness of the film of zinc sulfide 
the value of the absorption coefficient for each 
point on curve A can be determined. It is now 
possible to plot a luminous flux-time curve, elim- 
inating the effect of absorption. Since the bright- 
ness I is given by the relation 


where F is the flux in lumens per cc, x the thick- 
ness, and u the coefficient of absorption, then 
The variation of flux in time is represented in 


Fig. 8. This is not a simple exponential curve as 
required by the destruction theory. 


ZnS exposed to Polonium 
Luminous Flux /cc 


\ 


LUMENS / CC 


IN DAYS 


Fic. 8. Same material as Fig. 7. Luminous flux/cc. 
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Complexity of light emission 

It has been observed by several investigators 
that the light emitted by phosphorescent ma- 
terials consists of fluorescence and often more 
than one type of phosphorescence. This in itself 
precludes the possibility of observing a simple 
exponential decay curve of brightness. A quantita- 
tive determination of the relative importance of 
these factors in zinc sulfide is rather difficult. A 
more suitable substance for this purpose is the 
anhydrous calcium sulfate described above. 
Fairly transparent plates cut from crystals of this 
substance when irradiated with rays from radium 
show fluorescence and phosphorescence to about 
equal degrees. 

A plate of this material was exposed to 4.6 
mg/cm? of radium and its brightness measured at 
successive short intervals. After each measure- 
ment the radium was removed and the residual 
brightness measured immediately. The first set of 
measurements gives the combined fluorescence 
and phosphorescence while the second set gives 
the phosphorescence alone. The difference be- 
tween the two readings gives the fluorescence. 
A small correction was necessary for the decay of 
phosphorescence between the removal of the 
radium and the measurement of the residual 
brightness, an interval of about 30 seconds. 

In Fig. 9, curve A gives the luminescence of 
such a sample, B its phosphorescence, and C, 
which is obtained by subtracting the ordinates of 
B from those of A, its fluorescence. The lumines- 
cence curve, then, is a composite of a decaying 


= 

2 

3| = 

2 


TIME. IN HOURS = 


Fic. 9. Plate cut from a crystal of anhydrous calcium 
sulfate exposed to 4.6 mg/cm? of radium: A, luminescence; 
B, fluorescence; C, phosphorescence, 


fluorescence and a rising phosphorescence curve, 
both modified by an increase in the absorption 
coefficient of the material. The luminescence 
curves shown in Figs. 5-8 are all of the same type 
with variations which must be attributed to the 
relative intensities of their fluorescence and 
phosphorescence factors. 


DISCUSSION 
Light emission 

It has been shown that a simple destruction 
theory of light emission is inadequate to explain 
the brightness decay of the materials discussed in 
this paper even when the change in the absorp- 
tion coefficient of the material is taken into con- 
sideration. Light emission in phosphorescent 
materials by a, 8 and y-rays must therefore be 
considered as due to excitation. The nature of this 
excitation is not clear. The curves 5-8 show, 
however, that the process of light emission does 
not result in permanent changes such as colora- 
tion. 

It has been shown in the previous paper‘ that 
if each active center cannot be excited more than 
once, the number of such centers necessary to 
account for the luminous energy obtainable from 
a given quantity of zinc sulfide is greater than the 
number of molecules in the material. The con- 
clusion to be reached is that the emission of light 
does not destroy the luminescent centers. This 
conclusion is in agreement with the findings of 
Wolf and Riehl® that each luminescent center may 
be excited 1000 to 2000 times before being 
destroyed. 

The curves 5-8 also show that coloration alone 
does not account for the decay of brightness as 
was previously suggested by the present writer. 
A chemical decomposition, not necessarily con- 
nected with the process of light emission might be 
assumed to account for the destruction of the 
luminescent centers. It has been shown that the 
coloration reaches a saturation value. Chemical 
decomposition produced by the rays may be con- 
sidered a process of the same type and therefore 
reach sooner or later an equilibrium value. Since 
these two factors determine the brightness of a 
continuously irradiated material, its final bright- 


® Wolf and Riehl, Ann. d. Physik 11, 103-112 (1931). 
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ness would not necessarily be zero but would 
reach a constant value determined by the extent 
of its coloration and its degree of decomposition. 
This is in agreement with the experimental ob- 
servations on a number of materials where com- 
plete extinction of brightness was never observed 
even when the time was made very large.’ 


Light absorption 

Examination of absorption curves in Figs. 1-3 
shows that no simple explanation such as de- 
struction of molecules of glass or of a limited 
number of centers due to impurities dissolved in 
it, is adequate to explain the phenomenon. The 
first is eliminated by the fact that the curves are 
not simple exponentials, and the second by the 
failure of the samples to attain the same satura- 
tion value of absorption with different intensities 
of radiation. 

If the coefficient of absorption is considered as 
proportional to the number N of absorbing cen- 
ters, and the formation of these centers is similar 
to ionization and recombination of ions in air, 
then N would be given by the equation. 


dN/dt=kR—aN?, 


where R is the intensity of radiation and a is the 
recombination constant of the ions. The solution 
of this equation is 


N=C(e2@at—1)/ (e241). 


C, the saturation value of the number of absorb- 
ing centers, is a function of the intensity of radia- 
tion. A theoretical curve from this equation does 
not agree very closely with the experimental 
curve. If a limited number of potential absorption 
centers be assumed to exist in the material, then 


dN/dt=kR(m—N)—aN?, 


where N can never be larger than m. The solution 
of this is of the form 


where 8B = [4akRm+(kR)* }'. A theoretical curve 
obtained from this equation agrees well with the 
experimental curve of absorption. The constants 


7J. A. Rodman, Phys. Rev. 23, 4 (1925). Leon E. Smith, 
Phys. Rev. 28, 3 (1926). G. S. Gessner, Phys. Rev. 36, 2 
(1930). 
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m and C, however, as determined from the curve 
give a value for C which is larger than m and 
therefore inconsistent with the original assump- 
tion. It is also very difficult to explain on the 
ionization theory the permanency of coloration in 
glass at ordinary room temperatures. 


Decay of phosphorescence 


Similar difficulties are encountered in explain- 
ing by some such simple theory the decay of 
phosphorescence in luminescent materials after 
the removal of the exciting source. In Fig. 10 the 
continuous curve represents the decay of phos- 
phorescence in a calcium sulfate plate which had 
been previously irradiated with 4.6 mg/cm* of 
radium until maximum brightness was attained. 


ad CaSO, exposed to Radwm 
Decay of Phosphorescence 
A~ Experimental Curve 
Qe Theoretical “ 
a 
B- 
Je 
= 
\ 
Se 
— 
Time HOURS 


Fic. 10. Decay of phosphorescence for the same sample as 
in Fig. 9, 


If the luminescence is due to excited electrons 
returning to their normal energy levels, as is 
usually assumed, the brightness should decay in 
accordance with the equation® J = Ae~*', where J 
is the brightness, ¢ the time, and A and ¢ are 
constants for a given sample. If, however, 
luminescence results from the recombination of 
ions, the brightness should decay in accordance 
with the equation® J=B/(at+c)*, where B, a, 
and ¢ are constants for the given sample. Curves 
corresponding to these two equations of decay 
are given in dotted lines in Fig. 10. Neither of 
these agrees even approximately with the ex- 
perimental curve. 

These experiments seem to indicate that in 
phosphorescent materials irradiated by a, 8 
and y-rays, light emission, coloration, and chem- 


*H.S. Allen, Photo-Electricity, p. 222, 1925. 
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ical decomposition take place more or less inde- 
pendently of each other and are not easily ex- 
plainable as simple monomolecular transforma- 
tions. A number of simple exponentials with 
different constants might be found to fit each ex- 
perimental curve. Further investigation, how- 
ever, seems necessary to determine these and 
their physical significance. 

In conclusion, the author wishes to express his 
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the photometric measurements and to Dr, 
Horace C. Richards for assistance and criticism of 
the work, also to Corning Glass Company who 
furnished the glass specimens. This investigation 
was made possible by a special grant of the 
Faculty Research Committee of the University 
of Pennsylvania which is hereby gratefully 
acknowledged. 
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The Theory of Forced Double Refraction 


K. F. HeERZFELD AND R. H. Lee, Johns Hopkins University, Baltimore 
(Received May 10, 1933) 


An attempt is made to better understand the nature of 
the double refraction in crystals by calculating the forced 
double refraction induced in rocksalt and KCI crystals 
through elastic stress (photoelasticity). In a previous 
paper only the change in the Lorentz-Lorenz force had 
been taken into account. In this paper is also considered 


the “‘potential dip’’ which is due to the forced anisotropic 
arrangement of the ionic charges. The formulas are 
developed according to classical theory and to wave 
mechanics but they contain too many unknown constants 
to be checked against experiment. 


INTRODUCTION 


HE question how the anisotropy of non- 

cubic crystals results in double refraction 
is not yet solved. Ewald' who investigated the 
question first assumed that the anisotropy of 
the Lorentz-Lorenz force was alone responsible. 
This anisotropy results from the fact that in 
different directions from a particle the sur- 
roundings are different. He could give only rough 
estimates as insufficient data were known at the 
time. Bragg? and Hylleraas* calculated double 
refraction for carbonates, nitrates, alum, quartz, 
Hg2Cl2, Hg2Bre and succeeded in getting approx- 
imately correct values. They took into account 
not only the anisotropy of the Lorentz-Lorenz 
force but also that of the charges of the surround- 
ing atoms. They could adapt the values which 
they ascribed to the polarizability of the con- 
stituent ions to fit their results. 

If an isotropic, transparent medium is sub- 
jected to pressure perpendicular to two parallel 
surfaces it becomes double refractive. If one 
chooses as such a medium sodium or potassium 
chloride one has a cubic crystal made artificially 
double refracting by changing the structure by 
a known amount and has to deal only with ions, 
the refractivity of which is known from the 
investigation of the undisturbed state. For this 


'P. P. Ewald, Thesis, Munich, 1912; Ann. d. Physik 
49, 1 (1916). 

*W. L. Bragg, Proc. Roy. Soc. A105, 370 (1924); A106, 
346 (1924). 
*E. Hylleraas, Zeits. f. Physik 36, 859 (1926); 44, 871 
1927), 


reason Herzfeld‘ calculated the double refraction 
resulting from the anisotropy of the Lorentz- 
Lorenz force (as had been done before in an 
unpublished thesis by Heine®) and compared it 
with experimental results of Pockels* and Maris.’ 
No agreement was, however, reached. In a dis- 
cussion with Professor Ewald and Dr. Hermann 
of Stuttgart it was pointed out that the success 
of Bragg and Hylleraas depended on the fact 
that these two investigators took into account 
what Ewald called the potential dip. In the case 
of an ionic lattice the field from the electrical 
(excess) charges of the surrounding ions forms a 
potential dip at the place of the ion which we 
have just considered, a dip which changes the 
restoring force on an electron subjected to the 
influence of an outside electric wave. Hence the 
anisotropy of arrangement due to the pressure 
might result in an anisotropy of the potential dip 
and, therefore, in an anisotropy of the restoring 
force and of the refractive index itself. These 
effects are calculated here for a sodium chloride 
type lattice under the action of pressure. It 
turned out that in the undisturbed lattice the 
potential dip is zero to the order considered here. 


GENERAL THEORY OF ForRCED DOUBLE 
REFRACTION 


Following the notation of Herzfeld we consider 
a lattice of the rocksalt type compressed in the 


4K. F. Herzfeld, J.O.S.A. 17, 1 (1928). 
®W. Heine, Thesis, Munich, 1922. 

®*F. Pockels, Ann. d. Physik 39, 440 (1890). 
7H. B. Maris, J.0.S.A. 15, 194 (1927). 
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z direction. The outside electrical field might 
have either the z or the y direction. We call the 
electrical moment induced in a positive ion p)\+ 
and p,, respectively in the two cases. For the 
negative ions we have similarly p,;- and p,_. If 
we have N ions of each kind per cubic centimeter 
the total polarization is 


(1) 
The connection with the refractive index is 
(2) 


where D is the electric displacement. Further- 
more the moment ? of an ion is given by 


b= (3r/4r)(D+E’), (3) 


where EF’ is the field at the place of the ion due to 
all the other ions and is, in the undisturbed 
lattice, given by 


E' = —4nP+(4/3)aP, (4) 


where the second term on the right side of the 
equation is the Lorentz-Lorenz force. 3r/47 is 
the polarizibility of the ion and is given by 


r=) (5) 


When pressure is applied, the lattice is com- 
pressed in the z direction so that the lattice 
distance becomes do(1+ 6,,) in the z direction and 
ao(1+6,) in the x and y directions. That modifies 
the factor (4/3) in (4) as well as 7 in (3). 

Assume that in the parallel case the Lorentz- 
Lorenz force on a positive ion gives 


(6) 


where 
6= 64. 


The coefficients of and are different - 


because the positive ion considered is not 
arranged geometrically in the same way in re- 
spect to other positive ions as in respect to nega- 
tive ions; but for /’),_, the force on negative ions, 
ps and p_ are exchanged. 


CALCULATION OF LORENTZ-LORENZ FORCE 


FE’ is due to the action of all other ions, each 
polarized in the z direction by a homogeneous 
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field, on the ion considered, which we locate at 
the origin. The potential of a dipole is 
— p(d/dz)(1/r). Accordingly the resultant poten- 
tial at the origin is 


Vi= —p,(0/d2) ~ (1/r) — p-(8/d2) (1/r), (7) 


where the summations have to be extended, 
respectively, over all positive and negative ions, 
excluding the one at the origin. It follows that 


V,/dz) 


From this, Herzfeld has calculated 
—0.585, as,= —0.081, 
=1.248, — 2.255. (9) 


CALCULATION OF THE POTENTIAL Dip 


The potential due to the excess charges is 


(1/r)—e (1/r). (10) 


The electric field E’ is —edVe2/dz and is zero at 
the origin, because of the symmetry of the 
arrangement. The restoring force on an electron 
k in the z direction is defined by 


(—e)E.= (11) 


where z is the small deviation of the electron 
from the origin. Therefore 


PVs a? 
(12) 


oz? +r 02? 


Comparing with (8), one sees that no new coef- 
ficients are needed but that 


k,/e= —42rNela—B ]é. 


There, k and a—§ take simultaneously either the 
index | or the index |. That applies to positive 
ions. For negative ions the ions of like sign (i.e., 
contributing the coefficient a) have the opposite 
charge, k_ has the opposite sign. 

For an undistorted rocksalt lattice the charges 
contribute nothing to the restoring force. 

The task of calculating the influence of the 
additional restoring force on r is different 
according to the classical and the quantum 
theory. 
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CLASSICAL THEORY 
The restoring force on one electron without the potential dip, is 


k® (13) 
The restoring force in the distorted lattice is then 
k°+k Ap,’). (13’) 
The corresponding addition to each term of (5) is 
3am(v;?—v?) vj;?—v? v?) 402m 3am(v;—v?) rm(v;?— v*) 


where the sign is the sign of the charge of the ion. 


According to Herzfeld and Wolf, it is possible to represent the mol refraction of rocksalt and sylvin 
by three terms 


[(m?—1)/(n?+2) ]=(e®N/3am) ([fi/ (v1? — f2/(v2? — )). (S’) 


Of these, the third was attributed to the positive ion, the first two to the negative. If we change the 
numbering of the terms and call them (including the common factor outside the bracket) P, Ps, Ps 
designate the second term in (5’) by P; (which is by far the largest and is due to the negative ion), 
the first term P»2 and the last P, we can, upon introducing the values a and 8 from (9), write for 
n\\'—n,’, the difference due to the potential dip, 


4.007 pelt Py 


=— (15) 
(n?+ 2)? fi P, Vo? — ve P, v,?— py? 
since 
(n?—1)/(n?+2) 
Division of Eq. (15) by (m?—1)?/(n?+2)? and multiplication by (m?—1)?/2n gives us finally 
(: 
(n?—1)? 4.007 P, Py 
= — (16) 
2n fi P; 
| 
P, P, 4 


It should be noted that in deriving Eqs. (6), (15) and (16) no account was taken of the change 
in N, the number of molecules per cubic centimeter. To do so would have introduced a term in 8. 
For both kinds of crystal |6| is of the order of magnitude of 10-*. Hence & is negligible compared 
with 6. 

CALCULATIONS AND TABLES 


Numerical calculation of Eq. (16) has been made for three wave-lengths by using the formulae 
of Herzfeld and Wolf. The value of 6 was obtained from the paper by Herzfeld as were the terms 
calculated by him from Eq. (6) which we have added to our result to obtain the complete difference 
(ni; —n,). The experimental values for this difference were also taken from the same paper. 

Two tables are included, the first for sodium chloride, the second for potassium chloride. Column I 
gives the wave-length for which the calculations were made; column II gives the numerical values 
of the terms calculated by Herzfeld; column III gives those calculated here from Eq. (16) mul- 
tiplied by 107; column IV gives the sum of II and III, namely (,;—m,)10"; column V gives the 
corresponding experimental results; lastly, column VI gives the ratio of V to IV. 
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TaBLE I. Sodium chloride. Negative pressure of 1 kg per | Taste II. Potassium chloride. Negative pressure of 1 ke 
square centimeter. per square centimeter. 
I II Ill IV V VI II Ill IV V VI 
Last two Last two 
terms of terms of 
Eq. (6) Eq. (16) 107 (mj; Obs. Eq. (6) Eq. (16) (mj; —ma)10" Obs. 
x10" x10" Calc. Obs. Calc. «107 x 10° Calc. Obs. Cale. 
4358 1.382 —12.388 —11.006 2.197 —0.20 | 4358 — 12.98 — 10.43 — 23.40 —2.973 0.13 
5461 —11.419 —10.234 2.059 —0.20 | 5461 ~—12.67 9.70  —22.37 —2.900 0.13 
—11.245 —10.093 1.922 —0.19 | 5790 —12.56 — 957 —22.13 —2.895 0.13 


There is no agreement between experiment and theory except the order of magnitude. Experiments 
made in this physical laboratory recently by Mr. A. Smith and calculations by Joseph E. Mayer* 
seem to indicate that the results which Herzfeld and Wolf got for the position of the main absorption 
by evaluating the dispersion formula are wrong. As long as no better values are known, it is impos- 
sible to say how far this accounts for the failure of the calculations shown in the tables. 


QUANTUM THEORY 
Consider an ion in the undisturbed crystal. The Schrédinger equation of a dispersion electron is 
Ay + (V—E)po=0. 
We now add the perturbation function due to the potential dip, 
V! = (4N/2)e76[ (a; (aa — Ba) (x? +?) ] 
= F — 8); cos? 3+ 
= cos? —2,174r? J. (17) 


In the last line the sign is the sign of the ion. This perturbation function has cylindrical symmetry 
and accordingly the appropriate undisturbed wave functions are the ones that occur in the Stark 
effect. 

The unexcited state has, of course, spherical symmetry and is, therefore, not degenerate. Call the 
undisturbed wave function yo, then the ground state is not split up. 

We assume that the unperturbed wave function can be written 


— 
4n(1+m) ! 


) Py(cos F(j, 1) (18) 


where F(j, 1) depends on r only and is so normalized that 


fi F(j, 1) |*r’dr=1. (18’) 


The perturbation matrix is given by 


V(j, 1, m; k,l’, m’) = (19) 


From the form of V’ (17) it follows that (19) is different from 0 only if m’=m; l’=1 or Il’ =1+2. 
Furthermore, the transition probabilities in (5), f(¢, 1, m; j, I’, m’) are different from zero only if 


8 Joseph E. Mayer, J. Chem. Phys. 1, 270 (1933). 
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'=1+1 and either m’=™m (for light polarized parallel to z) or m’=m+1 for light polarized normal 
to z. Each member of the mol-refraction formula (5) will be changed by 


A = (- 4+), (20) 
where we have put 
~¥;| P05 9) (21) 


we have the same selection rule for the electric moment p as we stated before for f. 

It is true that there ‘will also appear formerly ‘‘forbidden”’ lines, but their amplitude is of the 
first order in V’, their intensity therefore quadratic. 

Now the application of the general theory of perturbations® together with the properties of V’ 
and the selection rules leads to the following formula: 


V(0;k, 0,0) V(k, 1, m;j, 1, m) 
A| p(0; 7, 1, m)|*=2p(0; j, 1, m) 0,037, 1, m) — Ww p(0; k, 1, m) 
k j 
V(0; Rk, 2, 0) 


k,2,0;j,1,m)}. (22 
j m)| (22) 


Here we have made use of the fact that the p's can be chosen real.'® From (22) we get therefore :!! 
f(0; 3) 
— 


i y?—y? (v;?— v?)? 


CVG, 1, m; 7, 1, m) — V(0; 0) ]p(0; 7, 1, m) 


W, — Wo, V(0; k, 0, 0) 
p(0; 7, 1, m) p(k, 0, 0; 7, 1, m) 


$23, 
(0: 1, m)p(k, 2, 0: j, 1, m) 
1, m)p(O; k, 1, m) p(0; 1, m) p(k, 2,0; 7, 1, m 
W.-W; Wo—Wi,2 


vivetv? 


1, m; 7, 1, m)p(0; 7, 1, m)p(O; Rk, 1, m) 


p2 
+ V(0; ———— p*(0; j, 1, m) 


~ V(0; k, 0, 0)p(0;j, 1, m) p(k, 0, 03, 1, m) 
—— V(O;k, 2, 0)p(0;j, 1, m) p(k, 2,05 j,1,m). (23) 
i ok 
Define now 


r2(j, 1, m) = l, m)r? F(R, l', m)r*dr, A= +F2rNe’*s. 


*See, for example, M. Born and P. Jordan, Elementare © J. H. Van Vleck, Electric and Magnetic Susceptibilities, 
Quanten Mechanik, p. 198, Eq. (22), Berlin, 1930. p. 273, Oxford, 1932. 


“The prime on the double sum means j #&. 
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We have then 


V(k,1,0;7,1,0)=A 0,231 r°(k, 137, 1), 
V(k, 1, 1;7,1, +1)=—A 1,373 r°(k, 137, 1), (24) 
V(0; k, 2, 0) =A(2/15)5!-4,007  72(0; R, 2). 


We have for light polarized in the z direction, m =0; for light polarized in the x-direction, m= +1 
and, as each of these corresponds to circularly polarized light 


=} p7(0; j, 1, +1)+p7(0; J, 1, —1)} = i, +1) =p°(0; J, 1, 0) =p,’, (25) 


the latter because of the ‘‘law of spectroscopic stability.” 

Therefore, if we form for (23) the difference between light polarized parallel and at right angles to z 
(designating this difference with A’), the second and third terms drop out. 

Furthermore, we have” 


p*(k, 2, 0; 7, 1,0) =2/5 p.7(k, 2; J, 1), 
2,037, 1)=31 2,037, 1, +1) +p°(k, 2, 057, 1, —1)} =1/10 p.7(R, 2; 7, 1). (26) 


p(k, 2,0;7, 1, +1) and p(0;j, 1, +1) have opposite signs because the increase of m by 1 is in the 
former case accompanied by a decrease, in the latter by an increase in /. The corresponding is true 
for m= —1. On the other hand p(k, 2, 0; 7, 1, 0) has the same sign as p(0; 7, 1, 0), because no change 
in m occurs. 

The final formula for the difference in refractive index due to the potential dip, apart from the 
Lorentz-Lorenz force, is: 


2—1)(n?+2) 0; 7, 1))-! 0; 7,1) f(0;k, 1)\! 
6n 5h i ak Vk 
(7,1; k, ‘) r(0;k, (27) 
(v;? — v?) (v4.2 — v?) v?(0; k, 2)v(7, 1; Rk, 2) 


Here a; is +1 if the absorption lines j and k both belong to the negative ion 
is —1 if the absorption lines 7 and k both belong to the positive ion 
is zero if the absorption lines 7 and k belong to different ions. 


While the first sum includes j=, this is excluded in the second sum. Comparing with the classical 
expression (16), one sees 
A’qu 1 44Ne?_ n?+2 Ne?r?(1; 1) 
fi (1; 1) =0,84 - fi. 
A’kl 3(m?—-1) Sh hv; 


As N=1/2r,?, where ro is the distance between two adjacent ions 


A’qu n?+1e? r°(131) 1 
~0,42 


— fy. 
—1 ro fo" hv; 


Now hy, is somewhat larger than the lattice energy 1.74 e?/ro, therefore this expression is probably 
<1, although of this order of magnitude. 
One sees from Tables I and II, that (27) should be about —1/10 of column III for rocksalt, —1 of 


column III for KCI to give the experimental values. Whether that is possible cannot be determined 


#2 Born and Jordan, reference 9, p. 160. 
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until more is known about the absorption spectrum and the atomic states in crystals, so that the 
different expressions can be calculated. The difference between the classical and quantum results is 
primarily due to the use of a more general dynamical model in the latter. When the matrix amplitudes 
involved in (27) are specialized to the harmonic oscillator, the quantum formula is the same as the 
classical one. 

We wish to thank Professor Dieke and Dr. M. Goeppert-Mayer for helpful discussion and advice, 


and particularly Mr. Serber and Professor J. H. Van Vleck of the University of Wisconsin for im- 
portant corrections. 
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Magnetic Rotation Spectrum of the Red Bands of Sodium 


W. R. FREDRICKSON AND CARL R. STANNARD, Syracuse University 
(Received June 30, 1933) 


The magnetic rotation spectrum of the red bands of 
sodium shows a large number of doublets. A doublet is 
interpreted as made up of an R-branch and a P-branch 
line having a common upper XK’ rotational level. The K’ 
values for the doublets are calculated and it is found that 
the rotational levels in the various initial vibrational states 
from 0-12 which give rise to magnetic rotation doublets are 
53, 45, 33, 14, 60, 52, 43, 31, 16, 58, 51, 40, 28 respectively. 


That just these levels give rise to magnetic rotation lines 
suggests that these levels are perturbed. Since a ‘II state 
is predicted to lie below the upper '= state we have perhaps 
here evidence of a perturbation between a 'Z and a *II state, 
To account for the magnetic rotation doublets it has been 
necessary to revise the vibrational numbering, increasing 
the initial state numbering by three. 


OME twenty-five years ago R. W. Wood! and 
his collaborators reported the magnetic rota- 
tion spectrum of the blue-green and red bands of 
sodium. The experimental procedure was to pass 
light through two nicols between which was 
placed a tube of absorbing vapor, the latter being 
arranged between the poles of an electromagnet 
in such a way that light passed through the tube 
in a direction parallel to the magnetic field. When 
the nicols were set for extinction, the magnetic 
rotation spectrum was produced upon the excita- 
tion of the magnetic field. 

F. W. Loomis? has shown that in the blue-green 
system most of the magnetic rotation lines cor- 
respond to the first lines in the R-branch and 
from this fact he has been able to make a very 
accurate vibrational analysis of this system. The 
blue-green bands are due to a 'II—'2 electronic 
transition and Loomis and Nusbaum’ have 
shown that the same electronic transition in 
other molecules shows the same magnetic rota- 
tion effect. Serber,‘ in a theoretical discussion of 
the problem, concludes that the magnetic rota- 
tion effect is associated with the Zeeman effect 
in such a way that lines showing large Zeeman 
splitting should also show large magnetic rota- 


1R. W. Wood, Phil. Mag. 10, 408 (1905); 12, 499 (1906); 
Wood and Hackett, Astrophys. J. 30, 339 (1909). 

2 F, W. Loomis, Phys. Rev. 31, 323 (1928). 

’Loomis and Nusbaum, Phys. Rev. 38, 1447 (1931); 
39, 89 (1932). 

*Serber, Phys. Rev. 41, 489 (1932). 


tion. The magnetic rotation lines in the 'II—' 
bands agree with this conclusion. 

The red sodium bands, however, are due to a 
1Y—'Z transition and should therefore show no 
Zeeman or magnetic rotation effects. Yet one 
finds that the red bands show more magnetic 
rotation lines than the blue-green and that the 
rotation in the red system must be greater since 
the lines can be obtained with weaker magnetic 
fields. The red rotation lines also differ from the 
blue-green in that they show no relation to the 
band origins. 


EXPERIMENTAL RESULTS AND INTERPRETATION 


The magnetic rotation spectrum was obtained 
by the same method as used by Wood with the 
exception that a long solenoid was substituted 
for the electromagnet. The entire spectral region 
from the D lines to A8000A was photographed 
with a one meter grating in a Rowland mounting 
and the region from the D lines to \7400A was 
photographed in the first order of a twenty-one 
foot grating in a similar mounting. Some two 
hundred lines were measured on the low dis- 
persion plates while over eight hundred were 
measured on the high dispersion plates. The great 
increase in the number of lines on the latter is 
mainly due to fact that what appears as a single 
broad line at low dispersion shows up as a group 
of lines at high dispersion. 

Fig. 1A shows the spectrum at high dispersion. 
The spectrum falls into two distinct sections. 
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Fic. 1A. High dispersion magnetic rotation spectrum of sodium. 


The first, typified by the middle strip, extends 
from the D lines to about 6800A and shows a 
large number of doublets. The second, designated 
by the bottom strip, shows groups of lines form- 
ing doublets. The top strip covers exactly the 
same region as the middle but the temperature 
of the absorbing vapor in this case is consider- 
ably less than for the middle strip. This top strip 
shows only two series of doublets, a and 6, which 
increase in doublet separation as the series 
progresses to the red. Most of the high-wave 
number components of doublet show negative 
rotation (see Table I) while the other component 
shows positive rotation which, according to 
Serber,‘ should correspond to an R and P branch 
line, respectively. 

The immediate interpretation of a doublet was 
that it consisted of an R and P line having a 
common upper K’ level. Two interpretations of a 
series seemed possible. The first being that a 
series arises from transitions (in absorption) to a 
common upper vibrational state from lower 
vibrational states (i.e., doublets similar to those 
observed in fluorescence series). If such were the 
case the doublet spacing should decrease slowly 
as the series progressed to the red while it was 
observed that the spacing rapidly increased. The 
second interpretation was that since these 
doublets, shown on the top strip, appeared when 
the vapor was relatively cold they must all be 


due to transitions from low final vibrational 
states to various upper levels. 

This latter interpretation could be tested in 
the following way: The doublet separation is 
R(K” —1)—P(K""+1) and therefore can be 
equated to the expression for the final state com- 
bination. Thus, 


doublet separation = A, F” = (4K"’+2)B” 


Using the B’s and D's obtained from the rota- 
tional analysis of the blue-green® and red® bands, 
we could calculate the value of K” which is 
numerically the value of K’ of the upper level 


12,0 11,0 


Cold vapor 


B- 


Hot vapor 


12,0 
Fic. 1B. Part of Fig. 1A showing that the 12,1 band 
appears only on the hot vapor plate. 


12,1 


5 Loomis and Wood, Phys. Rev. 32, 223 (1928). 
6 Fredrickson, Phys. Rev. 34, 207 (1929). 
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which is common to both members of a doublet. 
Then taking the P-branch component of a 
doublet together with the value of the B’s and 
D's and increasing the K value by one since 
for the P-branch K” = (K’+1), the origin of the 
band to which the doublet belonged could be 
calculated from the expression’ 


vp = (B’+B")K" 
—2(D’+ 


Since it is known that the R-branch turns at very 
low K values, the head of the band should be 
very close to the origin and thus the calculated 
origins should agree with the heads of the bands 
as measured in absorption. 

In Table I are listed the nine doublets observed 
in the top strip of Fig. 1. The first column gives 
the wave numbers of the lines, the second the 
rotation observed by Wood, the third the doublet 
separation, the fourth the calculated upper K’ 
level, the fifth the calculated origin of the band, 
and the sixth the observed absorption band. The 
agreement between the fifth and sixth columns 
is good,—only one doublet (No. 7) showing any 
disagreement. As seen from the last two columns, 
what appeared on the photographic plate as two 


TABLE I, Magnetic rotation doublets appearing on top 
strip of Fig. 1. The sign of the rotation is taken from 
Wood and Hackett's paper, Astrophys. J. 30, 339 (1909), 


Vibra- 


tional 

Sep- vy assign- 

ara- (absorp- ments 

Doublets Rot. tion K’  (cale.) tion) old new 

1 14973.5 36.2 60 15121.3 151215 1,0 4,0 
14937.3 + 

2 15125.1 — 31.9 52 15236.7 152368 2,0 5,0 
15093.2 + 

15273.3 265 43 15348.8 153493 3,0 6,0 
15246.8 + 

4 154239 19.5 31 15461.1 15461.7 4,0 7,0 
15404.4 + 

5 155661 — 10.3 16 15574.1 155744 5,0 8&0 
15555.8 + 

6 15537.2 — 35.0 58 15684.9 15684.0 6,0 90 
15502.2 + 

7 156854 — 31.3 51 15798.6 15794.8 7,0 10,0 
15654.1 + 

8 158238 — 24.7 40 15902.8 15904.0 8&0 11,0 
15809.1 + 

9 15981.1 17.5 28 160130 16012.3 90 12,0 
15963.6 — 


7 Mulliken, Rev. Mod. Phys. 2, 69 (1930). 
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Fic. 2. Square array of vibrational analysis. Circles 
represent the magnetic rotation data, crosses the absorp- 
tion data. 


separate series is really only one progression 
arising from a common vibrational state v’ =0. 

If this interpretation is correct, bands with 
final vibrational states other than zero should 
likewise show similar doublets. For example, if 
the rotational level K =28 in the upper vibra- 
tional state v’=12 gives rise to a doublet in the 
12,0 band, it should also give rise to almost 
identical doublets in the (12,1) (12,2), etc., bands. 
There is evidence that these doublets do appear 
as indicated in Fig 1B. Here the 12,1 band 
appears quite clearly on the hot vapor plate 
while it can scarcely be detected on the cold 
vapor plate. The bands with 2” other than zero 
which have been identified are shown schematic- 
ally in Fig. 2. Not many doublets are assigned to 
these bands but this failure of the bands to make 
their appearance is really what might be expected 
since the vibrational distribution follows a wide 
Condon parabola and, as seen from Fig. 2, one 
likewise finds only a few absorption bands. 

In extending the interpretation to the region 
to the red of A6800A one has difficulty because 
the lines appear in groups. By plotting the 
groups to a large scale it was observed that the 
lines in the various groups were somewhat 
similarly spaced. Taking a line representing a 
so-called center of gravity of a group, one obtains 
doublets whose doublet separation is practically 
constant and whose spacing between the suc- 
cessive doublets indicates that they belong to 
low 2” values. Following the procedure for the 
doublets mentioned above, the origins were cal- 
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TABLE II. Magnetic rotation doublets and absorption band heads in the region beyond 6800A. 
Sep- vH Sep- vH 
ara- Yo (absorp- VH ara- Vo (absorp- vH 
v’ Doublets tion K’ (calc.) tion) (calc.) |v’ Doublets tion K’  (calc.) tion) (calc.) 
0 1 14499.1 14501.2 
0 2 14241.5 31.2 53 143428 14344.9}1 7 136386 28.3 45 13701.5 13703.2 
14210.3 13610.3 
0 3 14088.6 31.5 53 14187.0 14190.6 14190.1/1 8 13489.2 27.9 45 13550.4 13556.2 13556.0 
14057.1 13461.3 
0 4 13938.8 32.1 53 14034.2 14038.0 14036.8|/2 1 146961 204 33 14735.0 14734.7 14734.2 
13906.7 14675.7 
0 5 13789.7 31.8 53 13882.9 13886.8 13885.1/2 2 14539.8 20.2 33 14577.6 14575.9 14577.9 
13757.9 14519.6 
0 6 13642.2 31.6 53 13732.6 13731.8 13734.8|2 3 143870 20.7 33 14423.3 14421.2  14423.1 
13610.6 14366.3 
0 7 13496.2 31.7 53 13584.3 13581.6 13586.2|2 4 14235.0 20.2 33 142708 ——— 14269.8 
13464.8 14214.8 
0 8 13351.6 31.6 53 134330 ——— 13439.0|2 5 140853 200 33 141206 141181 
13320.0 14065.3 
0 9 13182.9 31.9 53 13296.0 13292.4 13293.6|2 6 13933.6 19.8 33 139680 13967.4 13967.8 
13211.0 13913.8 
1 1 14545.6 28.0 45 14619.6 14618.0 14618.2}2 7 13786.2 19.8 33 13819.6 13818.2 13819.2 
14517.6 13766.4 
1 2 14387.0 27.9 45 14459.3 14461.6 14461.9}2 8 13636.6 19.5 33 13669.2 13665.0 13672.2 
14359.1 13617.1 
1 3 14235.0 28.0 45 143060 14308.7 14307.1)}2 9 13492.5 20.3 33 135233 ———  13526.6 
14207.0 13472.2 
1 4 14085.2 28.1 45 14153.5 14154.5 14153.8)3 0 15003.3 89 14 15008.6 15006.7 15008.0 
14057.1 14994.4 
1 5 13933.6 27.9 45 14000.4 14003.0 14002.1|3 1 148454 9.2 14 14850.1 14850.4 14850.2 
13905.7 14836.2 
1 6 13781.6 27.3 45 13847.0 13848.3 13851.8|3 2 146886 94 14 14692.8 14691.9 14693.9 
13754.3 14679.1 
TABLE III. Additional band heads from absorption plates. 
v (obs.) v (cale.) v’ v (obs.) v (cale.) v’ v (obs.) v (cale.) 
0 11 13008.3 13007.4 3 16 12670.3 12669.9 6 19 12619.4 12619.8 
> te. 12864.7 12866.8 3 17 12534.3 12537.7 7 16 13123.2 13123.5 
1 13 12846.2 12844.8 4 9 13760.0 13757.1 7 17 12993.1 12991.3 
1 14 12706.4 12707.5 4 11 13473.2 13470.9 7 18 12861.6 12860.8 
» 13239.1 13240.4 4 15 12918.5 12918.3 8 17 13101.5 13102.8 
2 12 13098.6 13099.8 4 16 12783.5 12784.4 8 19 12842.3 12843.6 
2 13 12960.2 12960.8 4 18 12522.3 12521.8 9 18 13087.5 13083.1 
3 & 13787.6 13788.2 5 11 13585.1 13584.7 9 19 12960.2 12954.4 
3 9 13641.5 13641.6 5 12 13442.6 13444.1 9 21 12706.4 12702.3 
3 10 13499.4 13498.8 5 16 12899.3 12898.2 10 20 12942.3 12937.5 
$ 12 13214.2 13285.8 5 17 12764.8 12766.0 10 22 12693.7 12689.0 
3. 13 13076.2 13076.8 5 19 12507.2 12506.7 12 20 13156.1 13155.2 
3 14 12942.3 12939.5 6 15 13145.4 13145.2 13 21 13140.8 13137.8 
ee. $s 12802.9 12803.8 6 16 13008.3 13011.3 13 22 13017.2 13014.5 
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culated. These origins, however, corresponded to 
nothing in the known absorption vibrational 
analysis but agreed quite well with what the 
absorption heads would be if we increased the 
numbering of the initial state vibrational levels 
by three. These doublets are listed in a manner 
similar to the preceding table in Table II1,—the 
vibrational assignments in the first two columns 
being the new assignments required to account 
for the magnetic rotation doublets. 


VIBRATIONAL ANALYSIS 


The suggested change in the vibrational assign- 
ments demanded a reinspection of the absorption 
bands. A change in the vibrational numbering 
was not altogether unexpected for it had been 
shown by Loomis’ that the formula for the final 
state vibrational levels, which he was able to 
determine very accurately from the magnetic 
rotation data, of the blue-green bands, failed to 
represent the red bands in their reported posi- 
tion.® 

The region beyond A6800A was _ therefore 
rephotographed. The new measurements are 
listed in the next to the last column in Table II 
and in Table III. The band heads agree quite 
- well with Loomis’ formula for the final state and 
the entire system can be represented by the 
formula 


v= 14,680.4+117.6(v’ +3) —0.38(v’ +3)? 
+3) +0.726(v"’ +4)? 
+0.0027 +3)?; 


the last three terms are the recent values. given 
by Loomis and Nusbaum"? for the final state of 
the green bands. The last column in Table II 
represents the wave numbers of the bands as 
calculated from the above formula. The last 
three columns should agree in value and it is 
observed that the agreement is fairly good. 


CONCLUSION 


In Fig. 2 we have the usual square array for 
representing vibrational data. The circles repre- 


§ Loomis, Phys. Rev. 31, 328 (1928). 
® Fredrickson and Watson, Phys. Rev. 30, 431 (1927). 
%” Loomis and Nusbaum, Phys. Rev. 40, 383 (1932). 
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sent the bands in which the origins could be 
calculated from the magnetic rotation doublets 
and the crosses represent the measured absorp- 
tion heads. The magnetic rotation doublets have 
not been recorded beyond A7400A because 
beyond that region we have only low dispersion 
plates and the magnetic rotation lines on these 
plates are very broad,—each line undoubtedly 
made up of many lines belonging to a number of 
bands. 

At the left in Fig. 2 are listed the K’ values of 
the various rotational levels of the various upper 
vibrational states which give rise to magnetic 
rotation doublets. It is interesting to see how the 
numbering seems to repeat itself at every fifth 
vibrational state. 

That just these levels should be the ones which 
give rise to magnetic rotation doublets is strong 
evidence that these levels are perturbed. We have 
previously mentioned that the red bands, being 
a 'Y—'Z system, should show no magnetic rota- 
tion because there exists in neither state a mag- 
netic moment along the internuclear axis. It has 
been considered possible, however, that the 
molecule may develop a magnetic moment at 
high rotational quantum numbers'' because of 
the uncoupling of the orbital angular momentum, 
but if such were the case we would expect all high 
levels to cause rotation. The presence of only an 
isolated doublet in a band practically limits the 
interpretation to that of perturbation between 
the upper 'S state and some state which has a 
magnetic moment. 

We may speculate as to what that perturbing 
state might be. A II state is predicted to lie 
below the upper 'S state'? and recently Van 
Vieck® has shown that perturbations between 
such states are possible. If the *II, or *II, are 
the perturbing levels we could account for the 
magnetic rotation by the presence of the mag- 
netic moment in these states. If the perturbing 
level is *IIp we would expect no magnetic rotation 
unless the state had developed a magnetic 
moment at high rotation due to the uncoupling 


" Jevons, Report on Band Spectra of Diatomic Molecules, 
p. 260. 

12 Mulliken, Rev. Mod. Phys. 4, 16 (1932). 

Van Vleck, presented at the Washington Meeting of 
the Am. Phys. Soc., April, 1933. Phys. Rev. 43, 1047A 
(1933). 
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of the angular and spin momenta as in case > or 
in Paschen-Back effect. It is found that the 
observed doublets do have large K’ values which 
js in agreement with the preceding requirement. 

The fact that the K’ value seems to repeat 
itself at every fifth state suggests that the per- 
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turbing state closely parallels the 'S state and 
has very nearly the same r, and w, values. 

The authors wish to thank Dr. W. W. Watson 
for his helpful discussions and also the Rumford 
Fund Committee for a grant which made possible 
the construction of the solenoid. 
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The Spectrum of Sulphur Dioxide 


TuncG-CuinG Cuow,* Palmer Physical Laboratory, Princeton University 
(Received August 2, 1933) 


A study has been made of the emission and absorption 
spectrum of sulphur dioxide with special emphasis on the 
region from 2000 to 2600A. It was found possible to run a 
discharge in flowing SO, with comparatively little dissocia- 
tion. The work on absorption was supplementary to that on 
emission and carried out at the same moderate dispersion. 
A large proportion of both the emission and absorption 


bands found fit into a scheme of vibrational levels. The 
upper levels indicate a frequency of about 375 cm™ or an 
integral multiple of that value. The lower levels confirm the 
1150 and 520 frequencies previously observed but leave the 
third frequency in doubt. The x, for both the 1150 and 
520 frequencies are small, of the order of one cm™ for the 
latter as determined from a progression of sixteen bands, 


WO years ago a systematic study of emission 
spectra of polyatomic gases was begun in 
this laboratory by Smyth’s' work on carbon 
dioxide. The present author undertook a similar 
study of sulphur dioxide but found it desirable to 
extend his original plan to include the investi- 
gation of the absorption spectrum as well as the 
emission spectrum. 

There was already a good deal of information 
available about the spectrum of sulphur dioxide. 
The absorption spectrum had been first studied 
in the region between 3882 and A2200 by Miss 
Lowater.” Henri’ later reported three systems of 
absorption bands; one weak system from 3900A 
to 3400, one moderately strong system from 3370 
to 2450 and one very strong system from 2350 to 
2000 and beyond. There was no evidence that 
these systems corresponded to three electronic 
transitions except the fact that the three portions 
were greatly different in absorbing power. 
Watson and Parker‘ measured the bands in the 
region from \3880 to 2630. They picked out 
bands (about 20 percent of the total) from this 
region and arranged them into progressions, using 
the fundamental frequencies 1369 (v3), 1164 (»;) 
and 610 (ve) cm™ for the lower state and 341, 387 


* Now at Peking University, Peiping, China. 

1H. D. Smyth, Phys. Rev. 38, 2000 (1931) and the 
literature given there. 

?F. Lowater, Astrophys. J. 23, 324, 338 (1906); 31, 311 
(1910). 

3V. Henri, Nature 125, 275 (1930). 

*W. W. Watson and Allan E. Parker, Phys. Rev. 37, 
1484 (1931). 


and 290 for the excited state. In view of the new 
interpretation of the infrared spectrum the 610 
cm frequency is probably incorrect. Nor do the 
frequencies of the upper state seem probable. 

The study of the emission spectrum was also 
started by Lowater.* Later studies made by 
Johnson and Cameron,® by Henri and Wolff,® by 
Martin’ and by Bair® show two different spectra. 
The one most easily obtained and most ex- 
tensively studied is definitely attributable to SO. 
Another set of bands, degraded toward the 
violet were found by Johnson and Cameron, and 
by Bair. These are apparently from SO, and in 
fact can be very simply correlated with some of 
Watson and Parker’s absorption bands as was 
pointed out by Chow and Smyth.® Other bands 
were found by Bair scattered over the region 
between 2200 and 3200A and for the most part 
degraded toward the red. 

In the present work measurements of the SO; 
bands, particularly those of shorter wave-length, 
have been repeated and extended both in 
emission and absorption. An attempt has then 
been made to arrange the observed wave- 
numbers in a rational scheme of vibrational 
levels. 


5 R. C. Johnson and W. H. B. Cameron, Proc. Roy. Soc. 
A106, 195 (1924). 

*V. Henri and F. Wolff, J. de Physique et le Radium 10, 
81 (1921). 

7 Emmett V. Martin, Phys. Rev. 41, 167 (1932). 

5 Bair, Astrophys. J. 52, 301 (1920). 

*T. C. Chow and H. D. Smyth, Phys. Rev. 38, 838 (L) 
(1931). 
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SPECTRUM OF 


EXPERIMENTAL 


The first thing attempted was the reproduction 
of the emission bands observed by Johnson and 
Cameron between 3105 and 2866A. This was 
unsuccessful but a great number of bands was 
observed below 2300 and attention was directed 
to them. 

For studying polyatomic spectra the first 
essential is a constant flow of gas. In the present 
experiment, sulphur dioxide was let into the 
discharge tube from the reservoir through an 
artificial leak and was drawn out by condensing 
it in a liquid-air cooled trap on the other side of 
the discharge tube. A three-stage Gaede diffusion 
pump was also operating to take care of the 
oxygen resulting from dissociation in the dis- 
charge tube. The amount of dissociation and loss 
of SOz resulting from passage through the dis- 
charge was found to be less than 4 percent. The 
sulphur dioxide was from an ordinary tank of 
SO2 but was distilled under vacuum several 
times before use and was stored over P2Os. 

After trying a number of types of discharge it 
was concluded that a weak glow discharge in a 
tube, such as is shown in Fig. 1, was most 
suitable. This was run at a pressure of about 0.5 
mm with a voltage of 920-1100 and a current of 3 
to 5 m.a. When such a discharge was running 
properly there was no sulphur deposit detectable 
on any part of the walls of the tube after a run of 
twenty hours. It was the positive column that 
was photographed. This filled the whole narrow 
glass tube with a bluish violet color. The quartz 
window was put as near as possible to the source 
to avoid reabsorption. 


gas inlet 
gas outlet 
LS mm 
l2cm 
-) 
at quartz 
mm window 
i “ = = 
= 


Aluminum anode 


graded seal 


—— gas ovtlet 


Fic. 1. Discharge tube. 
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The source proved too weak to give satis- 
factory plates with a Hilger E 1 spectrograph and 
consequently all our measurements are from 
plates taken with the medium size quartz Hilger, 
E 315. The spectra on these plates can be gener- 
ally described as consisting of the SO bands 
(4000-2500A) studied by Henri and Wolff® and 
by Martin,’ and a set of bands of entirely differ- 
ent character in the region 2600-2040A (38,000— 
49,000 cm~'). Fortunately the main parts of the 
two sets do not overlap. The measurements of 
band heads and of some of the maxima of the 
bands (where the heads were not sharp) are 
assembled in the table below. The discussion of 
these bands will be given later. 


Tue ABSORPTION SPECTRUM 


In this short wave-length region the only 
absorption measurements available were by 
Garret'® and were not sufficiently accurate to 
compare with our emission data. Therefore a 
study of the absorption spectrum was made. The 
continuous spectrum of a hydrogen discharge 
was sent through an absorption cell 23.5 cm long 
filled with SO, at a pressure of from 0.08 to 10 
mm. This was at room temperature. Some plates 
were also taken at 240°C with a 7.8 cm quartz 
tube in an electric furnace as an absorption cell. 
The band heads observed are listed in Table Il. A 
more thorough study of the absorption spectrum 
by using a bigger instrument was not carried out 
mainly because it was known that two other 
investigators were working on this problem. 
However, the plates obtained were very useful for 
comparison with those showing emission. In 
each case various checks led to the conclusion 
that the accuracy of the measurements is about 
5 or better. 

Before attempting an analysis of the emission 
bands it is necessary to establish their origin. 
They might conceivably arise from SOs, SO, 
SO*, Oo, Og*, Se or Set. Each of these possibilities 
has been considered in detail and all but the first 
discarded as unlikely. The strongest evidence for 
the assignment of the bands to SOs¢ is the absence 
of fine structure (at the dispersion used) and the 
many coincidences between the emission and 
absorption wave numbers. It is believed safe to 


” Scott Garret, Phil. Mag. 31, 505 (1906). 


TABLE I. Emission and absorption bands in sulphur dioxide. 
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Absorption Transition Emission Transition |Absorption Transition Emission Transition 
50302 Ag 45706 
50248 C; 45686 
49910 456072 bsCy 45665 
49550 Ay 45588 bBo, 
49477 Cs 45572 45560 
49133 B, 45500 biA3 
48781 A, 48789 A, 45455 boAy, 
48724 C3; 45414 
48375 B; 48387 B;? 45393 
482609 biAg 45367 
48355 45351 45350 
48121 45333 a,Co, a2Be 
48083 45314 b3By, byBs 
48005 A; 48009 A; 45199 a, B, 
48000 45188 
47965 C2 45119 5B 
47922 bsC4 45069 
47889 aA 5 45051 b; Be 
47776 45030 aA, 
47742 47747 44987 
47 707 44978 a, By 44977 a, Bo 
47636 aA, 47631 ayA 4 44953 2, bg As 
47623 Bz 44833 aA, 
47576 a,C3, b3B, 47586 b3B, 44798 
47489 b,A 3) 44063 
47391 44603 a 1A 0 
47345 aA, 44570 a,B,, 44573 a2B,, bsBy 
47246 A 2 2A 5 44543 bg Bo, 
47232 a,D3 47230 A 2?, a,B;, 4 44457 6 
47223 a,B;, b3A 4 44437 
47198 44362 
47110 b, Bo, aA, 44340 
47074 44305 44310 by, By 
47052 aC; 44293 1 44285 1 
47031 44236 
40968 b.A % 44199 ly b3A 0 
46874 B, 46871 B, 44121 
46865 40851 a,A 44100 
46829 a2B, 44085 
46818 a;Co, b3B; 44057 
46704 a B; 46710 4y a, B; 44039 
46697 43929 43925 b;A, 
46632 5B, 43918 
46588 aA 4 43824 
46570 43807 a2By?, 
46482 Co, eA 4 437 14 
46472 40474 a,By, yA be A & 43683 byA 
40464 43057 
40444 43550 
46425 43538 
46348 46350 43453 
46330 43407 43406 be A, 
46277 biA4 43340 
46259 43329 
46217 43310 by3B;, 
46190 byA 43145 43140 
46118 By 43137 
46104 a,A 2 46105 aA 2, 3A 6 43054 a3A 1 
46084 40094 43034 be Bo 
40077 a2B;, b;B, 40074 b3By 43018 
40060 42966 
45968 s 42954 
45923 42938 71 B, 
45888 42895 byA,? 
45859 42847 
45833 42759 
45817 420681 
45757 Ao, 45757 Ao 42058 a3Bo, 
45725 a,B, 45717 a, By, a4Ag 420634 
42589 
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TABLE I. (Continued.) 


Emission Transition Emission Transition Emission Transition Emission Transition 
42531 410600 40735 39412 
42525 byoBe 41573 40644 39372 
42513 41544 40598 39351 
42366 72Be 41 530 bgBo, 40534 39335 
42306 3A 0) bisBs 41422 40522 by Bo, 39224* 

42271 a,yB, 41398 40482 39137* 
42251 41380 40372 39087 
42171 7iBo 41346 40330 bipBs 39074 
42156 41313 bizBs 40278 39061 
42103 41231 40139 38977 yaBy 
42079 41108 40100 38956 
41994 y2A2 41089 40008 38919 
41923 41038 39973 38907 1 
41902 bg A, 41022 bBo 39930 38847 
41855 40939 39855* 38723* 
41817 40888 39825 38611" 
41802 40859 72Bo 39755 38534 
41723 40841 39730 38232 
41701 40822 bisBs 39541 bieBo 38213 
41661 40767 39522 bisBo 
41621 2B, 40743 39472 
*SO band. 


conclude that the bands under consideration 
arise from the sulphur dioxide molecule both in 
emission and absorption. 

The wave numbers given in Table I are in 
most cases the violet edge of bands degrading 
toward the red. In the cases of a few very diffuse 
bands the measurements are of intensity maxima. 
The first column gives the wave numbers meas- 
ured in absorption, the third those in emission. 
The second and fourth columns give the sug- 
gested transitions in terms of the levels of the 
upper and lower states shown in Fig. 3. These 
transitions are recorded only where the agree- 
ment between observed and calculated values is 
better than 10 cm~'. Nor are transitions recorded 
where the agreement is good but the absence of 


49477C, 


491328, 
49550A, 


_ > 48724C, 
_> 487804, 


|_| 


other bands in the progression indicates that it is 
accidental. The levels used are discussed below. 


ANALYsIS OF ABSORPTION DATA 


The most striking features of the absorption 
spectrum are seen more clearly in the densitome- 
ter records of the plates than in the plates them- 
selves. Such a record is reproduced in Fig. 2. The 
most natural assumption seems to be that the 
strong bands or groups of bands standing out so 
clearly on the short wave-length side of the 
spectrum correspond to transitions from the 
normal state of the molecule to various vibra- 
tional states of an excited electronic level. There 
are nine or ten of these groups with separations of 
about 370 cm~'. Closer inspection shows that 


Complete Absorption 


50298C, 


& 
a 
= 


| | | No Absorption 


Fic. 2. Densitometer record of part of plate showing absorption spectrum of SO, at a tem- 
perature of 240°C, pressure of 1.5 mm and path length of 7.8 em, 


641 
| 


642 TUNG-CHING 


these groups do not constitute a simple v’ 
progression as both their structures and sepa- 
rations vary in a complicated way. Various level 
schemes give a fairly satisfactory correlation of 
the band heads but none has been found that 
seems so much better than the others as to be 
established with certainty as the correct one. In 
any case the various level schemes devised have 
certain features in common. For the upper levels 
the vibrational frequencies used are of the order 
of 375 cm™ or integral multiples of that number. 
The separation of successive levels seems to 
increase for the first two or three and then to 
decrease again. This effect can be seen clearly in 
the upper level scheme of Fig. 3 which is the most 
satisfactory on the whole. On this scheme one 
vibration frequency is about 750 to 770 cm"! 


50,000}- Cs Bs 
“As 
8 
49,000 Ay 
3 Bs 
“As 
46, 
47,0 ' 
46,000} Be 
Ay 
45,000}- 
sooo. y 
1597_,, 
is 
7097 
7000+ 
6593 x 
6101 b 
6000+ 
5600 
5100 
5000} Dio 
4600 4595 
4093 
3450 a, b 
A2500 A2205 
2000 2069 ae 
AISS3 Al6TO 
OoL_A 


Fic. 3. Vibrational levels in SO. In upper electronic 
state C levels were not used in analyzing emission spectrum. 
In lower (normal) electronic state only levels marked A 
were used in analyzing absorption spectrum. 
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with the others uncertain, though the constant 
differences B,—A,—-,=1129 suggest this as an- 
other. An interesting point is the regular increase 
of the A,C,_; separation, the successive values 
being 23, 40, 57, 73 and 54 cm”! giving a second 
difference of 17 cm~ if the last value be omitted. 
However the C levels are necessary for only a few 
bands and most of their combinations with the 
lower levels are definitely missing or coincide 
with other bands. They are perhaps not real. 

A single lower level and the triple system of 
upper levels given in Fig. 3 explain the principle 
features of the strong absorption shown in the 
densitometer record. But if this record is ex- 
amined more carefully or the table of measured 
wave-lengths consulted it is clear that there are 
many bands that do not fit into such a scheme. 
Probably there are many more vibration levels in 
the upper state than have been used but it is also 
probable that other levels than that of zero 
vibration in the normal state contribute to the 
absorption. From the Raman!" and infrared® ® 
spectra we know that the vibration frequencies of 
the normal state are v;=1150, ve=525 and »; 
=1360 (or 13402). An examination of the 
measurements shows 1150 cm”! differences fre- 
quent but 525 and 1360 differences only oc- 
casional. Levels corresponding to 7," =1 and 2 at 
1150 and 2300 are well established and combi- 
nation levels at 1674 for v,)’=1, ve’’=1 and at 
2210 for = 1, v2’’=2 seem probable. A level at 
1555 wve’’=3 also gives several bands and is 
included in spite of the absence of v2’ =1 and 2. 

The intensities of the bands from the higher 
vibrational states are lower than those from the 
initial state as they should be, but there is no 
consistent variation with temperature between 
the plates taken at room temperature and 240°. 
However these experiments were not sufficiently 
extensive to be of much value. 

The predominance of vibrations involving »; in 
the lower state remains unexplained, particularly 
as ve and perhaps v3 states were observed in 
emission. It suggests the possibility that the 
observed 1150 differences should be attributed to 


uC, R. Bailey and A. B. D. Cassie, Proc. Roy. Soc. 
London, A137, 622 (1932); C. R. Bailey, A. B. D. Cassie 
and W. R. Angus, Proc. Roy. Soc. London, A130, 142 (1930). 
12K. W. F. Kohlrausch, The Smekel-Raman Effect, p. 172. 
8 Dadieu and Kohlrausch, Phys. Zeits. 33, 165 (1932). 
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the upper state, the 1128 difference already 
observed there being the interval between the 
second and third levels of this particular vi- 
bration. 


ANALYsIS OF EMISSION DATA 


Taking the same upper levels and lower 
frequencies used for absorption it is possible to 
put some 70 percent of the emission bands into a 
level scheme. As is to be expected many new 
levels appear in the lower state. These are shown 
in Fig. 3. The most striking features are long v2” 
progressions. In one case there are transitions 
from the By. upper level to 17 v2’’ levels having 
successive differences between 520 and 493 cm“. 
Other v2”’ progressions are interpreted as tran- 
sitions from Bo, A;, A; and B; to the same set of 
lower levels. 

The v,’’=1, 2 levels appear as in absorption. A 
levei at v,’’ =3 gives a number of bands but most 
of them correspond to other transitions also. The 
2;"’=4 level falls on the v2”’ = 9 level. 

A search was made for v3"" progressions and the 
levels i, C2, C3, were established tentatively as 
v;’=1, 2, 3. The same differences among bands 
of lower wave number led to the y;, 72, y3 levels 
interpreted as v’’=5, v3’"=1, 2, 3. However, all 
but one of the bands given by the c levels have 
other interpretations so that no ¢ transitions are 
recorded in the table. This is not the case with the 
7 levels but their interpretation is very dubious. 

Turning to other combinations, transitions to 
the as, i.e., 2 levels of ab- 
sorption are also present in emission. An ad- 
ditional level 6;=2502 is quite well established 
and is apparently v;+ 73. 

It may be of interest to note that the second 
differences in both the 2,’ and v2" progressions 
are very small. A similar result has been observed 
in COs. 

FLUORESCENCE 

Just as the number of lower levels involved is 

restricted in an absorption spectrum the number 
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of upper levels observed is restricted in a 
fluorescence spectrum. Mr. Lotmar in Géttingen 
has kindly communicated to us his preliminary 
results on fluorescence in SOs. They seem to be 
reconcilable with the level scheme suggested 
above but not an emphatic confirmation of it. A 
detailed discussion must await publication of 
Lotmar’s work."* 


THE ELECTRONIC STATES OF SOs 


In the early part of this paper the possibility 
was mentioned that there was more than one 
electronic transition involved in these spectra. 
The present work deals mainly with the shortest 
wave-length portion of the spectrum. Absorption 
here is stronger than that at longer wave-length 
and is separated from it by a gap. Furthermore no 
emission bands were found on our plates above 
2600A. It seems probable that there are two 
excited states, one having a », of about 46,000 and 
the other having a », of about that used by 
Watson and Parker, i.e., 31,500 cm. The 
publication of the data on fluorescence may 
throw further light on the question. 

The ultimate purpose of the experiments 
described in this paper was to throw light on the 
structure of SO. molecule in normal and excited 
states, its normal vibrations, translation proba- 
bilities, etc. But the results seem hardly to justify 
a discussion of these questions as yet. It is hoped, 
however, that they may be of assistance to others 
working in this field. 

In conclusion the writer wishes to express his 
deep gratitude to Professor H. D. Smyth who 
directed this work and helped to prepare this 
article. He also wants to thank Professor R. 
Ladenburg and Professor A. G. Shenstone for 
their help. 


4 Note added in proof: Lotmar’s paper appeared in the 
Zeits. f. Physik 83, 765 (1933) and is discussed in a letter 
by Smyth in this issue of the Physical Review. 


OCTOBER 15, 1933 


PHYSICAL REVIEW 


VOLUME 44 


On the Inversion of Doublets in Alkali-Like Spectra 
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Many of the doublets in the spectra of alkalies and alkali- 
like ions are inverted or anomalously narrow. This ‘“‘ten- 
dency to inversion"’ may be explained by the polarizability 
of the core, which may be formulated as the repulsion by 
the excited core configuration states of the levels of the 
optical doublet. A study is made of the energies and wave 
functions of these states for the spectra Na I to S VI, on 
two coupling schemes which allow approximate computa- 
tions to be made. On the assumption of LS coupling the 
separations of the excited core doublets turn out to be of 
the order of magnitude of the 2p x-ray doublet; half the 
excited core doublets are inverted and half of them normal. 
The two sets of doublets repel the levels of the optical 
doublet unequally, and, in the series of spectra considered, 
the differential repulsion is of the right sign to give an in- 
verting tendency of the same order, though somewhat less, 


than the observed inversion. The same problem is then 
treated by another approximation, by assuming the 2 
x-ray doublet energy large compared to the electrostatic 
interaction energy of the 3p and 3d electrons. The results 
with both models are proportional to the spin orbit inter- 
action energy of the “hole” in the 2p shell. The condition 
that the effects studied should tend to invert and not to 
broaden the optical doublets, and the further conditions 
that they should produce an actual inversion are considered, 
In most cases these questions may be answered from a con- 
sideration of the quantum numbers associated with the 
states capable of perturbing the doublet, and the ratio of 
optical and x-ray doublet separations. In this way we car 
understand the observed incidence, both of the tendency 
toward inversion and of the inversion itself. 


HE ordinary spin doublet formula may be 

derived from a consideration of the inter- 
action of the magnetic moment of the electron 
and the magnetic field, H=1/clEXv_], produced 
by its velocity relative to the nucleus. For a 
Coulomb field the energy of this interaction is 
given by 


Rheo®Z* j (7 +1) —1(1+1) —s(s4+1) 
2 
The doublet separation in wave numbers is 


with the level of greater j value lying higher. This 
formula is applied to the spectra of alkalies and 
alkali-like ions, with an effective nuclear charge 
substituted for Z to take account of the screening 
by the core electrons. The formula often fails, 
however, and the occurrence of inverted or 
anomalously narrow doublets is frequent in 
alkali spectra.' Thus the 5f doublet of Cs I is 
inverted, with Avy= —0.16; the 4f doublet has 
not been resolved. The d and f doublets of Rb I 
are believed to be inverted.? We quote in Table 


1 Bacher and Goudsmit, Atomic Energy States, 1932. 
2? Ramb, Ann. d. Physik 10, 311 (1931). 


I the 3d separations in frequency numbers for 
stripped atoms from magnesium through sulphur, 
as observed by Bowen and Millikan.’ The 
theoretical separations are calculated of the 
assumption of perfect screening, with 


Zess=(Z—10). 


TABLE I. 3d and 4d separations. The observed values are from 
Bowen and Millikan. 


3d 4d 
observed from theory observed from theory 
Mg Il —0.91 0.576 0.243 
Al Ill —1.70 2.92 —1.28 1.23 
Si lV 1.57 9.2 0.08 3.9 
PV 6.32 22.5 4.41 9.5 
SVI 37.75 46.6 19.7 


For an electron moving in a central static field 
the doublet separation is 


W 


where £, is the radial component of the electric 
intensity, and the bar denotes an average over 
the stationary state of the electron. Now for a 
neutral atom, and for all 7, E, is positive; this 
follows from Gauss’ theorem, since the charge 


§ Bowen and Millikan, Phys. Rev. 25, 301 (1925). 
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within any spherical surface is necessarily 
positive. Thus the model of an electron moving 
in any central static field necessarily gives a 
normal doublet, with the level of greater 7 lying 
higher. 

Two distinct quantum mechanical aspects of 
the core are neglected in the one electron 
model. The first of these arises from the circum- 
stance that the interaction of the valence electron 
and the electrons of the core is not completely 
represented by the field of the mean charge dis- 
tribution of the core electrons; the application 
of the exclusion principle to the formulations of 
these interactions introduces, as is well known,’ 
terms which correspond to the interchange of 
core and valence electrons. These terms give 
contributions to the energy of the two states of 
the optical doublet which are in general unequal. 
Their effect has recently been investigated by 
Johnson and Breit,® but is apparently too small 
to account for the observed inversions. 

The second, which seems to us to give greater 
promise of an explanation, involves a consider- 
ation of the excited states of the core. This aspect 
may be stated in various ways: as a polarization 
of the core, expressed by transitions of the inner 
electrons to excited states; as an admixture of 
wave functions for the excited states in that for 
the true stationary state described approximately 
by the definite configuration of closed shells and 
the quantum numbers of the valence electron; or 
as a repulsion of the levels of the ordinary optical 
spectrum by those of the excited core con- 
figurations. These are three statements of the 
same effect, but one quite different from that 
mentioned in the preceding paragraph. 

Strictly a level is said to be ‘“‘perturbed”’ by 
another of the same parity (oddness or evenness 
of |/|) and the same angular momentum, which 
lies close to it so that the nondiagonal element of 
the energy referring to the two states, in general 
small, is no longer negligible compared with the 
unperturbed energy difference. An example of a 
perturbation of this sort which anomalously 
inverts the doublets of certain terms is known 
in the spectrum of copper, in the configuration 
(3d)'° (np). Shenstone and Russell® have shown 

‘Fock, Zeits. f. Physik 81, 195 (1933). 


* Johnson and Breit, Phys. Rev. 44, 77 (1933). 
‘Shenstone and Russell, Phys. Rev. 39, 415 (1932). 


that the perturbing levels in this case are the 
inverted doublets of (3d)* (4s) (mp), and have 
given on this basis a complete explanation of the 
inversions observed. A similar explanation has 
been suggested at various times for the inverted 
doublets of alkalies.’ This case cannot really be 
an example of spectroscopic perturbation: the 
excited core states lie far from the optical ones, 
and the interaction energy is at best only a small 
fraction of this energy difference. The separations 
of states of the excited core configurations, how- 
ever, are relatively enormous, arising partly 
from a deep lying unpaired spin, and in some of 
the states the j values will be inverted, analogous 
to those of x-ray doublets. It might be expected, 
then, that the existence of these higher states can 
totally alter the doublets of the optical terms, 
and that their arrangement may be such as to 
give a net inverting effect. 

The atomic wave function for a true stationary 
state of an atom may be built up as a sum of 
wave functions for the approximately stationary 
states of definite configuration, and equal parity 
and angular momentum. In general nearly the 
whole of the true wave function is given by 
the term arising from one specific configuration, 
and this configuration is used (in addition to the 
quantum numbers m ,m,) to index the state; the 
coefficients in the sum of wave functions referring 
to other configurations dre small unless the 
conditions for a genuine perturbation are ful- 
filled. Thus we may write 


where y;,-is the approximate function for a state 
of definite configuration i and Vj;; is the matrix 
element of the electrostatic interaction energy 
corresponding to a transition between the con- 
figurations 7 and j. The two terms 7 and j repel 
each other by an amount 


Vi;!*/| 


Suppose that y,;, represents a doublet of the 
same orbital angular momentum L as the optical 
doublet, and neglect the normal doublet splitting 
in comparison with that of the excited one, e. 


7H. E. White, Phys. Rev. 40, 316 (1932); R. F. Bacher, 
Phys. Rev. 43, 269 (1933). 
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The shift in doublet separation is given, to first 
order in ¢ and v, by 


6E=AE 14 | Vi; | 2¢ /E,y?+2 Vi;v/ Eo. 


For definiteness we shall use this theory in an 
attempt to explain the data given in Table I, 
and in particular the inversion of the 3d doublet 
(j=3/2, 5/2) of Mg II, observed to have a 
separation of —0.9v instead of 0.56» as pre- 
dicted by tke spin theory. The normal Mg II 
atom consists of a single valence electron outside 
a completed 2 shell, so that a full description of 
the 3d configuration is (1s)? (2s)? (2p)® (3d). We 
shall consider only transitions of the p electrons, 
and need not include (1s)? (2s)? in designating 
the configurations. The lowest even levels with 
j=3/2, 5/2 arising from core excitation are in 
(2p)® (3s) (3p), but this involves a double jump 
from (2p)* (3d), so that a much more consid- 
erable contribution may be expected from 
(2p') (3p) (3d). Higher terms, with transitions 
of one of the 2p electrons to 4p, 5p, etc., will 
have similar but smaller effects. Unfortunately 
these states cannot be observed experimentally; 
the energy required to excite a core electron is 
much greater than the ionization potential of 
the valence electron, so that auto-ionization 
occurs. 

The number of possible combinations of the 
individual / and s vectors of an almost closed 
shell is equal to that for the missing electrons. 
Thus in the elements of Table I the number of 
levels arising from (2p)° (3p) (3d) is the same as 
that which would arise from (2p) (3p) (3d), and 
we are justified, insofar as multiplicities are 
concerned, in treating the atom as a _ three 
electron atom, 2p representing a “hole’’ in the 
2p shell. Two different coupling schemes for the 
excited core will be considered as limiting cases, 
although, since the electrostatic interactions and 
the 2p spin orbit interactions are of the same 
order of magnitude, the actual atom will be of 
some intermediate type. Calculations will be 
made first on the assumption of LS coupling, 
then according to a scheme in which the valence 
and excited electrons are coupled LS and then 
to the quantized j of the 2p hole. In the latter 
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case the levels fall into two groups, corresponding 
to the two members of the inverted x-ray 
doublet. It will be seen that the repulsion of both 
components of the optical doublet is the same jf 
only the direct integrals of V\; are considered, 
i.e., those which would appear even if the wave 
functions were not antisymmetrized. If exchange 
terms are included, however, there is a net effect, 
proportional on both models to the spin orbit 
interaction energy of the 2p hole. 


Errect or LS Excirep Core STAteEs 


Since the electrostatic energy commutes with 
Land S it has no matrix elements connected with 
transitions in which L and S change. Thus the 
whole effect of the excited core states, if LS 
coupled, comes from the 7)’s of the excited 
configuration. An enumeration of these relevant 
states may be made from the vector model. 
First couple the 3p, 3d electrons to form singlet 
and triplet P, D and F terms. One ?D (L=2) 
arises from each singlet and triplet of the two 
electron configuration, making six in all. There 
are no equivalent electrons and all states are 
allowed. 


TABLE II. Diagonal elements of + f2(le-s2). 


Parent (1, +51) + Se) 

term J=s/2 J =3/2 Separation 


The separations ¢ of the individual doublets 
are given by the matrix of the spin orbit inter- 
action energy. This matrix is not diagonal in the 
LS scheme but we may assume the diagonal 
elements large and calculate them for the *D’s 
formed by the addition of 2p to 'PDF, *PDF. 
A method for obtaining the spin orbit inter- 
actions in LS coupling directly from the matrix 
equations has been given by Johnson.* Table 
II contains the diagonal elements of 51) 


+ where 


*M. H. Johnson, Phys. Rev. 38, 1628 (1931). 


INVERSION OF 


for 2p, and & is the corresponding parameter 
associated with the ionic multiplet. It may be 
assumed that these parameters for the ionic 
multiplets are negligibly small compared with 
t. As has been pointed out by Shortley® the 
spin orbit parameter of a missing electron in an 
otherwise closed shell is negative. Those doublets 
arising from 'P, 'D, *F are thus seen to be in- 
verted, while the other three are normal. The 
orientation of the doublets with P and F parent 
terms could have been predicted from the vector 
model. 

The linear combinations of determinantal 
eigenfunctions to represent the various D 
doublets were obtained by using angular mo- 
mentum operators, a method given by Gray and 
Wills."° The valence and excited electrons are 
first combined, then coupled with 2p, so that 
correlations with the individual doublets shown 
by the vector model may be made. No constant 
factor arises in the V,; because of the presence 
of six electrons in the p shell. Specifying the 
quantum numbers of 3p (of which there are six 
possibilities) in the excited core wave functions 
automatically selects terms in the antisym- 
metrized function for the ordinary *?) which 
combine. There are 3! nonvanishing terms for 
each pair of determinantal functions in the 
three electron case, and this factor is cancelled by 
the normalizing factor 1/(3!)!. If the problem is 
done explicitly in all seven electrons the 7! 
terms in V’,; are taken care of by the 1/(7 !)! of the 
antisymmetric permutations. The coefficients of 
the first direct and interchange radial integrals 
of V,;, the matrix element of e?/r:2, where 1 and 
2 refer to 3d and the electron making the 
transition 2p—np, then, are found by multiplying 
the coefficients of determinantal functions which 
combine with the optical doublet by the proper 
angular integrals (tables of which may be con- 
structed from Slater's and summing. This 
result is independent of the m of the 7=5/2, of 
course, and is the same for both members of 
the doublet. Because of the antisymmetry in 
spin of the singlet the sign between the direct 
and exchange integrals is different according to 
whether the parent term of state j is a singlet 


*Shortley, Phys. Rev. 40, 185 (1932) 
” Gray and Wills, Phys. Rev. 38, 248 (1931). 
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or a triplet. To this order, then: 


3\3 1 
10 15 

3 1 

10! 15 


G 

3\! 1 

2 5 
7\3 2 

vat 
21\! 2 

= (—) (#26). 
10 5 


The shift of optical doublet separation may 
be calculated from the formula = | 
substituting the values of the separations from 
Table II. The terms in Fo? cancel, but the further 
terms in V,,? give a net result. Since the exchange 
integral G; is much smaller than the direct 
integral Fy we shall here record only the first 
nonvanishing terms (Table III). Whether this 


Taste III. 


bECP) = FoG/Ee*) bECD) = —'4(FoG/Ee) 

bECP) = bECF) = 5(FoG/Ec*) 

= —2(FoG/Ee*) és. 


result means inversion depends on the relative 
sign of Fy and G;, on their magnitudes, and the 
energy differences of the levels, but before mak- 
ing an application to an actual case we shall 
investigate the other coupling scheme outlined 
above. 


Errect or CoupLep ExciTED CORE STATES 


In this scheme the quantized j of 2 is coupled 
to the ionic multiplets of the outer two electron 
configuration. Since the total orbital angular 
momentum of the three electrons is no longer 
quantized there will be matrix elements Vj; 
between the optical *D and all the J =3/2, 5/2 
states of the excited core, of which there are 
altogether thirty-six. The appropriate angular 
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momentum operators were used to obtain the 
three electron eigenfunctions. Products of coef- 
ficients of combining terms of the excited core 
and optical state wave functions were multiplied 
by the appropriate angular integrals and summed, 
for each of the relevant excited core levels. The 
contributions to the energy of the mutual repul- 
sion are proportional to the squares of these 
sums. Since the levels from a single multiplet 
coupled to a particular j of 2p lie very close 
together, the inverting effect of such a set of 
levels is proportional to | Vij 23/2). 
The results are given in terms of the first direct 
and interchange integrals. 


Jep=3/2: | Vis| 23/2) = — FoGi, 
Jep=1/2: | | 23/2) = 3 FoGi. 


The coupling scheme here assumed implies that 
all states associated with a definite j of 2p are 
grouped together, and their spread of energies is 
neglected. The two sets of levels act in opposite 
ways on the optical doublet separation, but they 
are separated by e, the width of the 2p x-ray 
doublet. The sum of their effects is 


LOE =§ 


since the 23,2 group lies nearer the optical 
doublet. 


Now, since 
e= (3/2), 
— (FoGi/Eo*) 


The net effect in this case is due entirely to 
the difference in the sign between Fy) and G, for 
excited core levels with singlet and triplet parent 
terms, and is equal to that portion of the LS 
result which arises from the same cause. Because 
of the different enumeration and character of 
the states no extra contribution from the triplets 
occurs here. Owing to the quantization of the j 
of the 2p “hole” instead of the total orbital 
angular momentum the electrostatic interaction 
energy associated with the transition 2p—np 
must be calculated separately for the J/= 3/2, 5/2 
states arising from each component of a triplet, 
and the individual elements squared before 
summing. Thus the weighting of terms from 


triplet parents, which appears in LS coupling, is 
lost, and the whole effect comes from the dif- 
ference in symmetry between singlets and 
triplets. 


APPLICATION OF THEORY TO MG IT. EVALUATION 
OF INTEGRALS 


Fy and G;, are the first two radial integrals of 


f fue 


when 1/rj2 is expanded in spherical harmonics. 
In the evaluation of these integrals for a par- 
ticular case, Mg II, hydrogen like one electron 
functions were used, with screening constants 
derived by interpolation from those given by 
Pauling and Sherman." For the configuration 
whose energy lies closest to the optical state, 
(2p) (3p) (3d), the values of Z.:; were 7 for 2p, 
3.7 for 3p and 2 for 3d. The results were found 
to be surprisingly stable under a_ reasonable 
Variation of screening constants. 


Fyo= ef R,(2p)Ri(3p) R2?(3d) (1 dre 


1.73 volts 


where R is the radial factor of the one electron 
eigenfunction. 


X = —0.43 volt electron. 


The radial wave function corresponding to 
spins parallel and antiparallel to the orbital 
angular momentum are not identical. This 
difference is appreciable only for the 2 electron, 
since only there are the energies for the two spin 
orientations considerably separated. Thus |); 
for levels associated with and are 
different. | 12=2 is small but since 
it is associated with only one set of the per- 
turbing levels Vj;9v is weighted by a factor Eo/¢ 
in comparison with | Vij! ?. 

To ascertain a correction for this difference 


we may write 


"Linus Pauling and J. Sherman, Zeits. f. Krist. 81 
(1932). 


is 


if- 
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Fo Fo ve=f, G, 32 8, 
f= e pi(2p) 
0 0 


where pi=Ri(2p)32—Ri(2p)1/2. The relativistic 
radial wave functions for a Coulomb field have 
been given as series in r by Darwin": 


3 
R(2p)3/2= reo(i In a), 


2 


a) 
R(2p)1/2= 1 In 2cr 


a;? 37 or 3 
4 8 2 4er 


where aq, is the fine structure constant times the 
nuclear charge. ¢=Zes;/na. 
25 


der 


To apply this to our problem we must determine 
what effective nuclear charge to use in evaluating 
a;. It is not possible to do this rigorously in any 
simple way; we have contented ourselves with 
elementary considerations which fix the possible 
limits of Z, and give us an estimate of the cor- 
rection term. A mean for the parenthesis was 
estimated before integrating, and found to be 
approximately unity, so that the integral is 
merely a,?Fo/4. Since both 3d and 3p weight the 
difference R(2p)3,2—R(2p)i;2 far out from the 
nucleus the corresponding parenthesis in g is of 
the same sign, and is of the same order of mag- 
nitude. The overlapping coefficients are just as 
before, and the ratio of the correction to the 
effect previously considered is 


Fogt+Gif Ey Eo 


4 


The change in the optical 2D splitting predicted 
by the theory depends on the data assumed and 
on the coupling. «= 2219», the separation of ?P 
in Mg IV. The energy difference between the 
excited core and optical states is that of the 
transition 2-3), approximately 480,000» or 
60 volts, from the spectrum of Mg III. The 


"C. G. Darwin, Proc. Roy. Soc. A118, 673 (1928). 


effective charge which enters the v correction 
through must lie between 7 and 11, and about 
9 seems a reasonable choice for it. This leads 
to 1/3 as the value of the ratio above. On LS 
coupling the shift in the 3d doublet separation is 


bE = (4/3) FoG,/ Eo?) (1+ 4) 0.9», 


sufficient to invert the levels, although not to 
the extent of the 1.5» discrepancy indicated by 
observation. The jJ result is half this, and the 
actual Mg II atom must lie somewhere between. 
Further contributions toward inversion may be 
expected from transitions of a core electron to 
other p orbits. The direct integrals for transitions 
to 4p and 5p were calculated, and found to be 
0.9 and 0.56 volt electrons, respectively. The cor- 
responding exchange integrals are small, so that 
the computations cannot be rigorous. If they 
diminish in the same ratio as the Fy’s the com- 
bined inverting tendency of these configurations 
is about 0.4 times that of the first excited core 
level. We shall see, moreover, that all the 
inverting effects would be somewhat larger if 
more accurate radial wave functions were used. 
We are thus led to believe that the observed 
inversion of the 3d doublet may be wholly 
accounted for by the effect treated in this paper. 


CONCLUSIONS 


The effect of excited core states on the optical 
doublet separation is toward inversion if Fo 
and G, are of opposite sign. A general argument 
for the signs of these radial integrals can be made 
for the elements of Table I. For the nearest 
perturbing state, (2p) (3p) (3d), Fo and G, are 
opposite in sign. Only one wave function, that 
of 3p, has a node. In the direct integral the 
function is weighted by 2p, which has its 
maximum far in; for the exchange integral the 
important region is that of greater r, because of 
the overlapping with 3d. The node will always fall 
between these two weighted regions. While the 
largest contribution may be expected from this 
configuration similar effects arise from transitions 
to other p levels, the magnitude of whose cumu- 
lative shifting tendency it is difficult to estimate. 
More nodes occur in the mp wave functions, 
(n>3), but the principal contribution to G, will 
come between the first two nodes, and that to 


_| 
nd 
of 
ir- 
ts 
yn 
e, 
p, 
id 
le 
2 
n 
ij — 
€ 


650 


Fo before the first one, so that the difference in 
sign persists. In a non-Coulomb central field, 
such as Hartree’s, all the nodes are farther in 
than those of hydrogen-like functions, making 
|—G,| greater. It seems not likely that the 
accompanying diminution in Fy would entirely 
compensate this increase, so that the numerical 
answer should be larger than that calculated for 
the 2p-93>p transition in Mg II. 

For actual inversion there must be, in addition 
to the tendency toward inversion due to different 
signs of Fo and G;, an enormously larger spin 
orbit interaction energy for the x-ray levels 
than for the optical ones. This condition is 
satisfied for the first few elements of Table I. 

In the inversion of 4d the argument is anal- 
ogous. The overlapping for np-4d is in general 
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farther out and opposite in sign to that of 2p- np, 
but the shift, even where the condition on the 
separation ratio is satisfied, should be less than 
that for 3d. This is in accord with the observa- 
tions in Al III. The 5/ term of Cs I is analogous 
to the 4d of Al III. 

Thus it appears that both the tendency 
toward inversion and the inversion itself can be 
roughly predicted from a consideration of the 
quantum numbers associated with the states 
capable of perturbing the doublet, and the ratio 
of the optical and x-ray doublet separations. 

The writer is deeply indebted to Professor J. R. 
Oppenheimer, with whom this investigation was 
carried out, and to the Physics Department of 
the University of California for a Whiting Fel- 
lowship during the past year. 
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An Application of Probabilities to the Counting of Alpha-Particles 


NorMan I. Apams, JR., Sloane Physics Laboratory, Yale University 
(Received August 1, 1933) 


The Poisson law of probability is directly applicable to 
radioactive emissions only in certain simple cases. In 
general, when a radioactive series such as the thorium 
series is present the probability function is complex. It can 
be found readily, however, if the interval of time for which 
probabilities are calculated is chosen so as to be either very 


short or very long compared with the mean life of each 
radioactive substance involved. The probability function 
for a-particles from the thorium series is determined on the 
basis of an interval of five minutes and is illustrated 
graphically in a special case. 


HEN events are distributed individually 

and collectively at random, the proba- 
bility of any number of events occurring within a 
chosen interval is given by the Poisson law 
p(n) =e"e~*/n!.! Here nm is the number of events 
whose probability is required, « is the mean or 
expected number for the given interval and e 
= 2.71828, as usual. The rays emitted during the 
disintegration of a radioactive substance or of a 
mixture of effectively independent radioactive 
substances constitute a set of events distributed 
individually and collectively at random in time, 
at least in the case where the amount of radio- 
active material does not vary appreciably over a 
period of time sufficiently great to contain a 
large number of the basic intervals for which the 
probability of events is desired. Experimental 
verification of the Poisson law as applied to 
simple radioactive disintegrations has been ob- 
tained in the case of a-particles by Rutherford 
and Geiger? and in the case of 6-particles by 
Kovarik.’ 

Unfortunately, in many cases of radioactive 
disintegration the events are not distributed 
individually and collectively at random. Such a 
case has arisen in experiments now being con- 
ducted by Kovarik and the writer to make an 
accurate determination of the disintegration 


‘For an explanation of the criterion “individually and 
collectively at random" and a complete discussion of the 
Poisson law the reader is referred to Fry, Probability and 
Its Engineering Uses, Chapter VIII. 

* Rutherford and Geiger, Phil. Mag. [6] 20, 698-707 
(1910). 

* Kovarik, Phys. Rev. 13, 272-280 (1919). 


constant of thorium. Here a-particles are counted 
by the same method as in the recent new determi- 
nation of the disintegration constant of uranium. 
Since the a-particles counted do not come from 
thorium alone but from the entire thorium series 
containing six a-particle emitters, each dependent 
on the preceding, the events are not individually 
at random and the Poisson law does not apply. 

The general derivation of the probability law 
for the case just quoted appears to be very 
difficult. We can, however, deduce the form of 
the probability function if a simple restriction is 
imposed. Let us suppose that the interval of 
time for which the probability of m a-particles 
must be found is either very short or very long 
compared with the mean life of each of the six 
a-particle emitters in the thorium series. We 
assume that the thorium and its following 
products are in radioactive equilibrium, so that 
the amount of each substance present is effect- 
ively constant, the mean life of thorium being so 
great. Evidently the number of a-particles 
coming from any product x, say, whose mean life 
is long compared to the chosen interval, is a small 
fraction of the number of atoms of x present or 
produced during the interval and the distribution 
of these a-particles in time is in accord with the 
Poisson law. On the other hand, the number of 
a-particles coming from a product y whose mean 
life is short compared to the interval is equal 
simply to the number of atoms of y produced 
during the interval from the preceding substance 
in the radioactive series. Thus, if y follows x the 


* Kovarik and Adams, Phys. Rev. 40, 718-726 (1932). 
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a-particles are associated in groups of two, the 
distribution of the groups in time following the 
Poisson law. Similarly, if another product z of 
short life follows y in turn, the a-particles are 
associated in groups of three and so on. 

The substances in the thorium series which 
give rise to a-particles are shown in Table I. As 
five minutes is a convenient period for which to 
record counts experimentally, let us use this 
length of time as the basic interval in calculating 
probabilities. Examination of the table shows 


TABLE I, 

Substance Mean life 
Thorium 1.8(10)'° yr. 
Radiothorium 2.74 yr. 
Thorium X 5.25 day 
Thoron 78.6 sec. 
Thorium A 0.21 sec. 
Thorium C 87.4 min. 


that thorium, radiothorium and thorium C emit 
a-particles distributed individually and _ col- 
lectively at random, while thorium X, thoron and 
thorium A give rise to a-particles in groups of 
three similarly distributed. Thus, if « is the mean 
rate of emission for the entire thorium series, the 
probability functions for the a-particles coming 
from the “‘singles’”’ class and from the ‘‘triples”’ 
class are, respectively, 

€/6) "t/5e—e/6 

n;! (n,/3)! 


where n,=0, 1, 2, 3, etc., and n,=0, 3, 6, 9, etc. 
The factor } in the second expression is intro- 
duced to normalize the coordinates in the usual 
manner. 

It is now a simple matter to find the probability 
function for any number of a-particles in the 
chosen interval for the entire series. Consider, 
for example, five. This number may be composed 
of one triple and two singles, or of five singles, the 
total probability function for which is 


ps(2) + p.(5). 


NORMAN I. 
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Similarly, the probability function for six @ 
particles is 


p:(6) p.(0) + p.(3) p.(3) + p.(0) p;(6), 
and that for seven is 
ps(1) + p1(3)- p.(4) + p.(7). 


Evidently the general form of the probability 
function depends on whether the number of 
a-particles is one less than, equal to or one. 
greater than an integral multiple of three. If we 
denote by v a number of a-particles which is an 
integral multiple of three, we have, for the three 
cases just mentioned, 


(v/3—1)!2! 
(v/3—2)!5! 1!(v—4)! 
(€/6)%e/2)'" 
O!(y—1)! 


Ne, 2)3 
(v/3)"0! (v/3—1) 13! 


= 


(v/3—2) '6! 1!(y—3)! 
(€/6)%e = 
(v/3—2)!7! 1!(v—2)! 
+ } 
O!(v+1)! 


We may use these expressions separately if we 
choose, but it is simpler to employ the compos- 
ite function P(v) = p_i(v—1)+ po(v) + p4i(v +1) 


ne . 


ve 
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which applies to the coordinate interval y—1 to 
y+1, inclusive. After some manipulation we find 


vs 


Po) 
(y/3—x) !(3x)! 


3x €/2 

—+1+ ) 

€/2 3x+1 


(y/3—x) !(3x)! 


3x 
x +1+ ). 


3x+1 


Since the coordinates have been normalized, 
the actual probability of a number of a-particles 
within the coordinate interval y—1 to v+1 in- 
clusive is 3P(v) and the probability for an ex- 
tended interval to ve+1, say, is P(v). 
When « is large we may replace P(v) by a con- 
tinuous function P(m) and express the proba- 
bility as the integral /(""P(n)dn in the usual way. 

The probability curve P(m) for «=210 is 
shown in Fig. 1. For the sake of comparison the 


~ 
no 110160 100 200 210 220 230 240 250 


Fic. 1. 


curve p(n) =e"e~*/n! is indicated by the broken 
line. This curve applies to any case where 
a-particles are produced by a single substance or 
by a mixture of effectively independent sub- 
stances, such as the uranium series. It is evident 
that the standard deviation o for P(m) is con- 
siderably greater than for p(m). In fact, in the 
first case the value of o, obtained from the curve, 
is approximately 21, while in the second case it is 
known to be 210!= 14.5. This difference is easily 
detected experimentally. 
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The Exponential Law of Radioactive Disintegration 


ARTHUR E. Ruark, University of Pittsburgh 
(Received June 22, 1933) 


The decay law, N= Ne™*, is based on the assumption 
that the probability w of an atom disintegrating in unit 
time is constant. The question arises whether this as- 
sumption is necessary. It might happen that atoms of 
different ages have different deca:- probabilities; and w 
might depend on age in such a way that exponential decay 
curves would result. A theorem of Wolfke, however, can be 
interpreted in such a manner as to show that w cannot 
depend on the age of the atom. This result is generalized. 


Consider a series whose members A, B, etc., have decay 
probabilities wa, wa, etc. If all members decay expo- 
nentially, if all w’s are smooth and continuous functions of 
their variables and if each w is, at most, a function of the 
age of the atom, the ages of its ancestors and the ages of its 


‘descendants; then the w’s are constant; and so in this more 


general case, the assumption of constant decay probability 
which underlies the wave mechanical theory of disinte- 
gration is necessary as well as sufficient. 


I. THE PuysicAL MEANING OF THE 
Decay CONSTANT 


CHWEIDLER derived the radioactive decay 
law, 


N=Noe™, (1) 


by assuming that in a preparation of N atoms 
the number which disintegrate in time dt is 
Nddt, where \ is independent of the time. The 
quantum theory of disintegration may be said to 
explain this assumption. The wave-mechanical 
picture of nuclei as potential barriers and the 
customary statistical interpretation lead to the 
statement that the probability of decay in unit 
time is constant. Wolfke! has pointed out that the 
experimental facts admit an alternative interpre- 
tation. It is conceivable that an individual atom, 
born at time zero, may possess a mechanism 
whose operation will cause it to decay between ¢ 
and t+dt. Wolfke shows that if this is true, the 
exponential law of disintegration requires that a 
fraction \e~tdt of the atoms born at time zero 
shall decay in the interval dt. According to this 
picture, chance is involved in the assignment of 
lives to the atoms at the time of birth but after 
that the behavior of each atom is determinate 
until it decays, when chance must enter again to 
fix the life of the daughter atom and so on. This 
alternation between determinate and _inde- 
terminate behavior seems unsatisfactory. 


1 Wolfke, Phys. Zeits. 30, 899 (1929), 


The discussion by Wolfke leads to the more 
general question whether radioactive atoms show 
aging effects. In dealing with such a fundamental 
point it is well to rely as much as possible on 
experiment rather than intuition. Let us examine 
the meaning of \, considered as an empirical 
constant, on the basis of decay experiments. In 
any preparation the atoms which decompose in a 
time-element dt are of various ages and it is 
conceivable that atoms of, different ages may 
transform at different rates, in such a way that 
the composite decay curve is exponential. If we 
write w(t—a)dt for the fraction of the atoms born 
at time a and still present at time ¢ which will 
decay between ¢ and ‘+d, the experimental value 
of \ appropriate to the time / is the mean value of 
w for atoms of all ages, decaying in that interval. 
Thus time-constancy of \ does not prove that w 
is constant, although, of course, if w were a 
function of age, we should expect that \ would 
depend on the mode of preparing the specimen, 
so that variations of \, due to this cause, would 
have been discovered. However, careful research’ 
has failed to reveal such variations. Therefore 
the question of interest is, what forms of w are 
consistent with a constant value of \ for all 
preparations? This question may be answered 


2 Rutherford, Wiener Berichte, 120, 303 (1911); Poole, 
Phil. Mag. [6] 27, 714 (1914); Curie, Le Radium 7, 33 
(1910); J. de Physique 10, 328 (1929); Behounek, J. de 
Physique 4, 77 (1923). 


654 


LAW OF RADIOACTIVE DISINTEGRATION 655 


partially by reinterpreting the results of Wolfke. 
We define p(¢—a)dt as the fraction of the atoms 
born between a and a+da which decompose 
between ¢ and ¢+dt. Whereas Wolfke considers 
that the life of each atom is well determined as 
soon as it is born, we make no such assumption. 
It is immaterial which atoms will decompose in 
dt; all we wish to know is p(t—a). Assuming that 
the exponential law is obeyed by the whole 
preparation, Wolfke’s equations show that 
p(t—a) =e and that w is constant. All this 
is on the assumption that / is at most a function 
of t—a. Apparently no experiments have been 
performed to show whether the decay constant 
of an aggregation of radioactive atoms depends 
on the lives of their ancestors and unborn 
descendants. The possibility of dependence on 
the lives of the immediate parents could be 
tested by using a preparation containing RaD, E 
and F. RaE is removed and for a short interval, 
I,, the preparation is allowed to grow another 
crop of Rak. This is quickly removed and is 
allowed to produce RaF, which is removed. A 
small additional crop of RaF is permitted to 
accumulate in the RaE preparation, and is 
removed during a short interval J2. This crop of 
RaF consists of atoms having parents which were 
born in the interval J; and which disintegrated 
in the interval J». Investigation of a series of such 
preparations would show whether A for RaF is a 
function of the ages of the parent atoms; but this 
is unnecessary, for we shall now prove by 
analysis of ordinary half-life experiments that the 
probability of disintegration per unit time is 
constant. 


II. GENERALIZATION OF WOLFKE'S THEOREM 


The general proof is complicated so we shall 
present the special case of a series having two 
unstable members, A and B, and a stable end- 
product, C. Of the A-atoms born in the interval 
da, an arbitrary number f(a)da will be included in 
our preparation. Collection is begun at time @, 
and ceases at time @2. Of the A-atoms born at 
time a, a certain fraction p(a, b, c)db dc will 
disintegrate in the interval b to b+db, while 
their daughter B-atoms will give rise to C-atoms 
in the interval c to c+dc. It is assumed that pis a 


smooth and continuous function of its variables, 
since otherwise our analysis would fail. 

By hypothesis, p depends only on the lives of 
the A- and B-atoms, so it is of the form p(b—a, 
c—b). The number of A-atoms in the preparation 
at time fis 


for this is simply the number whose death occurs 
after time 


Let us now form a preparation of B-atoms. The 
number of A-atoms decaying in the interval db 
is 


ff c—b)dc da-db, 


a=ay 


and of these we incorporate the fraction g(b) in 
our B-preparation. The number surviving at 
time / is 


ff f(a)g(b) p(b—a, c—b)de da-db, 


and so if accumulation proceeds from }; to be, we 
have 


ay bo 
= S(a)g(b)p(b—a, 


If Na(t) = for all such preparations, 
then 


f ff f(a)g(b) p(b —a, c—b)dc db da 
ay by 7 (4) 
= f da 
a, 0 


Since f(a) and g(6) are arbitrary, 


c—b)dc=e#' p(b—a, c—b)dc. (5) 
t 0 


This equation represents all we can learn about p 
by studying B-preparations. Differentiating with 
respect to ¢, using (5) and differentiating again, 
we get a differential equation whose solution is 
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p(b—a, t—b) = F(b—a)e*B, (6) 


where F(b—a) is an arbitrary function. Now 
suppose that N4(t)=Na4(O)e4', for all A- 
preparations. An argument similar to that above 
gives us F, and using the fact that 


ff (7) 
a 


we get 
p(b—a, c—b) pe. 


It follows in an elementary way that ws =), 
and wzg=Xz, which are the results desired. 

I wish to thank Dr. Morris Muskat for helpful 
discussion. 
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Section II reviews the general properties of the solutions 
of differential equations with periodic coefficients, the so- 
called Hill and Mathieu equations. In Section III asymp- 
totic formulae for the eigenvalues of the Mathieu equation 
in the oscillatory region (Eq. (9)) are obtained by applying 
the Schrédinger perturbation theory to the problem of the 


physical pendulum. Section IV applies the W. K. B. 
method to the Hill equations. Implicit equations for the 
eigenvalues are obtained (§9, Eq. (28), (30) §10). They are 
valid for E< Vmax. In the introduction (Section I) a group 
of physical problems are mentioned, to which the results 
may be applied. 


I. INTRODUCTION 


N a number of physical problems differential 

equations with periodic coefficients occur. 
These include, for instance, the problems of the 
vibrations of the elliptic membrane (Mathieu's 
original problem), the diffraction of light around 
an elliptic cylinder,' the variations of the orbit 
of the moon due to the sun’s attraction (Hill's 
original problem), quantum mechanics of elec- 
trons in a crystal,? quantum mechanics of the 
physical pendulum,’ the rotation of molecules in 
a crystal,‘ the torsional vibrations of the CHe 
or CH; group in the ethylene or ethane molecule.® 

Unfortunately our knowledge of the Mathieu 
functions and of the corresponding eigenvalue 
problem is still rather incomplete, especially on 
the numerical side, although valuable work in 
this direction has been done recently by Ince 
and Goldstein.* Some advances can be made, in 
my opinion, by starting from the physical 
problems and by the use of the approximation 
methods suggested by them. So we shall apply in 
Section III the Schrédinger perturbation theory 
to the problem of the physical pendulum from 
the ‘oscillator’ side and obtain certain asymp- 
totic series for the eigenvalues and the Mathieu 


1 See P. Epstein, Enc. der Math. Wiss. Vol. 5, p. 507. 

?F. Bloch, Zeits. f. Physik 52, 555 (1929). Ph. Morse, 
Phys. Rev. 35, 1310 (1930). 

3E. U. Condon, Phys. Rev. 31, 891 (1928). 

‘L. Pauling, Phys. Rev. 31, 430 (1930). 

5H. Nielsen, Phys. Rev. 40, 445 (1932). 

®Ince, Proc. Roy. Soc. Edinburgh 46, 20, 316 (1926); 
47, 294 (1927). J. London Math. Soc. 2, 46 (1927). Gold- 
stein, Trans. Camb. Phil. Soc. 23, 303 (1927). 


functions themselves. In Section IV we shall 
apply the well-known W. K. B. method to a 
slightly more general problem and_ obtain 
asymptotic expressions also for certain Mathieu 
functions of fractional order. 

In Section II we will first give a more quali- 
tative discussion of the properties of the Hill 
equation, of which the Mathieu equation is a 
special case. For all details and further develop- 
ments, especially on the mathematical side, I 
refer to Whittaker and Watson,’ Ince* and to the 
recent monograph of M. J. O. Strutt.® 


II. GENERAL THEOREMS ON THE HILL EQUATION 


$1. The Floquet theorem 
We will write the Hill equation in the form: 


d*y/dx?+(A—gF(x))y =0, (1) 


in which \ and g are constants; F(x) = F(x+ 7) 
and F(x) is so normalized that: 


f F(x)dx=0 | F(x)|31. (2) 
0 


The Mathieu equation is the special case corre- 
sponding to F(x) =cos 2x. The Floquet theorem!” 
states that the general solution of (1) is of the 
form: 


y = (x) + Cre**u2(x) (3a) 


’ Whittaker and Watson, Modern Analysis, Chapter 19, 

5 Ince, Differential Equations. 

°M. J. O. Strutt, Lamé sche, Mathieusche und verwandte 
Funktionen in Physik und Technik, Springer, Berlin, 1932. 

It is analogous to the Bloch theorem in the theory of 
metals. 
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or of the form: 
y =e (x) +Co{ xui(x) +u2(x)} J. (3d) 


In here u; and wu: are periodic in x with the period 
mw; mo, and moe in (3a) are the logarithms of the 
roots of a quadratic equation :" 


s*—2bs+1=0. (4) 


From this it follows that and aj, o2 are 
either both real or conjugate complex. Therefore 
o, and og can only be either both real or both 
pure imaginary. Eq. (3b) corresponds to the case 
that the quadratic Eq. (4) has two equal roots, 
of which zo is then the log. Because the double 
root of (4) can be only +1, o must be either 
zero or +1. 


§2. Stable and unstable solutions 


The o; and o2 in (3a) are certain complicated 
functions of \ and g, the determination of which 
is the main problem in the theory. 

The solutions where the o are real are called 
“‘unstable” solutions because they become infinite 
with +x. So if, e.g., y represents the displacement 
of a body from its orbit and x the time (as in 
Hill's problem), in such a solution the body 
would go farther and farther away from its 
orbit. The solutions where the o are pure 
imaginary are called ‘‘stable’’ solutions. They 
remain oscillatory in infinity. 

Because ¢=f(A, g) we can draw curves in a 
\, g diagram on which ¢ is a constant. One can 
prove that an infinite set of curves arises for 
each value of o. The o imaginary curves fill 
regions which we will call the stable regions. 
The o real curves fill the rest of the diagram; 
they give the unstable regions (see Fig. 1). The 
solutions of type (3), contain one stable and one 
unstable solution. They correspond to the curves 
which separate the stable and unstable regions. 
Each stable region is bounded by a curve cor- 
responding to o=0 and one corresponding to 
o= +1. In the case of the Mathieu equation the 
stable solutions corresponding to these double 
root values of o are the Mathieu functions of 
even (o=C) and odd (¢=+7) order. The cor- 


The principal values of the log have to be taken. The 
coefficient 6 is real. For a proof of these statements see 


Note 1. 
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Fic. 1. Plot of regions of stable and unstable solution 
for y’+(A+g cos 2x)y=0. Regions of stable solution are 
shaded. Points used to plot curves taken from Strutt's 
monograph, page 24. 


responding solutions of the Hill equation might 
be called the Hill functions. 


§3. The extreme cases 


In Fig. 1 we have four distinct regions: 
A<—g, no stable solutions; —g<A<0, almost 
all points correspond to unstable solutions; 
0<d<g, regions of unstable solutions and regions 
of stable solutions both present with regions of 
unstable solution covering greater area; \>g, 
greater part of area covered by regions of stable 
solutions. 

The region \>¢g is called the rotator side since 
in a quantum mechanical problem it is the case 
when the total energy which corresponds to \ 
is always greater than the potential energy, the 
maximum value of which is g."* The region 


In the quantum theory of metals this region corre- 
sponds to the case of nearly free electrons. It seems 
paradoxical that there are still unstable or ‘‘forbidden” 
regions, although they are narrow. It becomes clear, when 
one observes that there the de Broglie wave-length fulfills 
the Bragg total reflection condition." 

See further L. Brillouin, Die Quantenstatistik, etc., 
Berlin, Springer, 1931. 
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—g<X<0 is called the oscillator side since, in a 
quantum mechanical problem, it corresponds to 
the case when the system oscillates about a 
potential minimum." 

The transition region 0<A<g may be called 
the rotating oscillator region. In here the very 
broad stable bands of the rotator side go over 
into the extremely narrow stable regions of the 
oscillator side. 


§4. The proper solutions 

The proper solutions are those stable solutions 
having a preassigned period in x. Since u; and ue 
have the period z, the solutions having a period 
yx where vy is an integer, are those for which 
o=2ri/v where |r| is any integer =v/2 or zero. 

The problem of finding these proper solutions 
and the corresponding values of X (the eigen- 
values) as function of g usually arises in quantum 
mechanics where the potential V(@) has a period 
27/vin an angle Then @ and correspond 
to the same position so ¥(@)=y¥(0+27) if yp is 
to be single valued. Then after we transform to 
x so that V(x)= V(x+~7) we must have: 


In the following chapters we will only be con- 
cerned with these eigenvalue problems. We will 
distinguish between one, two, three, etc., minima 
problems according to the value of »(=1, 2, 3, 
etc.). The general character of the eigenvalues in 
these different cases can be easily seen from the 
Floquet theorem. Take first y= 1; then r can be 
only zero and we get, because ¢=0 and we only 
want the stable solution, the Mathieu or Hill 
functions of even order. Because the other 
solution in (3b) must be discarded as being 
unstable, it is clear that the corresponding 
eigenvalues are not degenerated.” They are 
given by the curves marked with a circle in Fig. 1. 
For v=2, r=0 or +1; we get, according to the 
definition in §2, both the Mathieu or Hill 
functions of even and odd order; the eigenvalues 
are again not degenerated. They are represented 
by the lines marked with either a circle or a 
cross. The case v=3 gives something new (see 


“In the theory of metals this region corresponds to the 
case of nearly bound electrons. 

Except for certain discrete values of g, where the 
eigenvalue curves may cross. See Note I. 
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Fic. 2. Eigenvalue curves for a v=3 problem: 
cos 2x)y=0; y(x)=y(x+3r). Curves marked 
= represent eigenvalues degenerate for all values of g. 
Dotted curves represent the o= +i curves and are not 
eigenvalue curves for this problem. 


Fig. 2). The value r=0 gives again the even 
order Mathieu functions with non-degenerate 
eigenvalues; r= +1 and r= —1 correspond to the 
two distinct solutions of Eq. (3a) with o,= 27/3, 
o2= — 21/3. The eigenvalues are therefore clearly 
doubly degenerate. For v=4 (see Fig. 3), we get 
the even and odd order Mathieu functions 
(r=0, »=+2) and a doubly degenerate eigen- 
value corresponding to r= +1. 

And so we can go on. In general it is easy to 
draw qualitatively the eigenvalues as function 
of g or to connect the eigenvalues at the oscillator 


Fic. 3. Eigenvalue curves for a y=4 problem. Same as in 


Fig. 2 in all other respects. 
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side with those at the rotator side. In the v- 
minima problem for g=0, the eigenvalues are 
given by: 


(Sa) 


when @ is again 2ir/v; they are then all doubly 
degenerate except the lowest. On the extreme 
oscillator side (g>+), one finds easily (see 


Section III): 
+g=4(n+})g}. (5d) 


The eigenvalues here are v-fold degenerate. To 
connect the eigenvalues for one extreme of g to 
those at the other extreme we need the oscillation 
theorem.* This theorem states that the eigen- 
values are ordered in such a manner that if A» 
corresponds to a proper solution having m 


nodes in the interval 0=x<yvz then \,.>A,, if 
m>n. It is easy to see that the number of nodes 
m corresponds to the integer m of (5a) and (5d) 
by: 


So the v lowest levels on the rotator side must 
join to the lowest v-fold degenerate level of the 
oscillator side; the second group of vy rotator 
levels joins the next oscillator level, etc. This 
enables us to draw the connecting lines, if one 
still remembers that all the levels are in the 
whole range of g doubly degenerated except those 
corresponding to the even and odd order Hill 
functions. The curves corresponding to these 
cannot be drawn without further analysis 
because, as is remarked in Note I, those with the 
same may cross a finite number of times." 


III]. PERTURBATION THEORY ATTACK 
§5. The quantum theory of the physical pendulum 
The Schrédinger equation for a pendulum of length r and having a bob of mass m oscillating in 


the gravitational field of the earth is:" 


d*y,/d02+ sin? =0, (6) 


where @ is the angle between the vertical and the string of the pendulum. We introduce the new 
variables given below and expand the sin? }4 term in a power series. 


x=4r0; k=4nxmg}/r; 


pt=1/r. (7) 


The equation now takes a form allowing the Schrédinger perturbation theory to be applied quite 


readily. 


dy /dx*+[a—2k? W=0. (8) 


The perturbation parameter is u. The zero order equation is then the well-known equation for the 
harmonic oscillator, with the eigenvalues a,” =4k(n+}) and the eigenfunctions: 


Wn? = (2k/ HH, (e(2k) 4), 


where //, is the mth order Hermite polynomial. Without any trouble one can carry out the per- 
turbation calculation to the second order in the eigenfunction and to the third order in a.'® The 


6 Ince, Ordinary Differential equation, p. 246. This question of connection was considered by Nielsen (reference 5) 
where one finds a figure in the case v =6. 

'7For an example where they cross once, see Strutt, reference 9, Fig. 3. In the case of the Mathieu functions 
(Fig. 2) they don’t cross at all. 

18 A more likely problem is a dipole oscillating in a plane in a homogeneous electric field. 

1 The matrix elements involved in the perturbation formulae can all be obtained with the help of the integrals: 


where m+n must be even. If m-+-n is odd the integral vanishes. 


(2+) 


((m+n)/2+j—p)! 
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results are given below: 
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= 
0 


a,°=4(n+ 

a,’ = —3(2n?+n+1), 
1 1 

a,? = —— — (6n*+9n?+9n+3), 
4k 4! 
3 


4k? 


The coefficients in the development of the perturbed eigenfunctions in terms of y,°, i.e., the coef- 


ficients dj? in: 


i=0 


> 


have also been calculated, at least for j=1 and 2. Because they take up much space and because 
they will not be used later, we will not reproduce them here. The writer will be glad, though, to 
communicate the results by letter to anybody who is interested in them. 


§6. Discussion 


For small values of » or large values of k the 
above formulae are very convenient for the 
numerical calculation of the eigenvalues. For the 
first three eigenvalues we have: 


ao/k = 2 —0.25p — 0.0312" p?—0.0278+ 
=6 +++, (9) 
= 10 — 


where p=u/k. To compare these results with 
the exact calculations of Goldstein we must first 
transform (8) to the standard form 


d*y /dx* + (4a+-16q cos 2x)y=0, (10) 

so that 
4ap?=4+ap/k, (11a) 
8q =1/p?. (11d) 


For the comparison see Table I. The first column 
gives g, which determines p. From (9) one then 


calculates a,/k, which, substituted in (11a),’ 


gives our value of a tabulated in column 2. The 
third column gives the exact values of a according 
to Goldstein and his symbols for the eigen- 
functions, the even order Mathieu functions. We 
see that if g is large enough the error is very 
small. 


TABLE I. 
Our Goldstein 
q —a —a Function 
2 4 6.005 6.0605 n=() 
3 4.898 9.615 9.015 Ceo 
2 4 2.3373 2.3326 n=1 
3 4.898 4.9593 4.9791 
4 5.656 7.8450 7.8385 Sez 
5 6.325 10.839 10.838 
5 6.325 5.183 5.110 n=2 
10 8.942 18.371 18.481 
20 12.65 49.087 49.382 Cee 
40 17.94 115.854 116,202 
100 28.28 329.972 330.137 


One must not expect that taking the higher 
order perturbations into account will always 
improve the result. It is fairly certain that (9) 
are not convergent but asymptotic developments 
of the eigenvalues in the extreme oscillator side. 
Even should the series for a, and y, converge, 
one cannot be sure that they would converge to 
the required solution because one does not take 
the periodicity condition into account. For this 
reason one must expect that the series for ¥,(x) 
represents the even order Mathieu function 
fairly well only for —x<xyu'<+~-. For values 
of x outside this region we may continue the 
solution periodically. Furthermore one sees from 
Fig. 1 that on the extreme oscillator side the 
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stable regions are extremely narrow so that one may expect that the series (9) will also represent 
the eigenvalues of all other stable solutions to within their inherent error. In Section IV we will 
calculate these extremely small splittings. 


IV. THe W.K.B. SoL_utTIon 
§7. Introduction 


We will treat in this chapter the v-minima problem with the well-known W.K.B. method. The 
physical problem which it may represent is, e.g., the quantum theory of a wheel (symbolizing a 
molecule or a radical like CHe) of 2v spokes of equal length. The wheel is rotating about a fixed axis 
and the ends of neighboring spokes carry opposite but equal charges. The wheel is placed in an 
inhomogeneous electric field (symbolizing the crystal or interatomic forces). The Schrédinger equa- 
tion would be: 

d*y /d0?+ /h?)(E—V(0))y=0 (12) 


where V(6@) = V(@+27/v) and we will suppose that V(@) has in each period 27/» one maximum and 
one minimum and is even around each of them. We measure @ from one of these minima and take § 
increasing counterclockwise. We choose the zero energy so that 


f =0. 


To apply the W.K.B. method we must distinguish between the classical and non-classical regions. 
The 2 roots of E— V(@) =0 are: 


6;= —a+2(i—1)x/», = +a+2(i—1)x/r, (13) 


where a is the first positive root. The region between @; and 6,’ will be called the 7th classical region 
(E—V>O0) and the region 0,’ the ith non-classical region (EK V <0). 
In each classical region the solution of (12) will be oscillatory and be given by: 


6 


where b; abbreviates 2(¢—1)/v and p is the angular momentum [ 2/(/£— V) }!. In each non-classical 
region the solution of (12) will consist of a decreasing and an increasing exponential: 


6 
f pao) + exp pao), (15) 


where 
bi =2(i—}) p=(21(V—E) }}. 


§8. Application of the connection formulae and of the periodicity condition 


The relations between the constants A;, B;, C; and 6; are determined by the connection formulae 
of Kramers.”® Applying these to extend the solution in the ith non-classical region to the ith classical 


2H. A. Kramers, Zeits, f. Physik 39, 828 (1926); H. A. Kramers and G., P. Ittmann, Zeits. f. Physik 58, 217 (1929); 
A. Zwaan, Diss. Utrecht (1929). These formulae are collected in a paper by L. A. Young and G. E. Uhlenbeck, Phys. Rev. 
36, 1154 (1930). For the application of the method to the two identical minima problem see D. M. Dennison and G. E. 
Uhlenbeck, Phys. Rev. 41, 313 (1932) and to the problem of two unequal minima, Ta-You Wu, Phys. Rev. in press. 
From the standpoint of these papers I treat the case of an infinite number of identical minima, where one requires periodic 
solutions. This problem has also been considered by M. J. O. Strutt (Math. Ann. 101, 559 (1929) using an analogous 
method with incorrect connection formulae. 
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region and requiring the result to be identical with (14) we obtain: 


A;,=C.8 cos (y/2+6;:+ 7/4), B;=(C,/28) sin (y/2+6;+ 7/4), (16) 
where: 

8 =exp (2x pao) =exp pao), (17) 

6,’ 

y=(2x/h) pd0=(22/h) pdé. (18) 


We extend the solution in the 7th non-classical region to the succeeding (i+ 1)th classical region in 
the same manner and obtain: 


A, = (Ci41/28) sin (y/2—6i41+ 7/4), cos (y/2—bi41+ 2/4). (19) 
Eqs. (16) and (19) immediately lead to a recurrence relation between the phases 6;: 
2p? cot (y/2+46:+ 7/4) = (1/28?) te(y/2—6i41+ 7/4) (20) 
and to a set of linear homogeneous equations for the amplitudes C;* 
cos (y/2+6:+ 2/4) = (Ci41/28) sin (21) 
Since ¥ must have the period 27, we must require: 
(22) 


Therefore (21) are vy homogeneous linear equations in the vy unknowns C;,. In order that the solution 
be not identical zero, the determinant of the coefficients must vanish giving: 


cos (y/2+6,+7/4) = (26) Il sin (y/2—64+ 7/4). (23) 
k=1 k=1 


At least one of the C; will now be arbitrary, which is just what one needs to normalize the solution. 
From (20) and (23) the simpler form follows: 


cos (7+28,) = cos (y¥—28i). (24) 


k=l 


Because of (22) Eqs. (20) and (24) become »+1 equations for the v 6; and y, which may be solved for 
yin terms of 8 and integers (the quantum numbers) only. Since, if V(@) is given, y and 8 depend on 
E only, this will give an implicit equation for the eigenvalues of E. 


§9. Solution in special cases 
This solution of (20) and (24) for small values of v is given below. Introducing the abbreviations: 


¢1=2 tan~'(1/26?), ¢2=2 [26?/1+46*+ (1+48*+ 166%) 


we obtain: 

v=1: y=(n+})r—(-1)"¢1; (a) 

y=2: y=(n+})r+¢1; 46,= —46.=2—(—1)"z, (b) 

v=3: 45, =46.=46; = r—(—1)"z, (a) 


" They arise from equating the first relations (16) and (19). Equating the second relations will give the same set of linear 
equations because of (20). 
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v=4: y=(n+4)r+¢1; 46, = —46,=46; = —46,= (d) 
y=(n+}4)z; 6,=6;=0, 62 =6,=6,+ 7/2, (c) 
v=6: 46,= —46,.= +46; = —46,=46, = —46,= 27—(-—1)"z, (b) 


y= 46,=46,= (—1)"z; —63;=6;= (0) 
y 52 = —63;=6;= F(—1)"}¢2. (6) 
The equations marked (a) correspond to the even order Hill functions, those marked (5) to the Hill 
functions of both orders, those marked (c) to the doubly degenerate solutions (see Figs. 2 and 3, 


Section IT). 
TABLE II. Splitting of first three eigenvalues for v=4, V(0) = —g cos 48. 


n q a’; a’, a’; |n q a’; a's a’; 

0 1 10.49 10.75 11.01} 2 3 78.64 80.93 83.05 
2 15.26 15.27 15.28 3 112.03 112.65 113.27 
3 17.79 17.80 17.81 8 146.48 146.52 146.56 

1 2 42.19 43.63 44.97 | 3 5 143.47 147.48 152.25 
3 51.26 $2.15 53.04 8 194.43 194.96 195.49 
5 72.17 72.17 72.17 12.5 251.72 251.76 251.80 


Table II shows the splitting of the first three eigenvalues in the case v=4, V(0@)=—g cos 46. 
We set: 


16g=8n*Ig/h? and 
to approach the standard form (10). The a’ is connected with the a of (10) by: 
a=a'—4q. 


The middle column gives the doubly degenerate eigenvalues (see Fig. 3). 


§10. The general solution 


We wish to correlate our W.K.B. solution with the Floquet theorem. To obtain the same notation 
as in Section II we introduce the variable x= so that V(x) =V(x+7) and =¥(x+pz). 
As in the proof of the Floquet theorem (see Note I) we will start with a fundamental set of W.K.B. 
solutions y;(x) and ye(x). For y; we choose the solution for which the phase 6; in the first classical 
region is zero and which is then extended over the whole region without regard for the periodicity 
condition. For ye we choose the analogous solution with 5,= 7/2. With the help of the connection 
formulae (20) and (21) we can write the characteristic equation analogous to (4) immediately: 


s?—2bs+c=0, (25) 

where now 
26 = sin (y+ 7/2), (26) 
166%)?+ 6484 csc? (y+ 7/2) (27) 


One sees that, in contrast to the case where one starts with two exact solutions, the coefficient ¢ 
is now not unity. However, since 8 is always >1, c is always nearly one (1<c<1.04). The difference 
is clearly due to our approximations. In the following we will take c=1. Calling the log of the roots 
of (25) again wo; and mez we have then o; = —o2=<a and from (26) follows: 


sin (y+72/2) = —46?/(484+1) cos ior =sech D cos (28) 
where 


Tv 
D=log 28? f pdx+log 2. (29) 


we 
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For the proper solutions we must take o = 2ri/v, |r| =v/2 or zero and then (28) is an implicit equation 
for E.* Comparing these results with the exact solution in the case of the Mathieu equation one finds 
that (28) gives too wide stable regions if E is near Vax. We can still improve the result (28) in the 
following manner. From (28) follows 


=(n+})r—(—1)" sin“ (sech D cos ior). 
Now if 8 is large, sech D is very small and we may put the sin! equal to its argument. So 


(n+4)x—(—1)" sech D cos ion. (30) 


If we now use (30) for all values of D, we get considerable better results for E near Vax, i-e., DElog 2. 
Furthermore, (30) exhibits for E near Vinax the so-called quarter quantization of Kramers and 
Ittmann (reference 16), which must be expected.” 

The author wishes to acknowledge his indebtedness to Professor Otto Laporte of the University 
of Michigan who suggested the problem and to Professor G. E. Uhlenbeck whose advice was of 
great value in carrying the work to a successful conclusion. 


Note I. 


Because most books do not contain the complete proof 
of the Floquet theorem we will give it here. 

Let y:(x), ya(x) be a set of fundamental solutions of (1). 
Then because (1) is invariant under the transformation 
x’=x+n, yi(x+7) and y2(x+7) will also be solutions of 
(1). We must have therefore:** 


=diyi(x) +di2y2(x), 


ya(x +7) =daiyi(x) +da2y2(x). (a) 


We wish now to find a solution y=c;y,+¢2y2 so that 
y(x+7) =sy(x). Using (a) we find that ¢;, must fulfill the 
equations: 

+¢2d,2 =0, 


+€2(s —d22) =0. 


The necessary and sufficient condition that (8) has at least 
one solution is: 


s—dy, diz 
dy, 

One proves easily that: 
= A(x) /A(0) =1, 


when A(x) =yiy2'—yoy:' is the Wronskian, which is a 


constant as a consequence of (1). We get therefore for s the 
quadratic equation (Eq. (4) of Section I1) 


(di: +dy2)s+1 (y) 


There are two main cases: 

a. (y) has two single roots. We can then find two of the 
required solutions y(x) and this leads to (3a). 

b. (y) has one double root, which must be +1. Supposing 
d,» and d2; not both zero, we can then find only one set of 
ci, ie., one required solution y(x). To find the second 
solution substitute y=xy+z. We find for z: 


2” +(\—gF(x))s = — 2dy/dx. (8) 


y has the period r(s = +1) or 2x(s = —1). The homogeneous 
part of (5) has y as the only periodic solution. We will be 
able to find a solution of (6) with period x or 27 if and only 
if y is orthogonal to the right-hand side, which is obvious. 
So we are led to (3). 

There remains the possibility that di: and d2, are both 
zero, Then the c are arbitrary and all solutions are 
periodic with period + or 27. We get again (3a) but with 
o, =02=0 or +i, so that both solutions are stable. One can 
prove that this can only occur for discrete values of g. 
They are the points where the eigenvalue curves for the 
Hill functions may cross (see §4). 


* Strutt (reference 16) obtains by his method the equation which would result if one omits the log 2 in the definition of 
D. That this cannot be correct follows immediately from the fact that then for E = V,,ax all solutions would be stable, 
contrary to the known results (see Fig. 1; E = Vmax corresponds to \=g). 

*8 Quarter quantization is to be expected for the limits of the stable regions if E is near Vinax. Then ¢ =0, +7, and sech 
D=0.8=}5 so that (30) goes over in: y 

** The y,, y2 may be chosen real and then the d;; will also be real. 


OCTOBER 15, 1933 


PHYSICAL REVIEW 


VOLUME 44 


The Energy Levels of the Rare-Gas Configurations* 


GeorGE H. SHort.ey, Palmer Physical Laboratory, Princeton University 
(Received August 4, 1933) 


The secular equations expressing the electrostatic and 
spin-orbit interactions of the configurations p*p and p'd are 
found to be very conveniently handled if one neglects the 
interaction between levels which have different levels of the 
ionic doublet as parents, an approximation significant for 
all but the lowest p*p’s of neon and argon. The 10 levels of 
pp and the 12 of pd are each expressed in terms of six 
parameters, one of which corresponds directly to the ionic 
doublet splitting. It is found that the higher members of 
the p*p series in neon and argon, the lowest configurations 
of this type in krypton and xenon, and all eight members of 


the neon p*d series agree with these formulas to within the 
magnitude of the above-mentioned neglected interaction 
(which becomes relatively very small for the higher series 
members) in all cases where two configurations do not ob- 
viously strongly perturb each other. Because of the gen- 
erally very weak interaction between configurations in these 
spectra, better agreement between theory and observed 
data is found than has been the case for any other atoms, 
This calculation gives the eigenstates for these configura- 
tions exact specification in the jj-coupling scheme. 


HE structure of the rare-gas spectra is 
recalled by Fig. 1, which is a plot of the 
energy levels of neon I, the most completely 
analyzed of these spectra, as given by Bacher 
and Goudsmit. All known levels of the rare gases 
are included in the systems n’p*nl where n’ = 2, 3, 
+++ for neon, argon, --+. This makes these 
spectra remarkably clean, almost as clean as 
one-electron doublet spectra. As seen in Fig. 1, 
the configurations of each series rapidly divide, 
with increasing m, into two groups of levels. The 
upper of these groups approaches the ?P;/2 limit 
of the ion, the lower the ?P3,2 limit. For p*s the 
upper and lower groups each consist of two levels, 
and for neon this configuration is completely 
known up to 2p°1lls. For p°p the upper group 
consists of four levels and the lower of six; this 
configuration is completely known in neon up to 
2p°7p. For 1>1, the upper group consists of four 
levels and the lower of eight; the 2p°nd of neon 
are completely analyzed up to n=10, while only 
two levels of the lowest p'f are found. 
The configuration p's has been considered in 
detail by Laporte and Inglis and by Condon and 
Shortley,! who used the formulas? derived by 


* Presented at the Washington Meeting of the American 
Physical Society, April 28, 1933. 

1 Laporte and Inglis, Phys. Rev. 35, 1337 (1930); Condon 
and Shortley, Phys. Rev. 35, 1342 (1930). 

2 That these formulas are applicable with change of sign 


Houston for the configuration sp. We shall here 
consider the configurations p°p and p'°d. 


The configuration n’p’np 


The configuration pp contains 1 level with 
J=3,3 with J=2, 4 with J=1, and 2 with J=0. 
Hence the diagonalization of the Hamiltonian 
for this configuration involves the solution of a 
linear, a quadratic, a cubic, and a quartic 
equation. These equations contain as parameters 
six radial integrals. Given the values of these 
parameters, the solution of these equations is 
rather complicated; to obtain the values of these 
parameters from the observed energies and then 
check the self-consistency of the equations, as we 
propose to do, is much more complicated. A 
simplification valid in many cases can, however, 
be made. Just as our consideration of one 
configuration at a time is an approximation in 
which we neglect the interaction between con- 
figurations, so we may use an approximation in 
which we also neglect the interaction between the 
groups of levels of the same configuration which 


of the spin-orbit parameter and hence of Houston's X is 
readily shown from previous results [Shortley, Phys. Rev. 
40, 185 (1932) ]. That X was found to be negative in these 
comparisons then shows that the electrostatic integral 
Go(sp) is positive, as is perhaps to be expected [cf. Bacher, 
Phys. Rev. 43, 264 (1933)]. 
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Fic. 1. Scale in thousands of cm™. Each line represents in general a group of observed levels 
not “resolved” on this scale. 


have different levels of the ionic doublet as 
parents. This latter approximation is expected 
(cf. Fig. 1) to be increasingly good as we go up 
the series of configurations. 

Now the quantum number which distinguishes 
levels belonging to one parent from those 
belonging to the other is the j value of the p® 
group, which is also the j value of the missing 


electron.* Hence in order to split the secular 
equation according to parentage, it is necessary 
to obtain the matrix of the Hamiltonian in the 
jj-coupling scheme. 

The matrix of spin-orbit interaction is diagonal 
in jj coupling, the diagonal element having the 


3 Cf. §2 of Shortley, Phys. Rev. 43, 451 (1933). 
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value 
(A| Fatt 
(1) 


where the summation is taken over all electrons 
in the configuration.‘ This element is independent 
of the J and My, values of the state A. The 
matrix for n’p'np is the same as that for n’pnp 
with reversed sign of ¢’={,’p». The values of the 
elements are 


n' 
a (34, 3g): 
b (3g, 14): ¢ 
(2) 
c (%, 34): 
d (14, Ws): 


In the parentheses are given the j values first of 
the p* core and then of the added p electron. The 
letters a, b, c, d are introduced as a convenient 
abbreviation for the four possible combinations 
of j values. States labeled by a and b belong to 
the lower doublet level, states labeled by c and d 
to the higher doublet level of the ion. 

The matrix of electrostatic interaction is 
diagonal in LS coupling; the diagonal elements 
are given on page 193 of reference 4. With the 
transformation of reference 3 this matrix is found 
in jj coupling to have the value 


Upper levels 
(2c}=—Fot 


(lc)=—Fot 
(1d) = — Fy+ 
(0d)=—Fot+ ¢ 


‘ Notation as in §5 of Shortley, Phys. Rev. 40, 185 (1932). 
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— Fot+ 3a —F, 


2a 2b 2c 
2a | | 


2b —2F2+4G, 4G2 — F:—4G, | 
la 1b lc ld 
la 2-S'F; —2-5'F, 
1d 2-5'F, 0 0 
1d 10°F, 0 0 0 
Oa Od 
Oa 2G, 


In the notation here used the J value is given as 
an arabic numeral, followed by a letter which 
specifies the electronic 7 values according to the 
scheme (2). The matrices are independent of the 
value of My. 

Adding the spin-orbit interaction (2) and 
splitting these matrices according to the parentj 
values as indicated by the broken lines, we obtain 
the following equations for the energy levels: 


+4G, 
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Lower levels 
FP: 
2a’ 
2b’ 
1a’ 


1b’ 
(Oa) = — Fo— — 


— Fo- 4S’ 


The six unprimed levels of this group belong in 
our approximation rigorously to the quantum 
numbers assigned to them. The eigenstates for 
the primed levels are linear combinations of those 
for the corresponding unprimed levels. We 
arbitrarity denote the higher of the two levels by 
a’ and the lower by 0’ since the level a would lie 
above the level } if the electrostatic interaction 
vanished. 

Now the value of one of the parameters, ¢’, 
should in all cases in which our approximation is 
good be obtainable from the splitting of the 
parent doublet, which is just 3s’. But this 
splitting is known from the spectra of the ions 
only to an accuracy of one wave number. Hence 
we prefer to determine this parameter, along 
with the rest, from the data for each configura- 
tion, and then to check it against the doublet 
splitting of the ion later. 

The six linear expressions in (4) should de- 
termine the six parameters. From the values of 
the parameters we could then predict absolute 
positions of the other four levels. However, in 
order to make the values of the parameters less 
sensitive to the small perturbations which must 
exist both between the two parts of one con- 
figuration and between neighboring configura- 
tions, we prefer to fit by least squares the eight 
quantities obtained by adding to the above six 
levels the means of 2a’ and 26’ and of 1a’ and 10’. 
This leaves the separations 2a'—2b' and 1a’—1b' 
fo be absolutely predicted. 

The results of such a calculation are shown in 
Fig. 2 for neon, Fig. 3 for argon, and Fig. 4 for 
krypton and xenon. In these figures the empirical 
value for 16’ has been taken as the zero for each 
configuration; a break in the wave-number scale 
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+4G2+[(345+3¢ Fe)? + (2F2—4G2)?]}! 


(4) 


cont, 
+[(34¢—14 


+4Go. 


indicates the separation between the two groups 
of levels. 

The first question one asks concerns the 
validity of the neglect of interaction between the 
two groups of levels. If the parameters we have 
obtained are approximately correct, the magni- 
tude of this interaction may be found by using 
the interaction elements of the matrices (3) 
according to the second-order perturbation the- 
ory. In this way the Jargest interaction was found 
to have approximately the value represented by 
the length of the short black bar drawn under- 
neath each configuration. This gives an idea of 
the agreement to be expected in the comparison. 
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Fic. 2. (The lengths of the black bars represent the 
maximum interaction between the two groups of levels in 
each configuration.) 
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Fic. 3. (The lengths of the black bars represent the 
maximum interaction between the two groups of levels in 
each configuration.) 


Xenon I 6p 


Fic. 4. (See text for explanation of xenon rearrangement.) 
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Having obtained the values of the parameters, 
one can readily find the coefficients in the 
expansion 


v(Ja’) =¥(Ja)(Ja| Ja’) +¥(Jb)(Jb| Ja’), 
¥(Jb’) =¥(Ja)(Ja| Jb’) +y(Jb)(Jb| Jb’). 


Since | (@]b’) |*=| (bla’)|2, 
and_ | one may all 
|(Ja|Ja’)\* the purity of the levels Ja’ and 
Jb’. It is, so far as this work is concerned, 
purely accidental that the electrostatic inter- 
action between 1c and 1d vanishes identically and 
hence that these levels are to our approximation 
100 percent pure. To the same approximation the 
separation between these two levels represents 
just the doublet splitting of the mp electron in the 
central field due to the core. 


NEON I (cr. Fic. 2) 


The 2p°3p is taken from the work of Inglis and 
Ginsburg,’ who solved the complete secular 
equation in LS coupling using the levels of 
J=0, 2, 3, and the mean of the four levels with 
J =1 to obtain the parameters, leaving otherwise 
the J=1 levels to be absolutely predicted. 

The 2)°5p is the first of the series to which our 
approximation is at all applicable. These levels 
agree to within the rather large maximum 
interaction between the groups. 

For 2p°6p, the interaction is much smaller and 
the agreement correspondingly better. 


TABLE I. Neon I 2p'np. 


3p 5p 6p 7p 

Fy — 1737 — 385.8 —314.7 —291.9 
F, 157.7 31.0 10.4 18.3 
Go 750.5 112.6 55.8 49.0 
G2 44.8 1.9 3.0 -17 
4 403 528.6 520.2 521.0 

40 8.1 5.4 13.2 
(2a| 2a’) 0.921 0.977 
(2b| 2a’) —0.390 —0.208 
(1a} 1a’) 0.683 0.697 
(1b| 1a’) 0.730 0.716 
% Purity 2a’b’ 84.8 95.5 
© Purity 1a’b’ 46.6 48.6 


5 Inglis and Ginsburg, Phys. Rev. 43, 194 (1933). With 
regard to this configuration see also St. Rozental, Zeits. f. 
Physik 83, 534 (1933). 
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2p*7p is found to agree not at all; the reason is 
that the lower 8 group exactly overlaps the 
upper 7p and is therefore expected to interact 
strongly with it. 

The values of the parameters and the coeffi- 
cients in (5) are given by Table I. The significance 
to be attached to these values must be judged by 
comparison with Fig. 2. The value of {’ is to be 
compared with the value 521 obtained from the 
parent doublet splitting of 782 cm“. 


ARGON I (cr. Fic. 3) 


The p*p configurations of argon remain rather 
well separated from each other up to 3p°8p, the 
last of the series which is completely known. 8p 0a 
falls about 200 cm™ below 7p 1d and the lower 
group of 11 falls close to the upper group of 8p 
as indicated in the figure. There is definite 
evidence of perturbation on 8p. The 5f group 
overlapping 6p does not have a pronounced 
effect. 


TABLE II. Argon I 3p'np. 


5p 6p 7p 8p 
Fo —805.8  —585.6 —5214 —515.7 
2 75.7 18.7 9.8 5.5 
0 235.6 95.5 53.6 28.3 
Gy 2.0 6.4 2.5 2.5 
“g 977.2 952.3 953.5 953.3 
44.1 7.1 3.6 5.0 
(2a| 2a’) 0.927 0.986 0.968 0.999 
(2b| 2a’) —0.378  —0.168 —90.254 —0.045 
(1a| 1a’) 0.702 0.690 0.690 0.723 
(1b| 1a’) 0.712 0.724 0.725 0.692 
% Purity 2a’b’ 85.9 97.2 93.7 99.8 
% Purity 1a’b’ 49.3 47.6 47.6 52.3 


The constants are given in Table II. The value 
of ¢’ is to be compared with the 954 obtained 
from the ionic doublet splitting of 1431 cm™. 


Krypton I AND XENON I (cr. FIG. 4) 


Even the lowest p*p’s of krypton and xenon 
are amenable to treatment with our approxi- 
mation. For krypton the 2p°5p and 6p are 
completely known. The 2p°5p agrees to within 
the error of the approximation. The 6p shows 
definite evidence of outside perturbation which is 
undoubtedly mainly due to the lower group of 7p 
which lies at only 6200 on the 6 scale, rather 
than to the 5f levels plotted. 
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In xenon only the 5p°6p is complete. This in 
the first plot agrees very poorly with the 
calculations. It is definitely perturbed by the 
lower 7p group which overlaps as indicated. 
Although one cannot make a rigorous assignment 
of these overlapping levels to configurations one 
suspects on comparing this figure with the others 
that the highest J=1 of 7p might belong more 
exactly to 6p than the 1d assigned to it. If one 
makes this rearrangement one obtains the plot on 
the right which shows a decidedly better 
agreement. 

The constants for these spectra have the values 
given in Table III. ¢ for krypton is to be 
compared with the 3581 obtained from the ion. 
The doublet splitting for the xenon ion is not 
known. 


TABLE III. 


Krypton I 4p'np Xenon I 5p*6p 


5p 6p Given Rearranged 
Fo — 2984 —2028.4  —4764 —4682 
F; 213 34.8 321 198 
Go 671 194.1 748 622 
2 20 16.3 —25 —2 
rad 3632 3567.6 7028 7094 
246 49.0 599 374 
(2a| 2a’) 0.955 0.999 0.954 
(2b| 2a’) —0.297 —0.022 —0.299 
(1a) 1a’) 0.736 0.750 0.775 
(1b] 1a’) 0.678 0.662 0.632 
% Purity 2a’b’ 91.2 99.9 91.0 
% Purity 1a’b’ 54.2 56.2 60.1 


The configuration n’p’nd 


The diagonal elements of spin-orbit interaction 
for p'd in jj coupling are 


pind 
b (33,34): 
(6) 
(1, Spt Fa 
d (4, 34): 


where we have introduced a notation a, b, c, d 
similar to that used for p'p. Here the parameter 
t», which is the same as the former ¢’, accom- 
plishes the splitting of the configuration into two 
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The electrostatic energies in LS coupling have the values 


groups, states characterized by a and bd lying in the upper group, those characterized by c¢ and ¢ 
in the lower. 


3F=—F,—2F, 1F= — F,—2F,+90G; 
7 
1D = — Fy+7Fs = — Fy—7Fs 
Transforming to jj coupling gives the electrostatic matrices 
4a 
— Fot+ 4a —2F, 
3a 3b 3c 
3a 23¢ Fe +24G; — 45-6! F24+12-61G; 84-53 F_—12-5'G, 
3b F:+36G; 14-301F)—6- 3056; 
3c 84-54 F,—12-5!G; —14-30! F,—6-30°G; 
2a 2b 2c 2d 
2a 46 Fs 25-219 Fs 86-149 96-214 Fy 
2b 216 Fy 16-64 Fs Fy 
2c —% 14°F, 14-6! Fs : 0 0 
2d 34-214 Fy 144 Fe 0 0 
la 15 1d 
la — 28¢F,+12G, 146 F2+4G, | 
1b 14¢ F,+4G, Fo+4gG, — 146-54 | 
1d 206G, | 
Ob 
Ob (8) 


Just as in the case of p*p we find if we neglect the interaction between levels belonging to 


different parents the following formulas to express the twelve levels in terms of six parameters: 


(3c)=—Fot Spt ba +30G; 
(2c)=—Fot fa 
(2d)=—Fot+ 
(1d)=—Fot +205G, 


(4a)=—Fo-MSpt fa—2 Fa 


(9) 
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= — Fo— Mf +[( 


Fo— Mb +[( 


(0b) = — Fa. 


We take the six levels given by linear ex- 
pressions and the means of 3a’, 3b’; 2a’, 2b’; and 
ja’, 1b’; fitting the six parameters to these nine 
quantities by least squares. This leaves the 3a’-3b’, 
2a'-2b’, la’-1b’ separations to be absolutely pre- 
dicted. 

From Fig. 1 we see that our approximation 
should be applicable to all the p*d's of neon I. 
Since these are known completely from 3d 
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2% F2—6G3+ 5464)? +96(3G3—1¢ 


54fa— 10 +8465 (9) 
cont. 


34 (4Gi 4145 Fe)?! 


through 10d we obtain the long series of con- 
figurations plotted in Fig. 5. These all agree 
within the accuracy of the calculation® except 6d 
and 7d, which are clearly perturbed. The most of 
this perturbation is due to the fact that the upper 
group of 6d lies only 26 cm™ below the lower of 


® The extreme coarseness of the wave-number scale for 
the right half of the plot is to be noted. 
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Fic. 5. Except for 6d and 7d, which perturb each other, the agreement is within the accuracy, indicated by the lengths 
of the black bars, of the neglect of interaction between levels belonging to different parents; in the cases 8d, 9d, 10d 
the maximum discrepancy is only 0.20 cm™, this agreement is within the accuracy of the calculation, 
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TABLE IV. Neon I 2p'nd. 


3d 4d Sd 6d 7d 8d Od 10d 

Fo —374.23 —306.59 —285.78  —273.97 —269.61 —266.52 —264.38  —263.55 
F; 15.507 6.595 3.620 2.008 1.433 0.900 0.615 0.467 
CG 2.98 1.70 1.04 0.64 0.47 0.29 0.20 0.16 
Gs 0.002 0.026 0.016 0.020  —0.004 0.004 0.001 0.002 
tp 530.11 521.62 520.40 519.10 519.49 520.24 520.20 520.35 
ta 1.178 0.002 0.141 0.476 0.064 0.043 0.000 0.044 
(3a| 3a’) 0.957 0.959 0.961 0.970 0.955 0.962 0.96 0.96 
(3b| 3a’) —0.285 —0.281 —0.276  —0.241 —0.297 —0.276 —0.29 ~0.27 
(2a} 2a’) 0.413 0.399 0.406 0.444 0.41 0.41 0.40 0.41 
(2b| 2a’) 0.911 0.917 0.914 0.898 0.91 0.92 0.91 0.91 
(1a) 1a’) 0.615 0.638 0.659 0.697 0.682 0.677 0.678 0.67 
(1b| 1a’) 0.789 0.770 0.751 0.718 0.733 0.736 0.735 0.74 
% Purity 3a’b’ 91.6 92.0 92.4 94.1 91.2 92.5 92 92 

© Purity 2a’b’ 17.1 15.9 16.5 19.7 17 16 16 17 

&% Purity 1a’b’ 37.8 40.7 43.4 48.6 46.5 45.8 46.0 45 


7d. That there is further perturbation on the 
upper group of 7d is shown by the abnormally 
large 2c-2d separation, which should be just the 
7d electron doublet splitting. This further per- 
turbation may be attributed to the fact that the 
lower part of 9d lies 100 cm! above the upper of 
7d and perhaps to the 2p°10s which lies about 20 
above 7d. 

The constants used in these calculations are 
given by Table IV. The value of ¢, is to be 
compared with the 521 obtained from the ionic 
doublet splitting. We should from the values here 
given predict a doublet splitting of 780.4+0.2 
cm in comparison with the observed 782 cm“. 

Thus the present theory accounts very satis- 
factorily for the observed structure of these rare- 
gas configurations. The perturbations which 
occur seem to be small except when two con- 
figurations of the same series overlap. Since there 
can be no spin-orbit perturbation corresponding 


to the interaction which splits the configuration 
into two groups, in estimating the mutual 
perturbation of two configurations we should 
compare the distance of the nearest levels not to 
the overall configuration size but to the much 
smaller spread of each group of levels. This 
requires that two configurations lie very close in 
order appreciably to affect each other. 

In those cases in which our approximations are 
valid we have obtained the actual eigenfunctions 
for each level in terms of the jj-coupling scheme.’ 
From these eigenfunctions we shall calculate the 
Zeeman effect, Stark effect, and line strengths for 
these configurations. 


7 The phases in this scheme are chosen in the following 
definite way (cf. reference 3). By use of Wigner’s formula 
(Gruppentheorie, p. 206) for the addition of 1 and 1] to obtain 
L, and s* are added to obtain and s* to obtain j*, 
and j* and j* to obtain J, each addition being made with the 
correlation to / and ] in the order specified. Here i refers to 
the p* core, e to the added p or d electron. 
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The actinium series arises from an isotope U*, or from 
two such isotopes U*** and U**, genetically connected. 
From work of Hahn and Meitner, the first possibility is 
almost certainly right. Assuming as a working hypothesis 
that U*** is the only long-lived uranium isotope, equations 
are developed for determining its decay constant A, from 
data concerning radioactive minerals and substances. 


Branching ratio: 0.03 0.04 
Ay in 2.28 x 1.79 
in yr! 0.1514 0.1509 


The best value for the half-life of U*** is (4.58 +0.09) 10° yr. 
If U*** exists and has a half-life long compared with that 
of U2, the values for \, would apply to U***. The ages, 
insensitive to the value of the branching ratio, are: 
Karlshus bréggerite, 0.81 10° yr.; Wilberforce uraninite, 
1.04 10° yr. Since quantitative study of chemical alter- 
ation may throw light on the process and prove helpful in 


These equations also give the decay constant \, of U** and 
mineral ages. The results depend on the value of the 
actinium “branching” ratio, which lies between 0.03 and 
0.04. Computations are carried out using both of these 
extreme values. Four minerals are discussed but only two, 
Karlshus bréggerite and Wilberforce uraninite, yield 
reliable results. The mean values from these are: 


Branching ratio: 00.3 00.4 
Half life of U**, yr. 3.10 108 3.96 x 108 
Half life of U**, yr. 4.58 x 10° 4.00 x 10° 


determining decay constants and mineral ages, equations 
are developed to show the effect of uniform leaching. Rates 
of removal or addition of Pb, U and Th are assumed 
different, but each of them is supposed constant in time. 
The resulting equations are applied to Katanga pitch- 
blende, to illustrate the method. 


I. THE PRESENT STATUS OF THE ACTINOURANIUM 
HYPOTHESIS 


T is generally accepted that the actinium 

series contains at least one member which is 
an isotope of uranium. Actinium has always been 
found in uranium minerals when search for it 
has been made. Further, the ratio of actinium to 
uranium is not far from constant and its devi- 
ation from constancy may be reasonably ex- 
plained as a result of alteration of these minerals. 
Hahn and Meitner? found that protoactinium 
slowly accumulates in pure uranium compounds. 
Their observations show that UY is derived from 
uranium; the problem which remains is to 
determine the details of their relationship. 

The earlier view that the actinium series 
originates by branching of the uranium series 
was discredited by Aston’s* mass-spectrograph 


' Preliminary notes: Phys. Rev. 42, 903 (1932); 43, 205 
(1933) and 43, 781 (1933). 

2 Hahn and Meitner, Berichte deutschen Chem. Ges. 54, 
69 (1923). 

3 Aston, Nature 120, 224 (1927); 123, 313 (1929); 129, 
649 (1932); Proc. Roy. Soc. A140, 535 (1933). 


studies of common and radiogenic leads in which 
he showed that some uranium leads contain 
amounts of Pb®’ too large to come from any 
source other than the actinium series. 

This means we must adopt the hypothesis of 
Piccard,‘ that the actinium series contains one 
or more uranium isotopes, genetically inde- 
pendent of U** and U**, Piccard supported a 
scheme of disintegration involving two so-called 
actinouranium isotopes, related through inter- 
mediate 8-rayers in the same way as the isotopes 
U** and U** (Fig. 1). Similar schemes involving 
two genetically connected isotopes, U** and 
U**, have been suggested by several workers but 
have never been supported by strong evidence. 
Hahn® has made extensive experiments to detect 
B-radiations which would be expected to occur 
in the hypothetical sequence of disintegrations 
connecting U*** and U*. He found no evidence 
of such radiations. It is fair to conclude from his 


‘ Piccard, Arch. Science Phys. et Nat. 44, 161 (1917). 
See also T. R. Wilkins, Nature 117, 719 (1926); and Phys. 
Rev. 29, 352 (1928). 

* Hahn, Zeits. anorg. Chem. 147, 16 (1925). 
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Fic. 1. Connection of the uranium and actinium series. 


results that (1) either U?*® does not exist; or (2) 
if it does exist, and if its immediate products are 
B-rayers, isotopic with UX, and UX2, these 
products must be so short-lived as to escape 
detection in spite of careful search. 

In the light of our knowledge of the half-lives 
of UX, UXse, and UZ, the first-mentioned pos- 
sibility appears much more probable than the 
second. If it is correct, then U** is the parent of 
the actinium series. Let us assume, as a working 
hypothesis, that this is true. Then from existing 
data, it is possible to obtain an approximate value 
of the decay constant of U** and thereby improve 
our knowledge of that of U?**; the ages of certain 
minerals can also be obtained.* This done, we 
may consider how the results would be modified 
in the improbable event that both U*** and U" 
exist. Decay constants of both isotopes cannot 
be rigorously determined at present because 
there are more unknowns than equations, but 
if one of the isotopes has a much longer half-life 
than the other, the previous analysis will furnish 
a good approximation to the decay constant of 
the more stable isotope. There are indications 


6 The calculations are based on the orthodox radioactive 
families (Fig. 1). 


that if both isotopes exist, U*** will have a much 
longer life than U**. Fajans’ has pointed out that 
for isotopic a-rayers the half-life increases with 
atomic weight, RaF being the sole exception. 

We shall often speak, for convenience, as 
though the existence of a single actinouranium 
isotope, U*, were an established fact. The 
reader is asked, therefore, to bear the above 
remarks in mind. 


II. Notation; Discussion OF THE ACTINIUM 
BRANCHING RATIO 


The subscripts 1, 2, 3 and 4 refer to the ura- 
nium isotopes U**8, U**4, the hypothetical U** 
and U**, respectively, while the subscripts Pa, 
Ac, Th, and Pb refer to the elements of which 
they are the symbols. Our notation is: 


N,, No, etc—numbers of atoms of the respective 
isotopes of U, per gram of U now present. 

N—the sum of N;, Nz and N,. 

N’ we, N’207, N’20s—numbers of atoms of the respective 
isotopes of Pb, per gram of U now present. 

Ne2oe, Noos—numbers of atoms of these isotopes 
produced by radioactive disintegration. 

Np», Nern, Ny-—numbers of atoms of the respective 
elements, and of common lead, per gram of U now 
present. 

Ai, Ao, etc., Apa, Aac, Arh—decay constants of uranium 
isotopes, etc. 

feos, f207, foos—fractions of Np, (total lead) constituted by 
Pb”, etc. 

426, 2207, @205—fractions of Ny, (common lead) constituted 
by etc. 

U/Pb, U/Th—ratio of mass of U to mass of Pb, etc. 

W,, Ws2, etc.—atomic weights on the chemical scale, 
0 =16. 


The actinium branching ratio, R, is defined by 
R=acNae/E, (1) 


where E represents one-half the number of a 
particles emitted per second by all the isotopes 
in one gram of uranium. Kovarik and Adams* 
found that E=12,385 a-particles/sec. per gram 
of U, with an uncertainty of 1 percent. The use 
of the ratio 
= (2) 
AN, 


7 Fajans, Le Radium 10, 171 (1913); Radioelements and 
Isotopes, p. 31, McGraw-Hill Book Co., N. Y., 1931. 
* Kovarik and Adams, Phys. Rev. 40, 718 (1932). 
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HALF-LIFE OF ACTINOURANIUM 677 
rather than R, simplifies many equations, so we shall adopt it. The available values show a con- 
siderable spread and we can only state® that with high probability B lies between 0.03 and 0.04. 
Calculations will be made using both these limiting values. 


III. EQUATIONS FOR DETERMINING HALF-LIVEs OF URANIUM ISOTOPES, AND MINERAL AGES 


On the working hypothesis adopted here, that the actinium series contains only one long-lived 
uranium isotope, we have just enough data to calculate its decay constant and the ages of certain 
minerals. The relations involved are as follows. 

Since the intermediate members of the uranium and thorium series are so short-lived as to be 
present in negligible amounts, we may write 


Noose = Nile*—1); (3) Neos = (4) 


If the actinium series originates in U*,!° 


New = N,(er# = 1). (5) 
Also, 
(6) Br Nj. (7) 
The observed quantities of the lead isotopes consist of original and radiogenic lead, as indicated by 


the equation, 
N' 206 =d206N i+ Noo, (8) 


and two similar relations, (9) and (10), giving N’s07 and N's. From the definition of Avogadro's 
number, 


6.064 (11) 
and from the experiments of Kovarik and Adams, | 
4 = 2(3.9083 X per year. (12) 


These equations contain ten unknowns: Ny, Ne, N4, A1, As, Neos, Noor, Neos, Nx, and ¢." It is easy 
to eliminate all unknowns except Ay, 44 and ¢. With times in years, and decay constants in yr.~', the 
result is, 


3.9083 10" 
(13) 
6.064 10(1+B, 2) Wi, 
206A1 4206 
(eit—1)—B (et—1) = (J fun ~); (14) 
N A207 
fe 206 W 4206 
feos Won U 


® Western and Ruark, Phys. Rev. 42, 903 (1932); 43, 205 
(1933). 

“If the series originates in U*** and if U*® is short- 
lived compared to U2", we would write 


Noo: = —1); 


and Eqs. (9) and (10) would contain additional terms. 
" Walling found the half-life of U2*4 to be 342,500 years 
[Zeits. f. physik. Chemie B10, 467 (1930)]. Collie [Proc. 


Roy. Soc. A131, 541 (1931) ], concluded that this half-life is 
not less than 10° years, but it seems unlikely that it is 
much more than 10° years. The mass of U*"* in uranium can 
be neglected in comparison with the masses of U*** and 
U*®, although its activity, of course, cannot be neglected. 
For our purpose, it is unnecessary to know A, exactly. The 
value of Xp, was redetermined by Kirsch, who finds the 
half-life to be 1.8 10" yr. [Phys. Zeits. 31, 1718 (1930); also 
Kirsch and Lane, Proc. A. A. A. S. 66, 357 (1931)]. 
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In (14) and (15) we use the very close approximation, V/N,=1+B),/\s, obtained by neglecting Ny 
as compared with N, and N,. It is convenient to write T=10~%, L; = 10°, /2.3026, etc., so that T 
is the mineral age in billions of years, and to rewrite the equations in a form adapted to solution by 
approximations. Inserting the values of coefficients of a general nature and replacing exponentials 


by powers of 10, we have” 


L,=([1/(1+B/2) }0.06663(1+0.9874B(L,/L,)); 
10%7—1=(Npp/N) (feos —0.5600f20s) (1 +.B(Li/Ls)) +0.5746(Th/U) (10-9727 —1)(14+B(Li/L,)); (17) 


1.374BL,T 


(16) 


10471 (10417 — 1) —(Npp/N) (feos —1.374f207)(1+B(Li/Ls)) 


IV. APPLICATION OF EQUATIONS TO MINERAL 
DaTA 


Table I gives the best values now obtainable 
for the isotopic and atomic weights of lead, 
uranium, and thorium.” The uncertainty of 
these figures is 0.03 atomic weight units, with 
the exception of that for common lead, where the 
uncertainty is 0.007 unit. Table II shows relevant 
portions of Aston’s'* mass analyses of leads, and 
Table III gives analyses of three of the minerals 
from which these leads were derived. The 
analysis of Katanga pitchblende is not available. 


TABLE I. Isotopic and atomic weights. 


Atomic weight 


Element or isotope Chemical scale 


Lead, 206 205.98 
207 206.98 

208 207.98 
Ordinary lead, chemically 207.22 
Thorium, calculated 232.03 
Uranium, 234 234.04 
235 235.04 

238 238.04 

239? 239.05 

Uranium, natural mixture 238.04 


% To effect a solution, one neglects at first the terms 
containing B. By using a rough value of T in the thorium 
term, an approximate value for 10“7—1 is obtained from 
(17); this and the corresponding value of L,T are substi- 
tuted in (18), which gives Ly7. L,/L, is computed and is 
used in (16) to obtain a fair value of L;. T is calculated from 
L,T and the process is repeated. The convergence is good 
for uranium minerals and little would be saved by using 
approximate equations. A graph of x/(10*—1) facilitates 
the work. 

13 Western and Ruark, J. Chem. Phys. 1, 717 (1933). 

4 Proc. Roy. Soc. A140, 535 (1933). 

18 Calculated on the hypothesis that U*® is the only 
actinouranium isotope by using an approximate value of 
the half-life of this isotope. 


(18) 


Of the minerals described, the thorite is not 
satisfactory for determining uranium decay 
constants because of the susceptibility of thorium 


TABLE II. Isotopic constitution of ordinary and 
radiogenic leads. 


Katanga Wilber- 
Ordinary Norwegian __ pitch- force Norwegian 

lead _bréggerite blende uraninite thorite 

Chemical atomic weight: 
207.22 206.14" 206.048 206.195 207.904 

Percentage 206: 

27.75 86.8 93.3 85.9 4.6 
Percentage 207: 

20.20 9.3 6.7 8.3 1.3 
Percentage 208: 

49.55 3.9 (0.02) 5.8 94.1 
Other isotopes: 

25 0 0 0 0 


* Raade, Norway. Calculated from isotope-percentages, 
and isotopic weights of Table I. 

>» Hénigschmid and Birckenbach, Ber. 56, 1837 (1923); 
von Grosse and Kurbatow, quoted by Aston, Proc. Roy. 
Soc. A140, 535 (1933). 

° Baxter and Bliss, J. Am. Chem. Soc. 52, 4851 (1930). 

4 Fajans, Sitz. Heidelb. Akad. Wiss. 3 (1918); Hénig- 
schmid, Phys. Zeits. 19, 437 (1918); Zeits. f. Elektrochemie 
25, 91 (1919). 


TABLE III. Mineral analyses. 


Norwegian Wilberforce Norwegian 

bréggerite* —_uraninite® thorite® 
Percentage U 61.16 53.52 0.45 
Percentage Th 4.38 10.37 30.10 
Percentage Pb 8.02 9.26 0.35 
Th/U 0.0716 0.194 66.1 
Pb/U 0.131 0.173 0.78 


* Fenner and Piggot, Nature 123, 793 (1929). 

> Wells, reported by Baxter and Bliss, reference c to 
Table II. 

° Fajans, Zeits. f. Elektrochemie 24, 163 (1918); Phys. 
Zeits. 19, 438 (1918), 
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minerals to leaching and because of the small 
amounts of Pb** and Pb®? in the mineral. Decay 
constants computed from this mineral are listed 
in Table IV as a matter of record but are devoid 


TABLE IV. Calculated decay constants and mineral ages. 


Thorite Bréggerite Uraninite 
B=0.03 B=0.03 B=0.04 B=0.03 B=0.04 
Age (yrs. X 107°) 

0.34 0.805 0.809 1.036 1.041 
Ai (per yr. X 10°) 
0.1512 0.1514 0.1508 0.1515 0.1510 
life (yrs. X 107°) 
4.58 4.578 4.574 4.590 
12+1 2.62 2.07 1.94 1.51 
0.056 0.357 0.459 


of significance. The data on the uraninite and 
the bréggerite appear fairly trustworthy. It is 
well known that alteration of a uranium mineral 
increases the lead-ratio,© Pb/(U+kTh). The 
uraninite, described in detail by Kirsch and 
Lane,” has a lead-ratio slightly higher than that 
of some samples from the same vicinity but the 
difference is not great enough to affect our cal- 
culations materially. Fenner and Piggot!® de- 
scribed the bréggerite as being without evidence 
of alteration. Comparison of its analysis with 
Gleditsch’s'® analyses of five samples from the 
same mine shows that, while its uranium content 
isin agreement with the lowest value of uranium 
from these, both the thorium content and the 
lead content are only about ninety percent of 
the lowest values found for these quantities in 
the Gleditsch samples. The low thorium content 
may indicate alteration. If so, calculations based 
on the Gleditsch samples show that the maximum 
effect of. leaching on the values calculated for \, 
and ¢ from the Fenner-Piggot sample would be to 
raise the one, and lower the other about five 
percent. A; is insensitive to such alteration. We 
assume, accordingly, that the mineral is unal- 


“Here rn. By using Kirsch’s 
value of Arh, k=0.26. 

Kirsch and Lane, Proc. A. A. A. S. 66, 357 (1931). 

Fenner and Piggot, Nature 123, 793 (1929). 

 Gleditsch, Ark. Mat. Og Naturv. Kristiania 36 (1), 
184 (1919); Norske Videns. Akad. Oslo I, Mat. Nat. KI. 
No. 3 (1925). 
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tered and base our computations on the Fenner- 
Piggot analysis. 

Values calculated from the three minerals 
discussed above are given in Table IV. The 
agreement between the results from the brég- 
gerite and uraninite is satisfactory. It is unfor- 
tunate that the required data are not available 
for more good uranium minerals, to indicate 
whether this agreement is genuine.” 

Kirsch and Lane” have shown that the age of 
the Wilberforce uraninite is (1.03+0.10)10° 
years. 


V. Tue Errect oF LEACHING ON RADIOACTIVE 
AGE DETERMINATIONS 


In recent years techniques have been devel- 
oped for judging whether or not a mineral 
specimen is essentially unaltered.” The principal 
aims of such work have been the selection of 
unaltered specimens, the establishment of reli- 
able lead-ratios and the attempt to understand 
why leaching increases the lead-uranium ratio 
in uranium minerals, while it decreases the lead- 
thorium ratio in thorium minerals.” 

It is desirable to begin a quantitative study of 
leaching which, as our knowledge increases, may 
throw additional light on the process and prove 
an aid in determining the ages of certain minerals. 


2° The kolm of Griillhégen, Sweden has attracted much 
attention. The atomic weight of the lead is 206,013. 
Aston (Proc. Roy. Soc. A140, 535 (1933)) obtained a single 
mass-spectrum which showed Pb* and Pb*”’ but he states 
that this spectrum was unsuitable for photometry. Previ- 
ously several writers calculated the composition of the lead 
on the assumption that it contains only Pb and Pb*”, 
The results depend markedly on the values used for the 
isotopic weights. By using the weights in Table I, such a 
calculation gives 96.7 percent and 3.3 percent, respectively, 
for these two isotopes. However, these values are so 
sensitive to the isotopic weights that with the limits of error 
placed on the latter, the value computed for \4 can vary 
from about 0.015X10~* yr.-! to 3.6 10~* yr... Without 
a more complete mass analysis, the kolm cannot give 
information on the decay constant of actinouranium. The 
mineral age is 0.42X10* yr., with an uncertainty of 4 
percent. 

% Kirsch, Geologie und Radioaktivitat, pp. 132-180. Work 
of Ellsworth, Todd, Davis, and Hecht and Kérner is 
described by Holmes on pp. 324-8 and 365-6 of The Age 
of the Earth, Bulletin 80 of the National Research Council, 
Washington, D. C. (1931). 

2 Holmes, Phil. Mag. 1, 1055 (1926). 
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It is possible to improve our understanding of 
the effects of alteration by considering three 
hypothetical cases: 

(1) Initial alteration. In this case practically 
no radiogenic lead is present when leaching 
occurs so the effect of alteration is the same as 
that of a decrease in the initial concentrations of 
U and Th. 

(2) Recent alteration. Here the composition of 
the lead is almost unaffected. With a knowledge 
of the value of \4, isotopic analysis would make 
possible a correct age determination by equations 
(3), (5), (7), (8) and (9), even though the ratios 
Pb/U and Th/U were untrustworthy. This case 
may be realized in a rough way by some ura- 
ninites; for Ellsworth,” discussing uraninites of 
Villeneuve, says, ‘“The lead ratio fortunately 
increases only slightly so long as the mineral 
remains black due to an appreciable UO, 
content, but once the UO, is entirely oxidized, 
uranium is lost much faster than lead, and the 
lead ratio rises to a large and entirely misleading 
value.” 

(3) Uniform alteration. It may be that for 
some minerals, the net effect of alteration is 
nearly the same as if it had taken place at a 
uniform rate. The simplicity of this assumption 
recommends it as a first approximation. If, in a 
radioactive series, the subscripts p and d refer 
to the long-lived parent element and the end- 
product, respectively, and if the intermediate 
elements can be neglected, uniform leaching is 
defined by 


dN,/dt= (AptRp) Np; 
dNa/dt=r,N,—kaNa, (19) 


where k, and ka may be called alteration con- 
stants. If we put k=k,—ka, 


Na=[ho/(Ap +k) (20) 


Application of (20) to each of the radioactive 
series gives three equations containing t, ky —Rkpp 
and ky,—kp» as unknowns. Thus ¢ can be ob- 
tained directly but in order to use these equations 
in finding alteration constants, we need the 
analysis of the mineral, the isotopic composition 
of the lead it contains and the analysis of an 


2H. V. Ellsworth, Am. Mineral. 15, 455 (1930). 


unaltered specimen from the same deposit, jf 
there be such. 

As an illustration, the method may be applied 
to the Katanga pitchblende of Table II, by 
making some assumptions of uncertain validity, 
The lead was separated from a composite sample 
of altered material by Hénigschmid and Bircken- 
bach.* They made no mineral analysis, feeling 
that the lead-ratio would be without significance, 
However, data on pp. 367-368 of The Age of the 
Earth, show that for three altered specimens 
designated as J, J and K, the mean ratio of lead 
to uranium is 0.096, while for several unaltered 
specimens it is close to 0.080. We write (Pb/U)’ 
for the latter ratio to distinguish it from the 
value for the leached specimen. Thorium is prac- 
tically absent. Assuming these data are appro- 
priate to the material of Hénigschmid and 
Birckenbach, we can write equations which 
determine Ay, \y, Ru —kp» and ¢. Eq. (13) applies 
directly. For an unleached thorium-free mineral, 


205.98 N2os+206.98 Noor 
238.04 NV 


which may be written 


238.04 AL 
205.98 
206.98 
205.98 


Further, we have two equations obtained by 
applying (20) to the uranium and actinium 
series. The results are: \,;=0.151410~° per 
year, \y=2.8X10-* per year, ky—kp,=0.64 
x10-* per year and ¢=0.55X10° years. Only 
the value obtained for 4 can yield any test of 
the assumption of uniform alteration. In this 
case the character of the data is such that, even 
if \, were accurately known, no definite con- 
clusion could be drawn. It is believed that when 
\, is more accurately known and further mineral 
data are available, similar investigations may 
throw much light on alteration processes. 


* Hénigschmid and Birckenbach, Ber. d. D. chem. 
Gesellschaft, p. 1837 (1923). 
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VI. CoNcLUSIONS 


In view of our present knowledge, a value of 
three or four hundred million years for the half- 
life of the principal actinouranium isotope cannot 
be far from correct. Further study of the ac- 
tintum branching ratio and complete data on 
additional minerals are needed to yield a more 
accurate value. 

The value calculated for the half-life of U** is 
relatively insensitive to that of actinouranium. 
Consequently, the uncertainty in the former is 
no larger than the uncertainty in our knowledge 
of the rate at which uranium gives off a@ par- 
ticles. Judging from the small mean deviation 
in the counting experiments of Kovarik and 
Adams, a liberal estimate of all uncertainties 
allows us to state the half-life of U®* as (4.58, 
+0.09) X 10° years, assuming that U** does not 
exist. Kovarik and Adams* gave 4.52, 10° 
years. This was computed on the branching 
hypothesis and the value used for the atomic 
weight of uranium was 238.17, which is too 


high. Our value should be considered as replacing 
that of Kovarik and Adams, since we have 
employed their data on the rate of emission of a 
particles by uranium, and have taken account 
of U**, The value of the International Radium- 
Standards Commission,® 4.4 10° years, is not 
comparable with ours. The Commission states 
that no account was taken of the actinium series, 
because the half-life of actinouranium was not 
sufficiently well known. 

The U** content of uranium lies between 0.15 
and 0.35 percent, provided U*** does not exist. 
The wide variation given is largely due to uncer- 
tainty in our knowledge of the branching ratio. 
On the other hand, if U*** exists and is long-lived 
in comparison to U**, the same figures would 
apply to the former. A content of U*** ten times 
higher than the larger value given here would not 
explain the present atomic weight, 238.14, 
assigned to uranium. 


2° Rev. Mod. Phys. 3, 427 (1930). 
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Electro-Optical Kerr Effect in Gases 


C. W. Bruce, University of Virginia 
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The Kerr constant B of CO, has been measured as a 
function of density from 0.06 to 0.32 g/cm® and as a 
function of temperature from 7.2 to 34.6°C. The results are 
in agreement with the Langevin-Born theory. The absolute 
value of B for CO, has been measured and is 0.34 x 10~"° per 
atmosphere for wave-length 4550A and temperature 34.6°C. 
This is also the theoretical value of B for CO,. The probable 


error of the measured value is 2 percent. The absolute 
value of B for N2 measured by comparison with CO, js 
0.057 X10~" for the above wave-length and temperature, 
with a probable error of 3 percent. The corresponding 
theoretical value for N2 is 0.050 10~". The Kerr constant 
of H, is shown to be less than 0.01 x 10~" which is not in 
disagreement with the theoretical value of 0.005 x 10-™, 


I. INTRODUCTION 


N 1875 Kerr' discovered that if an optically 

isotropic substance be placed in an electric 
field it acquires the optical properties of a 
uniaxial crystal with the optic axis in the 
direction of the electric field. The law governing 
this phenomenon, proposed by Kerr, states that 
if m, and mz are the refractive indices for the 
components of the light vector vibrating parallel 
and perpendicular to the lines of force in the 
substance, respectively, their phase difference D 
in radians after passing through the field F, in 
absolute electrostatic units, is 


(1) 


where / is the length of the light path through the 
field, \ is the wave-length of the light in a 
vacuum and B is the so-called Kerr constant of 
the substance. 

The theory of the Kerr effect as worked out by 
Langevin,? Born’ and others connects the Kerr 
constant with the index of refraction, dielectric 
constant and the light scattering coefficient 
(depolarization constant).*: 

Langevin assumes that the molecules are 
electrically and optically anisotropic and, as a 
result of the orienting action of the electric field 


! Kerr, Phil. Mag. (4) 50, 337 (1875). 

* Langevin, Le Radium 7, 249 (1910). 

3 Born, Ann. d. Physik 55, 177 (1918). 

4See Beams, Rev. Mod. Phys. 4, 133 (1932). 

5 See Szivessy, Handbuch d. Physik Vol. XXI, Julius 
Springer (1929). 


upon the induced doublets in the molecules, the 
substance as a whole becomes doubly refracting. 
Born extended the theory to include the perma- 
nent moment of the molecule. Raman and 
Krishnan® have further modified the Langevin- 
Born theory so that it becomes applicable to 
liquids in which the contribution of the molecules 
to the field becomes large. The quantum me- 
chanics of the Kerr effect has yielded little more 
than these classical theories.’ 

From the Langevin-Born theory, for a non- 


polar gas 

(2) 
6=(1/45kT)[ (a; —a2) (b; — be) 

+ J, (2a) 


where B,, is the Kerr constant per molecule; J, 
the number of molecules per unit volume; , the 
ordinary index of refraction; ¢, electrostatic value 
of dielectric constant ;k, Boltzmann's constant; 7, 
absolute temperature; ;b2b;, moments induced in 
the molecule along its three principal axes of 
optical anisotropy by unit electrical force in the 
light wave acting along the three axes re- 
spectively; and a,dea3, moments induced along 
the same axes by unit electrostatic force. From 
(2) and (2a) 


B=B,,N « N(n*+2)*(€+2)?/n (3) 


6 Raman and Krishnan, Proc. Roy, Soc. A117, 1, 589 
(1927). 

7 See Van Vleck, Electric and Magnetic Susce ptibilities, p. 
366, Oxford University Press (1932). 
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if \ and 7 are held constant while the density is 
varied. It is also seen from (2) and (2a) that if V 
and \ are held constant and 7 is varied BT 
=constant, since, for a nonpolar substance, « and 
n are practically independent of the temperature. 

From a result of the theory of light scattering 
by molecules* and an assumption due to Gans’ it 
may be shown that for a nonpolar gas at low 


densities, 
(n?—1)(e—1) 3A 
B= (4) 


where A is the light-scattering coefficient (de- 
polarization constant) of the molecule. Formula 
(4) is sufficiently accurate for this experiment. 

The theory for the transmission of plane 
polarized light through the Kerr cell, which is the 
same as through a uniaxial crystal, shows that 
the intensity of the Kerr beam is 


IT=I, sin? (D/2), (5) 


where Jo is the intensity of light through the 
Kerr cell and D is the phase difference in radians 
between the ordinary and the extraordinary ray. 
Using the value of D from (5) in (1) gives 


B=(1/nlE?) aresin (I/Io)!, (6) 
which for most work in gases may be written 
B=(1/nlE*)(I/Io)!. (7) 
From (7) 
(8) 


if Ip9Z, and / are constant. Eq. (7) is the one by 
which B is determined from experiment. The 
error in using (7) instead of (6) in reducing the 
measures in this experiment is never greater than 
0.1 percent. 

Three tests may be applied to test the validity 
of the Langevin-Born theory. These are: (1) a 
comparison of the experimental value of the Kerr 
constant as a function of density with the 
theoretical value as given by Eq. (3), (2) the ex- 
amination of the constancy of BT when the 
density and the wave-length of the light are held 
constant while the temperature is changed and 
(3) a comparison of the absolute values of B 
given by Eq. (4) with the values given by ex- 


®See Marx, Handbuch der Radiologie, p. 781. 
* Gans, Ann. d. Physik 65, 97 (1921). 


periment using Eq. (7). Needless to say experi- 
mental values of m and ¢ are used in Eqs. (3) and 
(4). These three tests may be taken as the ob- 
ject of this paper. 


Il. EXPERIMENTAL ARRANGEMENT 


A modification of the method of Stevenson and 
Beams" is used. The novelty of their method 
consists in placing the polarizer and analyzer 
inside the Kerr cell containing the gas under high 
pressure. This method eliminates the effect of 
strains in the windows of the Kerr cell and allows 
higher field strength to be used. On account of 
difficulty with the nicol prisms under high 
pressure they were replaced by a double image 
polarizer and analyzer. Fig. 1 shows the optical 
system of the final arrangement. C, is a photronic 
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Fic. 1. Optical system and circuit of d.c. amplifier. 


cell used to correct for small variations in the 
light source F. The light from F is collimated by 
L,. The stop 5S; limits the size of the beam of light 
entering the polarizer P which is a calcite 
crystal, set so that the plane of polarization of the 
light makes an angle of 45° with the field between 
the plates K of the Kerr cell. S,; cuts out the 
extraordinary ray. The analyzer is a Cotton" 
prism made of crystal quartz. It is set parallel to 
the polarizer so that with no field on the Kerr cell 
there is a single beam of light coming from the 
analyzer, called hereafter the light beam. With 
the field applied light through the Kerr cell 


© Stevenson and Beams, Phys. Rev. 38, 133 (1931). 
" Cotton, Comptes Rendus 193, 268 (1931). 
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becomes elliptically polarized so that two beams 
of light leave the analyzer, the ordinary ray 
comparable in intensity with the light beam and 
the extraordinary ray of small intensity, here- 
after called the Kerr beam. The lens Lz brings the 
light to a focus at S;, which may be set to cut out 
the ordinary ray but to allow the Kerr beam to 
pass through the filter B to the evacuated 
potassium hydride photoelectric cell C:, used in 
conjunction with a d.c. amplifier to measure the 
light intensity. A rotating sector / driven at a 
high speed by a motor may be placed in the path 
of the light beam to cut down its intensity to 
approximately that of the Kerr beam so that by a 
slight shift of Ss; it could also be measured to 
obtain absolute values of the Kerr constant. The 
plates K are heid in position by Bakelite rings 
and held apart by Bakelite separators 4 mm 
thick. The plates are 90 cm long and 2.5 cm 
wide. Electrical connections with the plates are 
made through ordinary spark plugs, the bases of 
which are soldered into holes bored in the tube. 

Fig. 1 also gives the electrical circuit of the 
d.c. amplifier. By a separate experiment the 
deflection of the galvanometer was found to be 
proportional to the light falling on C, in the range 
used. The resistance R is 10" or 10"! ohms 
depending on the sensitivity required. This entire 
circuit was placed inside a metal box for shielding. 
The inside of the box was kept dry by P:O;. This 
circuit is a modification of the one given by 
DuBridge.” 

Fig. 2 shows the construction of the analyzer 


10000 


Fic. 2. Construction of end of Kerr cell. 


”® DuBridge, Phys. Rev. 37, 392 (1931). 
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end of the Kerr cell. Copper washers at C and 
lead foil at L form gas-tight seals. The window W 
is glass 3 cm thick and 3 cm in diameter. By 
the worm gear G the analyzer may be set 
accurately parallel to the polarizer. The barrel of 
the Kerr cell is a cold-drawn steel tube 3.8 cm 
inside in diameter with walls 1.25 cm thick. 

The high voltage applied to the Kerr cell was 
obtained by rectifying the current from a five 
kilowatt x-ray transformer. The high voltage 
circuit is essentially the same as that of Stevenson 
and Beams." The voltage applied to the plates 
of the Kerr cell was measured by a Wulf elec- 
trostatic voltmeter capable of measuring up to 
50,000 volts with an accuracy of 0.5 percent. The 
period of the quartz fiber of the voltmeter is such 
that it will detect a 60-cycle ripple in the voltage. 
Under the conditions of this experiment the 
ripple was too small to be detected. 

The gas pressures were measured to 0.1 
atmosphere by a carefully calibrated bourdon 
spring gauge used in conjunction with a dead 
weight piston gauge. Amagat’s™ data on CO, 
were used to turn pressures into densities. 

The gases used were from commercial cylinders 
but tests on CO, showed it to be better than 99.5 
percent pure. The purity of the hydrogen and 
nitrogen is not known. The gases were dried by 
allowing them to pass through P.O; before being 
used. 


III. EXPERIMENTAL PROCEDURE AND RESULTS 


The variation of the Kerr effect with density 
was determined without the rotating sector and 
was standardized by setting the analyzer slightly 
off the parallel position. In this way about 0.2 
percent of the light through the Kerr cell goes 
along with the Kerr beam to the photoelectric 
cell. This light measured with the field off 
furnishes a complete check on the line-up of the 
optical system as well as a check on the sensitivity 
of the amplifier. 

The variation of B with density and the abso- 
lute values were measured at a temperature of 
34.6°C held constant by placing the cell in an 
electrically heated water bath, the temperature 
of which was regulated by a mercury-in-glass 


18 Stevenson and Beams, Phys. Rev. 38, 133 (1931). 
'* Amagat, Comptes Rendus, 663 (1885). 
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Fic. 3. Variation of Kerr constant with density, 


thermostat. It is estimated that the temperature 
varied less than 0.01°C. 

Fig. 3 shows the variation of the Kerr effect 
with density in CO, compared with the variation 
as computed from Eq. (3). The constant of 
proportionality in (3) was so chosen as to make 
the theoretical curve agree with the experimental 
at one point other than at the origin. The squares 
represent points computed by Eq. (3). Values for 
e were taken from the data of Uhlig, Keyes and 
Kirkwood" and n from Phillips'® measures on the 
index of refraction for 4358A. The average 
deviation of the experimental points from the 
curve drawn from theory is less than 1 percent if 
one omits the first measure for which the Kerr 
effect is small. 

Fig. 4 shows the variation of the Kerr constant 
per molecule with density in CO,. Again the 
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Fic. 4. Variation of Kerr constant per molecule 
with density. 


'® Uhlig, Keyes and Kirkwood, J. Chem. Phys. 1, 155 


(1933). 
* Phillips, Proc. Roy. Soc. A97, 226 (1920). 
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curve drawn is a theoretical one computed from 
Eq. (3). The Kerr constant per molecule increases 
with density as predicted by theory. 

Table I gives measurements on the variation 
of the Kerr constant with temperature in CO,. 
BT is constant within the limits of error of the 
measurements, estimated to be approximately 0.5 
percent. 


TABLE I. Variation of Kerr constant with temperature. 


Temperature (absolute) 307.6 280.2 
Density, g/cm*® 0.1035 0.1035 
(Intensity of Kerr beam)! 1.36 1.50 
T (Intensity of Kerr beam)! or BT 418. 420. 


The absolute value of B for CO, was measured 
at 34.6°C for a density of 0.1775 g/cm* and 
wave-length 4550A. This value of B reduced to a 
density corresponding to a pressure of one 
atmosphere is 0.3410-'® and is probably not 
more than 2 percent in error. The most uncertain 
quantity is the wave-length of light used. A Zeiss 
C filter gave a narrow band of wave-lengths, the 
effective wave-length of which depends on the 
character of the source and the variation of the 
sensitivity of the photo-cell with wave-length. 
With a monochromator an intensity wave-length 
curve was obtained from which the effective 
wave-length was determined to be 4550A. 

For comparison with the previously measured 
value of B for COs, it is necessary to reduce the 
above value to wave-length 5890A and 17.5°C. 
From the data of Phillips'* on and from Eq. (4) 
the reduction gives 0.28 The theoretical 
value of B computed by (4) is 0.2810~-"*. The 
experimental value obtained by Szivessy® is 
0.24 10-'*. The agreement here between experi- 
ment and theory shows the essential correctness 
of the theory. From the variation in the values of 
the light scattering coefficient A as determined by 
different observers'’ it would appear that for non- 
polar gases the light-scattering coefficient might 
be determined more accurately from measures of 
the Kerr constant than from direct measures. 

The Kerr constant B was measured for N¢ in 
terms of that for CO, since the value of B for N2 


17 Wien-Harms, Handbuch der Exp. Physik, Vol. XIX. 
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is so small as to make the use of the rotating 
sector impractical. The Kerr constant for Ne» is 
one-sixth that of CO., or B(Ne) is 0.057 
for 4550A and 34.6°C. The probable error of this 
value is 3 percent. Formula (4) gives B(N2) 
=0.05010-'°. This measure was made at a 
density of 0.0915 g/cm’. 

An attempt was made to measure the Kerr 
constant of He. With a pressure of 82 atmospheres 
and a field of 70,000 volts per cm the effect was 
too small to be detected. This means that the 
Kerr constant of He is less than 0.01107, 


BRUCE 


which is what one would expect, the theoretical 
value given by (4) being 0.005 107-"”. 
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Errata.—Variation of the Principal Magnetic Susceptibilities of 
Certain Paramagnetic Crystals with Temperature 


(Phys. Rev. 41, 818, 1932) 


Dr. K. S. Krishnan of Dacca University, India, has recently pointed out an 
error in the tabulation of certain values in the above paper. In Table II, page 
826, the signs of the numbers in the last two rows should all be positive instead 
of negative. In Table V, page 830, under the column headed @, the last four 
figures should read 


79 instead of —79 
—99 —74 
102 —102 
—74 —99 


B. W. BARTLETT 
Bowdoin College, Brunswick, Maine 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important 
discoveries in physics may be secured by addressing 
them to this department. Closing dates for this 
department are, for the first issue of the month, the 


twentieth of the preceding month; for the second 
issue, the fifth of the month. The Board of Editors 
does not hold itself responsible for the opinions ex- 
pressed by the correspondents. 


The Ionization of the Noble Gases by Positive Alkali Ions.—A Correction 


Dr. Otto Beeck has recently informed me of an error in 
his calibration of the McLeod gauge used in the work on 
the ionization of krypton and xenon by positive alkali ions 
on which we collaborated;'!: * and also in the work on the 
ionization of neon and argon by myself.’ He suggests that 
I publish the corrections in this country. The error was 
caused by his mistaking a piece of coordinate paper with 
20 divisions per inch (used for the scale of the gauge) for 
one divided in millimeters. Thus an error of nearly 38 per- 
cent was made in the pressure readings. The values for the 
efficiencies of ionization (the number of electrons ejected 
from the gas per initial positive ion per cm path per mm 
pressure of mercury at 0°C) in krypton and xenon! ? and 
those in neon and argon*® must be corrected. 


NEON . 
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The corrected curves (see Fig. 1) for the efficiencies of 
ionization of the noble gases by the alkali ions are plotted 
with those of Beeck* for neon and argon; and those of 
Sutton,’ Frische,* and Nordmeyer’ for K* in argon. The 
unmarked curves in neon and argon are those of Beeck* 
measured in Germany with accelerating potentials of 0-600 
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volts. The dotted curves are mine* (corrected) measured 
with an apparatus of design different from Beeck's with 
accelerating potentials of 600-2000 volts. The curves for 
krypton and xenon are those of Beeck and myself? (cor- 
rected). The efficiency of ionization is plotted as the or- 
dinate while the accelerating potential in volts is plotted as 
the abscissa. The noteworthy feature of these comparative 
curves is the striking fit between the curves of Beeck and 
those of mine taken in neon and argon with apparatus of 
different design. The importance in this fit lies in the fact 
that it indicates that there is no measurable wide-angle 
scattering of the alkali ions in neon and argon except for 
the case of K* in argon. My apparatus* was designed to 
record all positive ions entering the ionization chamber, 
while that of Beeck* recorded only those striking the col- 
lector plate. Any bombarding positive ions scattered out- 
side the cone determined by the entrance slit of the ioniza- 
tion chamber and the collector plate would be lost in Beeck’s 
apparatus. Therc fore, if any wide-angle scattering occurred, 
my values for the efficiencies of ionization would have been 
lower than the previous ones. This is seen to be the case 
only with K* in argon. 

The correction of the pressures does not alter my state- 
ment concerning the “probability of ionization’ (the 
chance of an electron being ejected in a collision between a 
bombarding positive ion and a bombarded gas atom, cal- 
culated from the kinetic theory mean free path for the 
positive ions in the noble gases nearest them in the periodic 
table). The values of these ‘‘probabilities,”” computed from 
data taken with the same apparatus': ‘4 with 500 volts as 
the accelerating potential of the positives, are: Na* in 
neon, 0.025; K* in argon, 0.026; Rb* in krypton, 0.024*; 
Cs* in xenon, 0.019*, 

J. CarLisLE Movzon 

Duke University, 

Durham, N. C., 
August 23, 1933. 


* Corrected values, 

10. Beeck and J. C. Mouzon, Ann. d. Physik 11, 737 
(1931). 

2 O. Beeck and J. C. Mouzon, Phys. Rev. 38, 967 (1931). 

3 J. C. Mouzon, Phys. Rev. 41, 605 (1932). 

*O. Beeck, Ann. d. Physik 6, 1001 (1930). 

®°R. M. Sutton, Phys. Rev. 33, 364 (1929). 

®C. A. Frische, Phys. Rev. 43, 160 (1933). 

7M. Nordmeyer, Ann. d. Physik 16, 717 (1933). 
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The Possibility of Electric Charge 


(1) The difficulty in understanding the existence of cor- 
puscular electric charges has been expressed frequently by 
the statement that an electron would dissolve within a 
short time if there were no forces other than electric forces 
to keep it together. If one puts the question what would 
be the order of magnitude of the time in which the electron 
would dissolve in the absence of all non-electric forces, 
one finds that a continuous distribution of electric charge 
has no possibility of existence at all. The question what 
will an accumulation of charge do under its own forces 
has no answer whatsoever because a continuous charge 
distribution “doing’’ anything would contradict either 
Maxwell's equations or the rules of conservation of energy 
and momentum and the inertia of energy. Furthermore 
there is no possibility of inventing any set of mechanical 
forces that would make it possible for a charge distribution 
to exist according to these rules, except in symmetrical 
cases such as a single infinite plane distribution or a single 
spherical distribution of charge, but barring even the 
possibility of two charged particles. The conclusion to be 
drawn from this is twofold. Either one gives up the exist- 
ence of continuous charge altogether and supposes the 
electron to be a mathematical point, which would contra- 
dict the inertia of energy, or one gives up the definition of 
the charge and flux of charge which is contained in Max- 
well’s equations: p=div E, pv/ec=curl H—E/c. 

(2) Maxwell's equations are: AivM=P and AivM*=0, 
where M is the six-vector of the field, M* its dual, and P the 
four-vector of the flux of charge with the components pv/c 
and tp. An identical transformation leads to the equation 


—AivT*=[P, M]=F, 


where the tensor 7* is defined by Maxwell's tensions p, 
Poynting’s vector S and the energy density w in the form 


Ps Poy Pas 1S,/c 
Tt= Pyz Py tSy/¢ 
Pros Poy Pos 1S,/¢ 


(i/c)Sz (i/c)S, (i/c)S, —w 


and [P, M] has the components p(E+[(v/c)1]) = Lorentz 
force and (i/c)(v- F) = rate of work. The conservation rules 
of momentum and energy require, if F0, that there are 
non-electric effects with a tensor 7™ so that 


Ato(T*+T™) =0, or —AwT*= F=Awl™. 


Since F #0 means P #0, we would conflict with the conser- 
vation rules if we put 7"=0. So a charge does not “‘dis- 
solve" within a short time under its own electric forces but 
it cannot exist in any stage of stability or dissolution at all. 

(3) The general solution 7 of the equation Aiv7™ = 


—AivT* is 


where 7* is any solution of the homogeneous equation 
AwT*=0, belonging to a Maxwell field without charge. 
(a) The special case 7*=—T* has no interest since it 
would mean that the electrodynamic effects are just can- 
celled by some anti-electrodynamic effects at every point. 
(b) We try to take 7 as a ‘“‘mechanical” effect. Mechan- 
ical quantities may be defined as bound to the substance of 
the charge and moving with the variable velocity v in such 
a way that there is no flux of energy S® and hence no mo- 
mentum g® with respect to a local system in which the 
charge rests: S°=c*g*=0, The 10 components P™, S*, w™ of 
T™ in a system in which the charge moves with v can be 
expressed by the seven rest quantities 


by help of three components of v. In a given electromag- 
netic field the charges and their velocities v, 9, ¥, are de- 
fined in every world point by P=AivM. To bring a given 
electromagnetic field into agreement with the conservation 
rules by help of mechanical forces means to satisfy the four 
equations for 7™ 
AtvT™ = —AwT* (four components), 

where the right-hand side is given and the 10 components 
P», Sm, w" of T™ on the left are dependent on the 7 rest- 
functions p® and w®, which are still free. However, if we 
really want to consider 7 as mechanical forces bound to 
the substance of the charge we must assume that the rest- 
tensions P® in the charge will be some definite physical 
functions of the local value of the gradient of v and of 
higher derivatives of v. Since v is already given by the 
electrodynamic field, the six px" will not be free any more 
in their local values even if we leave the general form of 
their dependence on grad v still free. Hence we have only 
seven less six that is one function free to satisfy the four 
equations AivT™ = —AtvT* by help of mechanical forces 7, 

We see, neither can a continuous charge distribution 
exist by itself nor can it be enabled to exist by help of 
mechanical forces except for the case where the electric 
forces themselves represent a kind of mechanical system 
with vanishing flux of energy and barring the slightest 
motion of the charge in space and time, if there is not 
perfect plane or spherical symmetry, excluding the case of 
more than one electric particle. 

A. Lannt 
Ohio State University, 
August 27, 1933. 
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The Emission and Fluorescence Spectra of Sulphur Dioxide 


In a paper in this issue, Chow has discussed the emission 
and absorption spectra of SO». In the Zeitschrift fiir Physik* 
which came to our laboratory just after Chow's paper was 
sent in, Lotmar has presented his results on the fluorescence 
spectrum in the same region. 

Lotmar's fluorescence was excited by the Zn line at 
2100A (47,600cm™) which falls on the SO, absorption 
band having its violet edge at 47,623 according to Chow. 
This corresponds to excitation to the B, level in Chow's 
scheme of analysis. Since only part of any band will appear 
in fluorescence there may be a considerable interval be- 
tween the center of the fluorescence band which is meas- 
ured by Lotmar and the violet edge of the same band 
measured by Chow in emission. Furthermore this interval 
may vary for different bands, especially for different types 
of vibrational transition. Consequently a direct comparison 
of the wave numbers in Chow's list and Lotmar'’s list has 
little point. However, comparing the 47,600 value given by 
Lotmar for the fluorescence band corresponding to the 
47,623 (or 47,631 in emission) band of Chow, we may guess 
that Lotmar's wave numbers should run about 25 cm™ 
below Chow’s. 

Making this correction I have compared Lotmar'’s series 
with Chow's data in detail and find a fair agreement. Of the 
three 1370 (v; of Chow, », of Lotmar) progressions reported 
by Lotmar (a, 8, y-series) none had been reported com- 
pletely by Chow, but the first two and fifth and sixth of 
the a-series and the first two of the y-series are in his list. 
The £-series is not. 

Of the 1150 (», of Chow, », of Lotmar) progressions, the 
A, B and C series are almost completely represented by 
bands in Chow's list but by bands usually fitting better 
elsewhere in his analysis. The D, E, F and G series receive 
no support from Chow's work. 

Of the 520 (v2 of Chow, »; of Lotmar) progressions, most 
of the a-series has already been recorded by Chow, while 
the c-series has many agreements in spite of four gaps in 
the list of thirteen bands. The 4, d, e, f and g series receive 
little or no confirmation but the six members of the A-series 
all appear in Chow's list. 

A more significant and more striking comparison is that 


between the numerical values for the vibrational levels in 
the normal state as derived from Chow's and from Lotmar's 
data. They are given in Table I. 


TABLE I, Vibration levels in SOz. vy = 1150; vg=520; v3= 1370, 


Level Chow Lotmar| Level Chow _Lotmar 
1150 1150 521 520 
2, 2300 2300 2v2 1043 1040 
3450 3450 1553 1560 
4, 4600 4600 2069 2075 
5750 2578 2590 
6900 6v2 3086 31 10 

etc. 7v2 3590 3630 
4093 
etc. 
lv; 1354? 1370 | n+. 1670 1670 
2v3 2676? 2720 | 2205 2180 
3966? 4060 | v5 2502 2505 
5385 
etc 


It is clear that the only serious discrepancy is in the », 
progression where Chow's values were admittedly tenta- 
tive. I think the long progressions of » and v2 with very 
slow convergence are well established, particularly the ». 
The v3 progression depends almost entirely on Lotmar's 
work, with some confirmation in the a-series. Perhaps the 
most interesting agreement is the »,+¥; level. This level 
came out of Chow's analysis purely empirically yet is found 
to agree perfectly with Lotmar’s result and with the infra- 
red value of 2499.2 As Lotmar points out coupling between 
these two vibrations is reasonable and I think we may say 
that we know its amount quite accurately for this first 
combination. 

H. D. Smyru 

Palmer Physical Laboratory, 

Princeton, N. J., 
September 15, 1933. 


1 Lotmar, Zeits. f. Physik 83, 765, 785 (1933). 
? Dadieu and Kohlrausch, Phys. Zeits. 33, 167 (1932). 


The Isotopes of Hydrogen by the Magneto-Optic Method. The Existence of H* 


In order to test further the presumption that the minima 
observed in the magneto-optic method as developed by 
Allison! depend upon the isotopes of the positive ions pres- 
ent, we examined water solutions containing approximately 
2 and 4 percent of heavy hydrogen (H?). A solution of HCI 
in ordinary water gives minima at 15.74 and 15.85 (Allison 
units); the latter, which is presumably due to H!, is con- 
siderably stronger. These values agree with those reported 
by Allison in his earlier work in which he concluded the 
existence of H*. In the heavy water solution we find the 
intensity of the 15.74 minimum almost equal to the one at 
15.85 and in addition a third minimum at 15.65. The posi- 
tion of this minimum corresponds to that expected of H’, 
and its intensity is somewhat less than the 15.74 minimum 
(H?) in ordinary water. 


We have checked the HCI results with HBr. The minima 
in the solution in ordinary water are at 12.97 and 13.09 and 
the new minimum in the heavy water solutions is at 12.84. 

The heavy water was supplied to us by Professor G. N. 
Lewis, and he and Dr. Spedding are carrying out further 
spectroscopic investigations of the existence of H?®. 

WeENDELL M. LATIMER 
Herpert A. YOuNG 
Department of Chemistry, University of California, 
Berkeley, California, 
September 25, 1933. 


1 Allison and Murphy, J. A. C. S. 52, 3796 (1930). 
Allison, Ind. Eng. Chem. Anal. Ed. 4, 9 (1932). See Phys. 
Rev. 37, 1178 (1931) for halogen acid values. 
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The Raman Spectrum of Gaseous Carbon Dioxide 


In a recent paper! it was ventured, that of the six new 
lines observed by Langseth and Nielsen,’ probably not more 
than two are members of the Raman spectrum of gaseous 
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CO,. This conclusion was based wholly upon the predicted 
Raman spectrum of the molecule since the corresponding 
energy levels had not as yet been observed in the infrared. 

More recently* Professor E. F. Barker has succeeded in 
determining the positions of three of the five relevant levels 
from direct measurements in the infrared. The discovered 
values are in excellent accord with the calculated ones 
referred to above. 

In the light of these results one may expect to find the 
Raman lines in question in the immediate neighborhood of 
the Raman transitions (AJ =0) implied by the energy level 
diagram given in Fig. 1. 

ArTHUR ADEL 

University of Michigan, 

September 25, 1933. 

' Adel and Dennison, Phys. Rev. 43, 716 (1933). 

* Langseth and Nielsen, Zeits. f. physik. Chemie B19, 
427 (1932). 

’ Barker and Wu, to be published in the near future. 


The Allison Magneto-Optic Effect 


In an investigation of a transmission line and spark gap 
system similar to that used by Professor Fred Allison, a few 
things have been noted which seem of some significance. 
The trolley arrangement has been slightly modified to form 
a real two-wire transmission line between the supply con- 
denser and the terminating coil. In order to eliminate the 
long period oscillatory current from the line, and to obtain 
the high-frequency oscillations occurring in the spark gap 
discharge, the line and the condenser spark gap system 
have been made separate circuits with the line connected 
directly across the spark gap terminals through two small 
series capacitances. The capacitance on the low voltage 
side of the gap is paralleled by a high resistance. In this 
arrangement the only voltage applied to the line is that 
appearing across the gap terminals. 

A preliminary investigation of the coils at the line 
terminals reveals the presence of high frequency oscillations 
ranging from some centimeters to a few meters in wave- 
length. A change in trolley position produces a marked shift 
in energy distribution between these high frequency com- 
ponents. This shift or selective effect is apparently from the 
low frequency to the high frequency side of the band as the 
line length between gap and coils is increased. These 
oscillations are due to the capacity between the electrodes 
of the gap and possibly also to phenomena of the Langmuir- 
Tonks type. In such a discharge one would naturally expect 
a broad band of frequencies. 

It is futile at present to speculate upon the complete 
significance of these oscillations, but it is hoped that some 
correlation between them and the Allison minima can be 


obtained when more data are available. An investigation is 
now in progress in which an attempt will be made to cali- 
brate the trolley system in terms of predominating fre- 
quencies and determine whether or not there is selective 
absorption occurring in the cell. 

It might be worth while to note that an arrangement 
similar to the above was tried out by the writer in Pro- 
fessor Allison's laboratory at Auburn this summer. The 
minima were apparently made considerably sharper, that 
is more sensitive to trolley shift. Several analyses of un- 
knowns were made with this modification and the isotopes 
of lead redetermined. The readings were in excellent agree- 
ment with those of Bishop, Lawrenz and Dollins.' It is 
perhaps needless to state that for this modification the 
nicols must be used in the parallel position and the coil 
fields assistinp. 

I am much indebted to Professor Fréd Allison for the 
privilege of using his apparatus and for his courtesy in 
allowing me to rearrange it in any desired manner, I have 
also profited a great deal by many discussions with Pro- 
fessor J. W. Beams. 

This research is made possible by grants from the 
National Research Council, the Research Corporation, and 
the University of Virginia. 

L. B. SNoppy 

Rouss Physical Laboratory, 

University of Virginia, 
September 30, 1933, 


1 Bishop, Lawrenz and Dollins, Phys. Rev. 43, 43 (1933). 
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Production of Neutrons by High Speed Deutons 


In a recent letter to the Physical Review' we described 
briefly an apperatus for the acceleration of positive ions and 
the production of neutrons by bombarding beryllium with 
helium ions of from 600,000 to 975,000 electron-volts energy. 
Since then we have bombarded targets of lithium chloride 
and beryllium with deutons of energies up to 900,000 
electron-volts, and have obtained in both cases yields of 
neutrons several hundred times as great as the yield 
previously obtained by bombarding beryllium with helium 
ions. 

The same apparatus was used for accelerating the ions 
and for detecting the neutrons as was described in the 
previous letter. A sensitive electroscope, the inside walls of 
which were coated with paraffin, was entirely enclosed in a 
lead cylinder of 5 cm wall thickness to shield it from x-rays. 
The neutrons supposedly penetrate the 5 cm of lead and 
eject recoil hydrogen particles from the paraffin on the 
inner walls of the electroscope. It is the ionization due to 
these recoil particles that gives a measure of the intensity 
of neutrons passing through the chamber. The lead cylinder 
with the electroscope inside was placed in such a position 
that the center of the electroscope was 13 cm from the 
center of the target, and ina direction from it perpendicular 
to the incident ion beam. 

In making the measurements of neutrons produced by 
deuton bombardment, it was necessary to make several 
comparison measurements in order to establish the fact 
that in each case the observed effect was due to neutrons 
rather than to y-rays, and that the disintegration was pro- 
duced by the deutons rather than protons. The hydrogen 
used to produce the deutons contained only 5 percent of 
the heavy isotope, so even in the so-called deuton measure- 
ments, the greater part of the ions was protons. Each 
target bombarded with deutons was later bombarded with 
100 percent protons for comparison. No effect was observed 
with protons, except a slight one in the case of lithium, 
which is discussed below, and it was therefore concluded 
that the large effects were produced by deutons. Measure- 
ments were also made with deutons, with the paraffin 
removed from the electroscope. With both lithium and 
beryllium the raté of discharge without paraffin was less 
than half that with paraffin, indicating that the greater 
part of the ionization observed was due to recoil particles 
ejected from the paraffin by neutrons, rather than to 
+-tays. The sensitivity of the electroscope to 600,000 volt 
x-rays filtered through 5 cm of lead was measured and 
found to be slightly greater without paraffin than with 
paraffin. 

Curves of the efficiency of production of neutrons from 
lithium chloride and from beryllium as a function of 
voltage are shown in Fig. 1. The total ion current during 
all the measurements was 10 microamperes. The back- 
ground, indicated by a straight line, is the rate of discharge 
of the electroscope due to x-rays plus its natural leak, and 
is the same in all cases in which no neutrons or y-rays are 
produced. It has been determined by bombarding a brass 
target with helium ions, a beryllium target with protons, 
and a brass target with protons. 


The lithium chloride target, bombarded with 10 micro. 
amperes of protons at 800,000 volts gave a small effect, 
about twice the background, which was decreased to less 
than half that when the paraffin was removed from the 
electroscope, indicating that the radiation was neutrons, 
It is possible that this is a double process in which the 
a-particles produced by proton bombardment in turn dis. 
integrate lithium with the production of neutrons. It has 
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already been shown by Curie and Joliot? that a-particles of 
long range are capable of producing neutrons from lithium, 
and the effect here observed seems to be of the right order 
of magnitude to be accounted for in that way. However, 
the effect is to be investigated further. 

The simplest hypotheses as to the nature of the two 
disintegrations with deutons are: 


n!, 


Assuming tentatively that the mass of the neutron is 
1.0065, and taking as values of the other atomic masses 
referred to 


Li? 7.0146 
Be® 9.0155 


B' 10.0137 H? 2.0136 


He* 4.0022 


we find that in the case of the disintegration of beryllium 
the neutron comes off with a kinetic energy of about 
9X 10* electron-volts. In the case of lithium, the kinetic 
energy of the neutron plus that of the two @-particles comes 
out to be about 16 ¥ 10* electron-volts. It is probable that 
these calculations will be greatly modified by subsequent 
experimental data, but they serve as a rough approximation 


! Crane, Lauritsen and Soltan, Phys. Rev. 44, 514 (1933). 
2 Curie and Joliot, J. de Physique 4, 278 (1933). 
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to show that the neutrons produced by this method have a 
rather high kinetic energy. 

It is of interest to make a comparison of the intensity of 
neutrons produced by bombarding beryllium or lithium 
with deutons with that obtained by bombarding beryllium 
with a-particles from a strong polonium source. One divi- 
sion in our electroscope corresponds to about 70 recoil 
hydrogen particles. At our maximum deuton current 
(30 microamperes) and at 900,000 volts, about 1000 recoil 
particles per minute would be obtained, which is of the 
order of one hundred times the intensity obtained by 
means of the strongest polonium a@-particle sources now in 
use. It is therefore possible, by the use of deutons, to obtain 
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an efficiency in the artificial production of neutrons which 
is comparable with the efficiency of production of a-parti- 
cles by proton bombardment. 

We are greatly indebted to Professor G. N. Lewis, who 
kindly furnished the heavy hydrogen used in these ex- 
periments. 

H. R. 
C. C. LAuRITSEN 
A. SOLTAN 
Kellogg Radiation Laboratory, 
California Institute of Technology, 
September 30, 1933. 
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